- 5 名前：デフォルトの名無しさん mailto:sage [2009/07/19(日) 23:56:35 ]
- >>2
こんなか？
double [][]a = fileInMatrix(inputfolder,matrix_in); //行列用の読み込み
double []b = fileInVector(inputfolder,vector_in); //ベクトル用の読み込み
double []x = fileInVector(inputfolder,vector2_in); //ベクトル用の読み込み
double []solution=new double[a.length];
for(int i=0;i<a.length;i++){
solution[i]=0;
for(int j=0;j<x.length;j++){
solution[i]+=a[i][j]*x[j];
}
}
for(int i=0;i<solution.length;i++){
solution[i]-=b[i];
}
fileOutVector(outputfolder,vector_out,shutsuryoku,solution); //ベクトル用の書き込み
double euqlid=0;
for(double d :solution){
euqlid+=d*d;
}
System.out.println("e-norm:"+Math.sqrt(euqlid));
double max = solution[0];
for(int i=1;i<solution.length;i++){
if(max<solution[i]){
max=Math.abs(solution[i]);
}
}
System.out.println("max-norm:"+max);
|
 |
