|
| |
 精品推荐 |
 |
|
| |
|
|
|
|
数值计算程序大放送-线性代数方程组
|
日期：2008年3月31日 作者： 查看:[大字体
中字体 小字体]
|
{
d=p;js[k]=j;is=i;
}
}
}
if (d+1.0==1.0)
{
free(js);
printf("err**fail\n");
return(0);
}
if (is!=k)
{
for (j=k;j<=n-1;j++)
{
u=k*n+j; v=is*n+j;
p=ar[u]; ar[u]=ar[v]; ar[v]=p;
p=ai[u]; ai[u]=ai[v]; ai[v]=p;
}
for (j=0;j<=m-1;j++)
{
u=k*m+j; v=is*m+j;
p=br[u]; br[u]=br[v]; br[v]=p;
p=bi[u]; bi[u]=bi[v]; bi[v]=p;
}
}
if (js[k]!=k)
{
for (i=0;i<=n-1;i++)
{
u=i*n+k; v=i*n+js[k];
p=ar[u]; ar[u]=ar[v]; ar[v]=p;
p=ai[u]; ai[u]=ai[v]; ai[v]=p;
}
}
v=k*n+k;
for (j=k+1;j<=n-1;j++)
{
u=k*n+j;
p=ar[u]*ar[v]; q=-ai[u]*ai[v];
s=(ar[v]-ai[v])*(ar[u]+ai[u]);
ar[u]=(p-q)/d; ai[u]=(s-p-q)/d;
}
for (j=0;j<=m-1;j++)
{
u=k*m+j;
p=br[u]*ar[v]; q=-bi[u]*ai[v];
上一篇:IIS 6.0配置HTTP压缩
下一篇:微软怎么设计Exchange 2003的体系结构
|
| 相关文章： |
|
|
|
| 相关软件： |
|
| |
|
