偏微分方程利用c++表达式模板与共轭梯度法实现

上传者: xishan555 | 上传时间: 2022-06-03 05:46:40 | 文件大小: 28KB | 文件类型: RAR
Vector.h 是向量类,包含生成向量及各种操作符重载 Template.h 是各种表达式模板的集合,包含必要的向量加减乘法操作符重载。 Matrix.h 是AX=b中,关键A数组生成的类,这里我用了数组压缩技术,即把带状数组A压缩,使程序执行更有效率。 cgexpr.cpp是Main函数,包括使用三种时间差分即 Explicit,Implicit 和 CrankNicolson。 运行程序时需要在Command line里按如下格式输入10个指令: cgexpr hx hy tend tao a iterations eps residual.txt result.txt 其中cgexpr是主函数文件名,hx,hy,是有限差分对应的x,y大小,tend是时间长度,tao是时间差分对应的ht大小,a是使用哪种时间差分格式:0是Explicit,0.5是CrankNicolson,1是implicit.iterations 是一个时间段里循环的最大次数,eps是你设定的误差。residual.txt 和 result.txt 分别是误差和最后结果输出。 这个程序实现的偏微分方程是: @u/@t= (delta)u 你可以根据你需要计算的偏微分方程,修改Matrix.h中对应的m,n,t的表达式即可。具体表达式需根据你的方程推出。

文件下载

资源详情

[{"title":"( 9 个子文件 28KB ) 偏微分方程利用c++表达式模板与共轭梯度法实现","children":[{"title":"Ass4","children":[{"title":"CG&PDE","children":[{"title":"cgexpr.cpp <span style='color:#111;'> 4.88KB </span>","children":null,"spread":false},{"title":"cgexpr <span style='color:#111;'> 35.50KB </span>","children":null,"spread":false},{"title":"cgexpr.o <span style='color:#111;'> 30.38KB </span>","children":null,"spread":false},{"title":"Template.h <span style='color:#111;'> 3.26KB </span>","children":null,"spread":false},{"title":"readme.txt <span style='color:#111;'> 856B </span>","children":null,"spread":false},{"title":"Vector.h <span style='color:#111;'> 1.57KB </span>","children":null,"spread":false},{"title":"Makefile <span style='color:#111;'> 289B </span>","children":null,"spread":false},{"title":"Matrix.h <span style='color:#111;'> 1.84KB </span>","children":null,"spread":false},{"title":"cgexpr.cpp~ <span style='color:#111;'> 4.88KB </span>","children":null,"spread":false}],"spread":true}],"spread":true}],"spread":true}]

评论信息

  • zys134679 :
    有点小贵,但是还是值得的
    2015-10-12

免责申明

【只为小站】的资源来自网友分享,仅供学习研究,请务必在下载后24小时内给予删除,不得用于其他任何用途,否则后果自负。基于互联网的特殊性,【只为小站】 无法对用户传输的作品、信息、内容的权属或合法性、合规性、真实性、科学性、完整权、有效性等进行实质审查;无论 【只为小站】 经营者是否已进行审查,用户均应自行承担因其传输的作品、信息、内容而可能或已经产生的侵权或权属纠纷等法律责任。
本站所有资源不代表本站的观点或立场,基于网友分享,根据中国法律《信息网络传播权保护条例》第二十二条之规定,若资源存在侵权或相关问题请联系本站客服人员,zhiweidada#qq.com,请把#换成@,本站将给予最大的支持与配合,做到及时反馈和处理。关于更多版权及免责申明参见 版权及免责申明