?? main.cpp
字號:
#include "svd.h"
#include "stdio.h"
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 a[12]={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[12]={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[16],v[9],c[12],d[12];
for(i=0;i<16;i++)
u[i]=0;
double eps=0.000001;
i=dluav(a,4,3,u,v,eps,5);
printf("\n");
printf("i=%d\n",i);
printf("\nMAT U Is:\n");
for(i=0;i<=3;i++)
{
for(j=0;j<=3;j++)
printf("%e ",u[i*4+j]);
printf("\n");
}
printf("\n");
printf("MAT V IS:\n");
for(i=0;i<=2;i++)
{
for(j=0;j<=2;j++)
printf("%e ",v[i*3+j]);
printf("\n");
}
printf("\n");
printf("MAT A Is:\n");
for(i=0;i<=3;i++)
{
for(j=0;j<=2;j++)
printf("%e ",a[i*3+j]);
printf("\n");
}
damul(u,a,4,4,3,c);
damul(c,v,4,3,3,a);
printf("\nMAT UAV IS:\n");
for(i=0;i<=3;i++)
{
for(j=0;j<=2;j++)
printf("%e ",a[i*3+j]);
printf("\n");
}
printf("\n\n");
printf("第二個例子\n");
i=dluav(b,3,4,v,u,eps,5);
printf("i=%d\n",i);
printf("\n");
printf("MATUS IS:\n");
for(i=0;i<=2;i++)
{
for(j=0;j<=2;j++)
printf("%e ",v[i*3+j]);
printf("\n");
}
printf("\n");
printf("MAT V IS:\n");
for(i=0;i<=3;i++)
{
for(j=0;j<=3;j++)
printf("%e ",u[i*3+j]);
printf("\n");
}
printf("\nMAT B IS:\n");
for(i=0;i<=2;i++)
{
for(j=0;j<=3;j++)
printf("%e ",b[i*4+j]);
printf("\n");
}
printf("\n\n");
printf("MAT UBV IS:\n");
damul(v,b,3,3,4,d);
damul(d,u,3,4,4,b);
for(i=0;i<=2;i++)
{
for(j=0;j<=3;j++)
printf("%e ",b[i*4+j]);
printf("\n");
}
printf("\n");
}
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -