24小时热门版块排行榜    

查看: 748  |  回复: 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

金虫 (小有名气)

【答案】应助回帖

如果找不到加速的方法,可考虑分块计算或试用使用其它软件。
4楼2022-09-13 10:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0857 资源与环境 285分 +3 未名考生 2026-03-09 3/150 2026-03-10 23:25 by rock_sdust
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 一志愿211,0860总分286食品工程专业求调剂 +8 林林Winnie 2026-03-05 8/400 2026-03-10 18:44 by 清风月
[考研] 2026考研求调剂-材料类-本科211一志愿985-初试301分 +10 虫友233 2026-03-07 10/500 2026-03-10 17:10 by Demonsssss
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +7 小小番的茄 2026-03-09 7/350 2026-03-10 16:48 by ztnimte
[考研] 0703求调剂 +3 jtyq001 2026-03-10 3/150 2026-03-10 16:45 by ztnimte
[考研] 08工科 +5 li李乐成 2026-03-06 5/250 2026-03-10 14:51 by 在风落中
[考研] 一志愿清华深研院材料专硕294分,专业课111分,本科中南大学材料,有六级,有工作经验 +4 H14528 2026-03-04 4/200 2026-03-10 11:03 by dj19940201
[考研] 293一志愿华东理工 0817化学工程与技术 调剂 +5 fjj0912 2026-03-07 5/250 2026-03-09 09:13 by 30660438
[考研] 297求调剂 +3 胡达灵 2026-03-05 5/250 2026-03-08 23:13 by 清风月
[考研] 求调剂 +3 拾柒12。 2026-03-08 3/150 2026-03-08 17:33 by bingxueer79
[考研] 化学工程求调剂 +12 化工人999 2026-03-04 12/600 2026-03-08 11:17 by 我的船我的海
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 求调剂 +3 snw石 2026-03-07 3/150 2026-03-07 20:01 by houyaoxu
[考研] 诚求调剂,323分 有专利/科研/实习经历 +6 橙子cyx 2026-03-06 8/400 2026-03-07 19:06 by 新篇章DFSS
[考研] 304求调剂 +4 52hz~~ 2026-03-05 5/250 2026-03-07 15:47 by lature00
[考研] 085701环境工程 求调剂 +7 xiiiia 2026-03-04 7/350 2026-03-05 09:49 by 热情沙漠
[考研] 293求调剂 +3 是乐渝哇 2026-03-04 3/150 2026-03-04 23:12 by wutongshun
[考研] 307求调剂 +4 xuyuan111 2026-03-04 5/250 2026-03-04 20:39 by kakakapanpan
[考研] 材料专硕346求调剂 +3 旺一下 2026-03-04 3/150 2026-03-04 16:26 by sslc1985
信息提示
请填处理意见