24小时热门版块排行榜    

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

金虫 (小有名气)

【答案】应助回帖

引用回帖:
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的回帖
查看全部 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

金虫 (小有名气)

【答案】应助回帖

如果找不到加速的方法,可考虑分块计算或试用使用其它软件。
4楼2022-09-13 10:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 0857环境调剂 +4 熠熠_11 2026-03-10 4/200 2026-03-10 22:47 by 星空星月
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +7 小小番的茄 2026-03-09 7/350 2026-03-10 16:48 by ztnimte
[考研] 材料工程专硕调剂自荐信,初试323 +6 李白26 2026-03-07 6/300 2026-03-10 16:16 by 无际的草原
[考研] 材料专硕调剂 +7 慕辰123 2026-03-05 10/500 2026-03-10 14:22 by peike
[考研] 研究生招生 +3 徐海涛11 2026-03-10 6/300 2026-03-10 12:43 by 徐海涛11
[考研] 求调剂 一志愿苏州大学,0856化工323分 | 本科应化 | 有专利/竞赛/科研助手经历 | +7 橙子cyx 2026-03-06 9/450 2026-03-10 12:03 by kiokin
[考研] 环境调剂 +5 晓看天暮看云 2026-03-09 5/250 2026-03-10 10:46 by 星空星月
[考研] 考研一志愿长安大学材料与化工309分请求调剂 +6 dtdxzxx 2026-03-06 8/400 2026-03-10 08:31 by Linda Hu
[考研] 313求调剂 +4 Yyt杨1 2026-03-07 5/250 2026-03-09 13:48 by macy2011
[考研] 求调剂,数一英一274分 +4 小菲会努力 2026-03-08 4/200 2026-03-09 12:40 by 一定上岸哟_
[考研] 材料调剂 +4 xxxcm 2026-03-08 7/350 2026-03-09 10:11 by EBSD
[考研] 337求调剂 +3 睡醒,。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
[考研] 085701环境工程专硕求调剂 +5 幼稚园大班打手 2026-03-08 5/250 2026-03-09 07:07 by houyaoxu
[考研] 一志愿211 化学305分求调剂 +3 0703杨悦305分 2026-03-05 3/150 2026-03-07 21:15 by yinhuanshun
[考研] 2026调剂】考试A区0703化学类323分 诚求接收 +3 卷柏卷柏 2026-03-05 4/200 2026-03-07 17:00 by yuandd_2001
[考研] 一志愿南京大学资源与环境专业283分求调剂 +8 追云?? 2026-03-04 8/400 2026-03-07 09:23 by 斩魂滴兔子!
[考研] 306求调剂 +7 Bahati 2026-03-05 7/350 2026-03-06 22:11 by 星空星月
[考研] 材料考研339求调剂 +3 Karry*^_^* 2026-03-04 3/150 2026-03-06 00:32 by wutongshun
[考研] 325求调剂 +5 学家科 2026-03-04 5/250 2026-03-04 14:04 by kakakapanpan
信息提示
请填处理意见