我有一个代码，可以根据笛卡尔空间中的节点坐标创建一组方程。然后我使用自动sym生成和一系列嵌套的匿名函数创建基于这些的jacobian。我知道这不是执行操作的最有效方法，但我唯一的另一个想法是在运行时编写和实现.m文件。这是我使用的方程式和方法的样本。
[joeuser@myhost:/tmp]$ nvcc -I/opt/cub -o a a.cu
nvcc warning : The 'compute_20', 'sm_20', and 'sm_21' architectures are deprecated, and may be removed in a future release (Use -Wno-deprecated-gpu-targets to suppress warning).
[joeuser@myhost:/tmp]$ ./a
0, 4000
1, 1000
2, 2680
[joeuser@myhost:/tmp]$ nvcc -I/opt/cub -o a a.cu -gencode arch=compute_30,code=sm_30
[joeuser@myhost:/tmp]$ ./a
0, 4000
1, 1000
2, 512
然后我执行优化，根据负载应用程序fe将从@ eq1输出的所有值设置为0。此过程需要200秒，其中大部分用于mupadmex。救命？我需要将时间顺序减少5或6 ......
同样，我的目标是自动创建jacobian，因为变量a b和c将改变，并影响@ eq1中的方程数