錯誤不在sor，在這里：
創(chuàng)新互聯(lián)-專業(yè)網站定制、快速模板網站建設、高性價比潞城網站開發(fā)、企業(yè)建站全套包干低至880元,成熟完善的模板庫,直接使用。一站式潞城網站制作公司更省心,省錢,快速模板網站建設找我們，業(yè)務覆蓋潞城地區(qū)。費用合理售后完善，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。環(huán)境是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動態(tài)綁定屬性country
p1.country = 'China'
#獲取name
print(p1.getName())
