CVODE中雅可比的影响

Feisy

已于 2022-03-14 04:49:12 修改

阅读量434

点赞数

分类专栏： SUNDIAL 文章标签： sundials

于 2022-03-14 04:48:30 首次发布

本文链接：https://blog.csdn.net/feisy/article/details/123470071

版权

SUNDIAL 专栏收录该内容

27 篇文章 4 订阅

订阅专栏

基于cvRoberts_dns那个例子进行测试

发现：没有提供Jac的情况下速度更快一点。

这个是不是跟模型有关，还是都是这样的表现？

有提供Jac:

At t = 2.6391e-01      y =  9.899653e-01    3.470564e-05    1.000000e-02
    rootsfound[] =   0   1
At t = 4.0000e-01      y =  9.851641e-01    3.386242e-05    1.480205e-02
At t = 4.0000e+00      y =  9.055097e-01    2.240338e-05    9.446793e-02
At t = 4.0000e+01      y =  7.158017e-01    9.185037e-06    2.841892e-01
At t = 4.0000e+02      y =  4.505360e-01    3.223271e-06    5.494608e-01
At t = 4.0000e+03      y =  1.832299e-01    8.944378e-07    8.167692e-01
At t = 4.0000e+04      y =  3.898902e-02    1.622006e-07    9.610108e-01
At t = 4.0000e+05      y =  4.936383e-03    1.984224e-08    9.950636e-01
At t = 4.0000e+06      y =  5.168093e-04    2.068293e-09    9.994832e-01
At t = 2.0790e+07      y =  1.000000e-04    4.000397e-10    9.999000e-01
    rootsfound[] =  -1   0
At t = 4.0000e+07      y =  5.202440e-05    2.081083e-10    9.999480e-01
At t = 4.0000e+08      y =  5.201061e-06    2.080435e-11    9.999948e-01
At t = 4.0000e+09      y =  5.258603e-07    2.103442e-12    9.999995e-01
At t = 4.0000e+10      y =  6.934511e-08    2.773804e-13    9.999999e-01

Final Statistics:
nst = 542    nfe  = 754    nsetups = 107    nfeLS = 0      nje = 11
nni = 751    ncfn = 0      netf = 22     nge = 570

use time = 16

没有提供jac:

At t = 2.6391e-01      y =  9.899653e-01    3.470564e-05    1.000000e-02
    rootsfound[] =   0   1
At t = 4.0000e-01      y =  9.851641e-01    3.386242e-05    1.480205e-02
At t = 4.0000e+00      y =  9.055097e-01    2.240338e-05    9.446793e-02
At t = 4.0000e+01      y =  7.158016e-01    9.185043e-06    2.841892e-01
At t = 4.0000e+02      y =  4.505209e-01    3.222826e-06    5.494759e-01
At t = 4.0000e+03      y =  1.832217e-01    8.943516e-07    8.167774e-01
At t = 4.0000e+04      y =  3.898091e-02    1.621669e-07    9.610189e-01
At t = 4.0000e+05      y =  4.936971e-03    1.984450e-08    9.950630e-01
At t = 4.0000e+06      y =  5.170103e-04    2.069098e-09    9.994830e-01
At t = 2.0812e+07      y =  1.000000e-04    4.000395e-10    9.999000e-01
    rootsfound[] =  -1   0
At t = 4.0000e+07      y =  5.204927e-05    2.082078e-10    9.999480e-01
At t = 4.0000e+08      y =  5.184946e-06    2.073989e-11    9.999948e-01
At t = 4.0000e+09      y =  5.246212e-07    2.098486e-12    9.999995e-01
At t = 4.0000e+10      y =  6.043455e-08    2.417382e-13    9.999999e-01

Final Statistics:
nst = 522    nfe  = 713    nsetups = 104    nfeLS = 36     nje = 12
nni = 710    ncfn = 0      netf = 18     nge = 558

use time = 9