由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Computation版 - 再请教Re: 请教大家c++ array
相关主题
如何用变步长的差分法解PDE?科普,谱元法
求助:DGBTRS scientific computing还是Fortran90好阿
请教一个关于PDE数值解方面的问题lapack++ and blitz++
what is second-order upwind?直接用NR的源程序需要注意什么吗?
谁介绍一本差分法解偏微分方程的入门书?在线等:紧急求助Matlab高手(ODE问题)
[转载] 请教PDE有限差分法的CFL ConditionFast solver of linear equation system
关于超大型矩阵的算法help on nonlinear programming solver
寻支持有限元和差分法的画图软件请大虾推荐一个fully implicit solver.
相关话题的讨论汇总
话题: array话题: mtl话题: 301话题: using话题: sor
进入Computation版参与讨论
1 (共1页)
x*****u
发帖数: 3419
1
首先谢谢大家的回帖,受益良多。
在下还有一件极幼稚的事情不明,还望大家帮忙解释解释。
不知道大家用过MTL(The Matrix Template
Library)没有。当我使用MTL里面的Matrix时,发现运行速度降低。于是我做了一个比较

发现MTL array的执行效率,低于c++ array。对于我使用的这个算例,竟有四倍之差!
*************************
the dimension is 301*301
using c++ dynamic array...
using SOR ...
the relax factor = 1.8
time in the solver(sec): 44.91
interations 11669 ~~~~
Norm of error < 0.000159678
*****************
the dimension is 301*301
using MTL array...
using SOR ...
the relax factor = 1.8
time in the solver(sec):
c******m
发帖数: 599
2
这个mtl的矩阵是parse的么?
对于稀疏矩阵,特别是差分法生成的带状矩阵,用完整的矩阵进行运算是绝对不可取的

【在 x*****u 的大作中提到】
: 首先谢谢大家的回帖,受益良多。
: 在下还有一件极幼稚的事情不明,还望大家帮忙解释解释。
: 不知道大家用过MTL(The Matrix Template
: Library)没有。当我使用MTL里面的Matrix时,发现运行速度降低。于是我做了一个比较
: ,
: 发现MTL array的执行效率,低于c++ array。对于我使用的这个算例,竟有四倍之差!
: *************************
: the dimension is 301*301
: using c++ dynamic array...
: using SOR ...

1 (共1页)
进入Computation版参与讨论
相关主题
请大虾推荐一个fully implicit solver.谁介绍一本差分法解偏微分方程的入门书?
再问个excel solver &VBA问题[转载] 请教PDE有限差分法的CFL Condition
[求助]ODE数值解不converge关于超大型矩阵的算法
Help! 有谁用过LSODE/LSODA?寻支持有限元和差分法的画图软件
如何用变步长的差分法解PDE?科普,谱元法
求助:DGBTRS scientific computing还是Fortran90好阿
请教一个关于PDE数值解方面的问题lapack++ and blitz++
what is second-order upwind?直接用NR的源程序需要注意什么吗?
相关话题的讨论汇总
话题: array话题: mtl话题: 301话题: using话题: sor