24小时热门版块排行榜    

查看: 744  |  回复: 3
【悬赏金币】回答本帖问题,作者chenye1717将赠送您 18 个金币

chenye1717

新虫 (初入文坛)

[求助] 【请教/求助】Mathematica中如何加速这个常规的三维For Loop语句? 已有1人参与

mathematica中如何加速这个常规的三维for loop语句?

for[i = 1, i <= 1000, i++,
for[j = 1, j <= 1000, j++,
for[k = 1, k <= 1000, k++,
if[x[]^2 + y[[j]]^2 + z[[k]]^2 >= r^2,
funval[[i, j, k]]=
g[x[], y[[j]], z[[k]]], funval[[i, j, k]]= 0.0]]];

我的目标是将有一个在三维空间内(在球体外部,x[]^2 + y[[j]]^2 + z[[k]]^2 >= r^2)构造table (名为funval)存储g函数的离散值。因为g函数比较复杂,所以导致计算速度很慢。
尝试过outer[g, x[], y[[j]], z[[k]]],效果不好,而且好像不能加入“位于球体外”这个限制条件。
请问有加速的方法吗?谢谢大家。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyh700600

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
if[x[]^2 + y[[j]]^2 + z[[k]]^2 >= r^2,
g[x[], y[[j]], z[[k]]

将:
x[]
改为:
x[]
2楼2022-09-10 09:47:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyh700600

金虫 (小有名气)

【答案】应助回帖

引用回帖:
2楼: Originally posted by zyh700600 at 2022-09-10 09:47:41
if[x[]^2 + y[]^2 + z[]^2 &gt;= r^2,
g[x[], y[], z[]

将:
x[]
改为:
x[]

CODE:
if[x[]^2 + y[[j]]^2 + z[[k]]^2 >= r^2,
g[x[], y[[j]], z[[k]]

将:
x[]
改为:
x[[i]]

3楼2022-09-10 09:49:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyh700600

金虫 (小有名气)

【答案】应助回帖

如果找不到加速的方法,可考虑分块计算或试用使用其它软件。
4楼2022-09-13 10:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenye1717 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
[考研] 求调剂 一志愿苏州大学,0856化工323分 | 本科应化 | 有专利/竞赛/科研助手经历 | +7 橙子cyx 2026-03-06 9/450 2026-03-10 12:03 by kiokin
[考研] 调剂的同学,走过路过,不要错过 +11 likeihood 2026-03-06 16/800 2026-03-10 11:38 by Djdjj12
[考博] 2026博士申请 +6 起泡酒 2026-03-08 6/300 2026-03-10 10:26 by 何由彻
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
[考研] 295复试调剂 +4 简木ChuFront 2026-03-09 4/200 2026-03-09 21:49 by wangjihu
[考研] 材料工程085601调剂求老师收留 +5 强木木木 2026-03-07 5/250 2026-03-09 20:09 by yuningshan
[考研] 307求调剂 +3 辛仁豆腐 2026-03-08 5/250 2026-03-09 14:09 by macy2011
[考研] 070300化学求调剂 +5 扑风铃的猫 2026-03-08 10/500 2026-03-09 13:58 by macy2011
[考研] 086000生物与医药319分求调剂 +4 Tolkien 2026-03-07 4/200 2026-03-09 13:49 by macy2011
[考研] 083000环境科学与工程调剂 +5 加油呀fxy 2026-03-07 6/300 2026-03-08 21:07 by sdgdsdgd
[考研] 0703化学求调剂学校 +3 一志愿京区211 2026-03-04 3/150 2026-03-07 21:35 by yinhuanshun
[考研] 一志愿211 085600 280数二英二求调剂 +3 月山斜 2026-03-06 3/150 2026-03-07 18:56 by houyaoxu
[考研] 306求调剂 +7 Bahati 2026-03-05 7/350 2026-03-06 22:11 by 星空星月
[考研] 287求调剂 +3 看看我. 2026-03-05 6/300 2026-03-06 10:40 by Iveryant
[考研] 275求调剂 +4 大爆炸难民 2026-03-06 5/250 2026-03-06 09:21 by guoweigw
[考研] 085602 293分求调剂 +3 SivanNano. 2026-03-05 3/150 2026-03-05 19:41 by caszguilin
[考研] 334求调剂 +3 Trying] 2026-03-03 4/200 2026-03-05 16:39 by zhukairuo
[考研] 325求调剂 +5 学家科 2026-03-04 5/250 2026-03-04 14:04 by kakakapanpan
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
信息提示
请填处理意见