请c/c++高手看一下关于sor程序的问题

错误不在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();

}

SOR迭代法 求MATLAB程序

迭代法 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 hashlib.sha512 为什么加密结果和shadow不一样

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

python类的题目求帮忙看下

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())

本文名称：sor函数python sorted函数python
分享URL：http://scgulin.cn/article/dopeccs.html