24小时热门版块排行榜    

查看: 747  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328,0703考生求调剂,一志愿为东北师范大学 +3 观素律 2026-03-09 4/200 2026-03-10 15:58 by houyaoxu
[考研] 0860求调剂(272分) +3 lllllcsjsj 2026-03-05 4/200 2026-03-10 15:29 by circleffyy
[考研] 一志愿天大化工(085600)调剂总分338 +5 蔡大美女 2026-03-09 5/250 2026-03-10 14:44 by ruiyingmiao
[考研] 求调剂 一志愿苏州大学,0856化工323分 | 本科应化 | 有专利/竞赛/科研助手经历 | +7 橙子cyx 2026-03-06 9/450 2026-03-10 12:03 by kiokin
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
[考研] 求调剂,数一英一274分 +4 小菲会努力 2026-03-08 4/200 2026-03-09 12:40 by 一定上岸哟_
[考研] 337求调剂 +3 睡醒,。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
[考研] 0817化学工程与技术312分求调剂 +7 T123 tt 2026-03-04 7/350 2026-03-08 11:18 by 我的船我的海
[考研] 材料调剂 +7 ounce. 2026-03-05 15/750 2026-03-08 09:39 by wangjihu
[考研] 301求调剂 +5 一二LV 2026-03-07 5/250 2026-03-07 22:20 by 18137688336
[考研] 求调剂 +4 呼呼?~+123456 2026-03-06 4/200 2026-03-06 23:11 by L135790
[考研] 复试调剂 +7 呼呼?~+123456 2026-03-05 10/500 2026-03-06 16:07 by @飒飒飒飒
[考研] 287求调剂 +3 看看我. 2026-03-05 6/300 2026-03-06 10:40 by Iveryant
[考研] 材料考研339求调剂 +3 Karry*^_^* 2026-03-04 3/150 2026-03-06 00:32 by wutongshun
[考研] 化工282求调剂一志愿211 +5 NA0912 2026-03-05 6/300 2026-03-05 20:10 by 2735147993
[考研] 301求调剂 +4 李LJR 2026-03-04 4/200 2026-03-05 13:31 by Iveryant
[考研] 267调剂求助 +5 聪少OZ 2026-03-04 5/250 2026-03-05 09:38 by kakakapanpan
[考研] 320材料一志愿西工大专硕英二数二 有两年光伏工作经验 +5 锤某人 2026-03-04 8/400 2026-03-05 09:24 by oxidpl
[考研] 085601 材料305分求助 +4 泡泡邮件 2026-03-03 6/300 2026-03-04 13:00 by 无际的草原
[考研] 281求调剂 +3 Y?l?h 2026-03-04 3/150 2026-03-04 10:51 by 汪!?!
信息提示
请填处理意见