int ginv(double* a,int m,int n,double aa[],double eps,double u[],double v[],int ka)
{
int i,j,k,l,t,p,q,f;
i=muav(a,m,n,u,v,eps,ka);
if(i<0) return(-1);
j=n;
if(m<n) j=m;
j=j-1;
k=0;
while((k<=j)&&(a[k*n+k]!=0.0)) k=k+1;
k=k-1;
for(i=0;i<=n-1;i++)
for(j=0;j<=m-1;j++)
{
t=i*m+j;
aa[t]=0.0;
for(l=0;l<=k;l++)
{
f=l*n+i;p=j*m+l;q=l*n+l;
aa[t]=aa[t]+v[f]*u[p]/a[q];
}
}
return(1);
}
main()
{
int i, j;
double a[4][3] = { {1.0, 1.0, -1.0}, {2.0, 1.0, 0.0},
{1.0, -1.0, 0.0}, { -1.0, 2.0, 1.0}
};
double b[3][4] = { {1.0, 1.0, -1.0, -1.0}, {
2.0, 1.0,
0.0, 2.0
}, {1.0, -1.0, 0.0, 1.0}
};
double u[4][4], v[3][3], c[4][3], d[3][4];
double eps;
eps = 0.000001;
i = muav(a, 4, 3, u, v, eps, 5);
printf("\n");
printf("EXAMPLE(1)\n");
printf("\n");
printf("i=%d\n", i);
printf("\n");
printf("MAT U IS:\n");
for (i = 0; i <= 3; i++) {
for (j = 0; j <= 3; j++) {
printf("%13.6e ", u[i][j]);
}
printf("\n");
}
printf("\n");
printf("MAT V IS:\n");
for (i = 0; i <= 2; i++) {
for (j = 0; j <= 2; j++) {
printf("%13.6e ", v[i][j]);
}
printf("\n");
}
printf("\n");
printf("MAT A IS:\n");
for (i = 0; i <= 3; i++) {
for (j = 0; j <= 2; j++) {
printf("%13.6e ", a[i][j]);
}
printf("\n");
}
printf("\n\n");
printf("MAT UAV IS:\n");
trmul(u, a, 4, 4, 3, c);
trmul(c, v, 4, 3, 3, a);
for (i = 0; i <= 3; i++) {
for (j = 0; j <= 2; j++) {
printf("%13.6e ", a[i][j]);
}
printf("\n");
}
printf("\n\n");
printf("EXAMPLE(2)\n");
printf("\n");
i = muav(b, 3, 4, v, u, eps, 5);
printf("i=%d\n", i);
printf("\n");
printf("MAT U IS:\n");
for (i = 0; i <= 2; i++) {
for (j = 0; j <= 2; j++) {
printf("%13.6e ", v[i][j]);
}
printf("\n");
}
printf("\n");
printf("MAT V IS:\n");
for (i = 0; i <= 3; i++) {
for (j = 0; j <= 3; j++) {
printf("%13.6e ", u[i][j]);
}
printf("\n");
}
printf("\n");
printf("MAT B IS:\n");
for (i = 0; i <= 2; i++) {
for (j = 0; j <= 3; j++) {
printf("%13.6e ", b[i][j]);
}
printf("\n");
}
printf("\n\n");
printf("MAT UBV IS:\n");
trmul(v, b, 3, 3, 4, d);
trmul(d, u, 3, 4, 4, b);
for (i = 0; i <= 2; i++) {
for (j = 0; j <= 3; j++) {
printf("%13.6e ", b[i][j]);
}
printf("\n");
}
printf("\n");
}
