{
}
}
u=i*n+k;v=i*n+is[k];
p=a[i][k];a[i][k]=a[i][is[k]];a[i][is[k]]=p;
free(is);free(js); return (1);
} int invGJ(double **a,int n)
int *is,*js,i,j,k,l,u,v; double d,p;
is=(int *)malloc(n*sizeof(int)); js=(int *)malloc(n*sizeof(int)); for(k=0;k<=n-1;k++) { { }
if(is[k]!=k)
for(j=0;j<=n-1;j++) {
u=k*n+j;v=is[k]*n+j;
p=a[k][j];a[k][j]=a[is[k]][j];a[is[k]][j]=p;
d=0.0;
for(i=k;i<=n-1;i++) for(j=k;j<=n-1;j++) { }
l=i*n+j;p=fabs(a[i][j]); if(p>d)
{d=p;is[k]=i;js[k]=j;}
if(d+1.0==1.0)
free(is);free(js);printf("error not inv\n"); return (0);
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说教育文库测量平差编程(11)在线全文阅读。
相关推荐: