错误不在sor，在这里：
创新互联-专业网站定制、快速模板网站建设、高性价比潞城网站开发、企业建站全套包干低至880元,成熟完善的模板库,直接使用。一站式潞城网站制作公司更省心,省钱,快速模板网站建设找我们，业务覆盖潞城地区。费用合理售后完善，10多年实体公司更值得信赖。
for(j=0;jm;j++)
{
for(l=0;ln;l++)
{
a[j][l]=rand();
b[j][l]=rand();
c[j][l]=rand();
d[j][l]=rand();
e[j][l]=rand();
f[j][l]=rand();
u[j][l]=0.0;
}
}
怀疑你的dmatrix函数写的有问题，代码写得很乱，不知道在干什么。估计是内存分配有问题。我这里的运行时错误是访问了非法地址0。环境是VS 2005.
我把你的程序稍微修改了一下，可以看得比较清楚错误的地方。
#include stdlib.h
#include stdio.h
//#include "nrutil.h"
#include time.h
#include math.h
#include stddef.h
#define MAXITS 1000
#define NR_END 1
#define FREE_ARG char*
#define EPS 1.0e-5
void nrerror(char error_text[]);
void sor(double **a,double **b,double **c,double **d,double **e,double **f,double **u,int jmax,double rjac)
{
int ipass,j,jsw,l,lsw,n;
double anorm,anormf=0.0,omega=1.0,resid;
printf("enter sor\n");
for(j=2;jjmax;j++)
for(l=2;ljmax;l++)
anormf+=fabs(f[j][l]);
for(n=1;n=MAXITS;n++)
{
anorm=0.0;
jsw=1;
for(ipass=1;ipass=2;ipass++)
{
lsw=jsw;
for(j=2;jjmax;j++)
{
for(l=lsw+1;ljmax;l+=2)
{
resid=a[j][l]*u[j+1][l]+b[j][l]*u[j-1][l]+c[j][l]*u[j][l+1]+d[j][l]*u[j][l-1]+e[j][l]*u[j][l]-f[j][l];
anorm+=fabs(resid);
u[j][l]-=omega*resid/e[j][l];
}
lsw=3-lsw;
}
jsw=3-jsw;
omega=(n==1ipass==1?1.0/(1.0-0.5*rjac*rjac):1.0/(1.0-0.25*rjac*rjac*omega));
}
if(anormEPS*anormf)
return;
}
nrerror("MAXITS exceeded");
}
void nrerror(char error_text[])
{
fprintf(stderr,"Numerical Recipes run-time error...\n");
fprintf(stderr,"%s\n",error_text);
fprintf(stderr,"...now exiting to system...\n");
exit(1);
}
double **dmatrix(long nrl, long nrh, long ncl, long nch)
{
long i, nrow=nrh-nrl+1,ncol=nch-ncl+1;
double **m;
m=(double **) malloc((size_t)((nrow+NR_END)*sizeof(double*)));
if (!m) nrerror("allocation failure 1 in matrix()");
m += NR_END;
m -= nrl;
m[nrl]=(double *) malloc((size_t)((nrow*ncol+NR_END)*sizeof(double)));
if (!m[nrl]) nrerror("allocation failure 2 in matrix()");
m[nrl] += NR_END;
m[nrl] -= ncl;
for(i=nrl+1;i=nrh;i++) m[i]=m[i-1]+ncol;
printf("dmatrix succeeded!\n");
return m;
}
main()
{
double **a,**b,**c,**d,**e,**f,**u,rjac;
int j,l,jmax;
double pi=3.1415926;
long m=10,n=10;
a=dmatrix(1,m,1,n);
b=dmatrix(1,m,1,n);
c=dmatrix(1,m,1,n);
d=dmatrix(1,m,1,n);
e=dmatrix(1,m,1,n);
f=dmatrix(1,m,1,n);
u=dmatrix(1,m,1,n);
printf("%x, %x, %x, %x, %x, %x, %x\n", (int)a, (int)b, (int)c, (int)d, (int)e, (int)f, (int)u);
scanf("%d",jmax);
rjac=cos(pi/jmax);
srand((unsigned)time(NULL));
printf("init begin...\n");
for(j=0;jm;j++)
{
for(l=0;ln;l++)
{
a[j][l]=rand();
b[j][l]=rand();
c[j][l]=rand();
d[j][l]=rand();
e[j][l]=rand();
f[j][l]=rand();
u[j][l]=0.0;
}
}
printf("init finished!\n");
sor(a,b,c,d,e,f,u,jmax,rjac);
for(j=0;j10;j++)
{
for(l=0;l10;l++)
{
printf("%f\t",u[j][l]);
}
printf("\n");
}
getch();
}
迭代法 matlab实现代码如下
function [x,n] = jacobi(A,b,x0,eps,varargin)
if nargin ==3
eps = 1.0e-6;
M = 200;
elseif nargin3
disp('输入参数数目不足3个');
return
elseif nargin ==5
M = varargin{1};
end
D = diag(diag(A)); %%求A的对角矩阵
L = -tril(A,-1); %%求A的下三角矩阵
U = -triu(A,1); %%求A的上三角矩阵
B = D\(L+U);
f = D\b;
x = B*x0+f;
n = 1;%迭代次数
while norm(x-x0)=eps
x0 = x;
x = B*x0+f
n = n+1;
if(n=M)
disp('Warning:迭代次数太多,可能不收敛!')
return;
end
end
运行效果如下：
扩展资料：
迭代法的收敛性判别
收敛性判别条件
SOR迭代法收敛的充分必要条件是ρ(λω)1，ρ(λω)与松弛因子ω有关。ρ(λω)与ω的关系以及SOR方法收敛的条件有如下定理。
定理1：(Kahan)对任意的A
，设其对角元皆非零，则对所有实数ω，有：ρ(λω)≥ ω-1。
推论：如果解Ax=b的SOR方法收敛，则有ω-11，即0ω2。
定理2：(Ostrowski-Reich)设A
，A对称正定，且0ω2，则解Ax=b的SOR方法收敛。
参考资料来源：百度百科-逐次超松驰迭代法
python的话直接使用crypt就行计算的时候要加上salt值。
cat /etc/shadow | greb root
得到密文。 密文由3部分组成，以”$”分隔，第一部分为ID，第二部分为盐值，第三部分为加密密文
密文示例：
$6$D0xsORq3b7GGsbYv$7L8myKFcJf1FqnvtzD.Zcbx9PeHnwGH1YJm/xDeEKmABPIxzUC7mQmEiIEHoHypXeJ969Tfp4c47WQgjkT0YH0
$id$salt$encrypted
采用方法 盐值 加密密文
$6$D0xsORq3b7GGsbYv #这一段就是盐值，在加密的时候连同前面的ID一起加密
import crypt
print crypt.crpyt("password","$6$D0xsORq3b7GGsbYv")
就可以得到shadow中的密文了。
其实shadow中的密文并不是直接由sha512得出的，还有一些其他操作,例如使用base64加密了salt和hash，可以看下这里man page
class Person(object):
name = "张三"
def __init__(self):
self.__age = 18
def set_age(self, age):
self.__age = age
#使用装饰器@classmethod修饰getName()方法
@classmethod
def getName(cls):
return cls.name
#使用类名访问name属性
print(Person().name)
#生成对象p1
p1 = Person()
#使用p1访问name属性
print(p1.name)
#给p1动态绑定属性country
p1.country = 'China'
#获取name
print(p1.getName())
