24小时热门版块排行榜    

查看: 745  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 265求调剂 +8 小木虫085600 2026-03-09 11/550 2026-03-10 19:44 by 绿色低碳化学
[考研] 一志愿211,0860总分286食品工程专业求调剂 +8 林林Winnie 2026-03-05 8/400 2026-03-10 18:44 by 清风月
[考研] 调剂 +3 呵唔哦豁 2026-03-10 3/150 2026-03-10 17:43 by Nuyoah678
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +3 18266118446 2026-03-09 3/150 2026-03-10 15:50 by houyaoxu
[考研] 0860求调剂(272分) +3 lllllcsjsj 2026-03-05 4/200 2026-03-10 15:29 by circleffyy
[考研] 哈工大材料324求调剂 +5 闫旭东 2026-03-10 7/350 2026-03-10 14:41 by ruiyingmiao
[考研] 310 070300化学求调剂 +4 扑风铃的猫 2026-03-08 5/250 2026-03-09 19:46 by 30660438
[基金申请] 第二个面上项目的获批几率如何? +3 wawm9296 2026-03-04 5/250 2026-03-09 13:00 by xujie4132
[考研] 中科大材料299求调剂 +10 DAIjiayo 2026-03-05 16/800 2026-03-09 10:45 by 斩魂滴兔子!
[考博] 26申博自荐 +9 烷基八氮- 2026-03-06 11/550 2026-03-08 23:20 by gaobiao
[考研] 303求调剂 +8 forgman95 2026-03-05 10/500 2026-03-08 12:41 by 蓝莓都给你
[考研] 301求调剂 +5 一二LV 2026-03-07 5/250 2026-03-07 22:20 by 18137688336
[考研] 第一志愿上海大学,专业化学工程与技术,总分288,求调剂 +3 1829197082 2026-03-07 3/150 2026-03-07 19:14 by houyaoxu
[考研] 276求调剂 本科天津工业大学 一志愿北京化工大学0855 +4 路lyh123 2026-03-06 6/300 2026-03-07 11:37 by shdliugang
[考研] 材料工程调剂 +6 zbcm_zbcm 2026-03-05 6/300 2026-03-06 15:36 by @飒飒飒飒
[考研] 求调剂 +5 danyyyy 2026-03-04 5/250 2026-03-06 08:16 by Leeding1356
[考研] 材料与化工304求调剂 +7 邱gl 2026-03-05 10/500 2026-03-06 01:11 by 卷柏卷柏
[考研] 325求调剂 +5 学家科 2026-03-04 5/250 2026-03-04 14:04 by kakakapanpan
[考研] 一志愿西交大材料学硕 346 求调剂 +3 zju51 2026-03-04 3/150 2026-03-04 13:01 by zhukairuo
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
信息提示
请填处理意见