xishan555的专栏

具体细节请参加我的共轭梯度法求解偏微分方程无并行版本。这个程序是建立在上一个基础上的升级。不同的是这里我用的是五点法而不是压缩的矩阵A。

程序用红黑高斯-塞德尔迭代法及有限差分实现偏微分方程： −(delta)u + k^2*u = f 设定有限差分hx,hy,误差大小，迭代次数，即可计算各个节点的数值。 如需计算其他偏微分方程，只需改动迭代部分的表达式即可。 程序内有输出部分，可以将结果直接按文件输出。如不需要，删除相关语句即可。

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的表达式即可。具体表达式需根据你的方程推出。

Strassen 算法 C++ 实现 任意矩阵相乘 用Command line 输入编制好的两个矩阵，输出相乘的结果矩阵。如果需要手动输入，删除程序里的相关语句，添加输入命令即可。