OSDN Git Service

[バックアップ]NyARToolkit for Java
[nyartoolkit-and/nyartoolkit-and.git] / src / jp / nyatla / nyartoolkit / core / NyARVec.java
index 8bc84d7..7c43f02 100644 (file)
@@ -156,121 +156,6 @@ public class NyARVec
                return -s;\r
        }\r
 \r
-       // /**\r
-       // * arVecTridiagonalize関数の代替品\r
-       // * a,d,e間で演算をしてる。何をどうしているかはさっぱりさっぱり\r
-       // * @param a\r
-       // * @param d\r
-       // * @param e\r
-       // * @param i_e_start\r
-       // * 演算開始列(よくわからないけどarVecTridiagonalizeの呼び出し元でなんかしてる)\r
-       // * @return\r
-       // * @throws NyARException\r
-       // */\r
-       // public static void vecTridiagonalize(NyARMat a, NyARVec d, NyARVec e,int\r
-       // i_e_start) throws NyARException\r
-       // {\r
-       // NyARVec vec,vec2;\r
-       // double[][] a_array=a.getArray();\r
-       // double s, t, p, q;\r
-       // int dim;\r
-       //\r
-       // if(a.getClm()!=a.getRow()){\r
-       // throw new NyARException();\r
-       // }\r
-       // if(a.getClm() != d.clm){\r
-       // throw new NyARException();\r
-       // }\r
-       // if(a.getClm() != e.clm){\r
-       // throw new NyARException();\r
-       // }\r
-       // dim = a.getClm();\r
-       //\r
-       // for(int k = 0; k < dim-2; k++ ){\r
-       // vec=a.getRowVec(k);\r
-       // // double[] vec_array=vec.getArray();\r
-       // NyARException.trap("未チェックパス");\r
-       // d.v[k]=vec.v[k];//d.set(k,v.get(k)); //d->v[k] = v[k];\r
-       //\r
-       // //wv1.clm = dim-k-1;\r
-       // //wv1.v = &(v[k+1]);\r
-       // NyARException.trap("未チェックパス");\r
-       // e.v[k+i_e_start]=vec.vecHousehold(k+1);//e->v[k] = arVecHousehold(&wv1);\r
-       // if(e.v[k+i_e_start]== 0.0 ){\r
-       // continue;\r
-       // }\r
-       //\r
-       // for(int i = k+1; i < dim; i++ ){\r
-       // s = 0.0;\r
-       // for(int j = k+1; j < i; j++ ) {\r
-       // NyARException.trap("未チェックのパス");\r
-       // s += a_array[j][i] * vec.v[j];//s += a.get(j*dim+i) * v.get(j);//s +=\r
-       // a->m[j*dim+i] * v[j];\r
-       // }\r
-       // for(int j = i; j < dim; j++ ) {\r
-       // NyARException.trap("未チェックのパス");\r
-       // s += a_array[i][j] * vec.v[j];//s += a.get(i*dim+j) * v.get(j);//s +=\r
-       // a->m[i*dim+j] * v[j];\r
-       // }\r
-       // NyARException.trap("未チェックのパス");\r
-       // d.v[i]=s;//d->v[i] = s;\r
-       // }\r
-       //\r
-       //\r
-       // //wv1.clm = wv2.clm = dim-k-1;\r
-       // //wv1.v = &(v[k+1]);\r
-       // //wv2.v = &(d->v[k+1]);\r
-       // vec=a.getRowVec(k);\r
-       // // vec_array=vec.getArray();\r
-       // NyARException.trap("未チェックパス");\r
-       // t = vec.vecInnerproduct(d,k+1)/ 2;\r
-       // for(int i = dim-1; i > k; i-- ) {\r
-       // NyARException.trap("未チェックパス");\r
-       // p = vec.v[i];//p = v.get(i);//p = v[i];\r
-       // d.v[i]-=t*p;q=d.v[i];//q = d->v[i] -= t*p;\r
-       // for(int j = i; j < dim; j++ ){\r
-       // NyARException.trap("未チェックパス");\r
-       // a_array[i][j]-=p*(d.v[j] + q*vec.v[j]);//a->m[i*dim+j] -= p*(d->v[j]) +\r
-       // q*v[j];\r
-       // }\r
-       // }\r
-       // }\r
-       //\r
-       // if( dim >= 2) {\r
-       // d.v[dim-2]=a_array[dim-2][dim-2];//d->v[dim-2] =\r
-       // a->m[(dim-2)*dim+(dim-2)];\r
-       // e.v[dim-2+i_e_start]=a_array[dim-2][dim-1];//e->v[dim-2] =\r
-       // a->m[(dim-2)*dim+(dim-1)];\r
-       // }\r
-       //\r
-       // if( dim >= 1 ){\r
-       // d.v[dim-1]=a_array[dim-1][dim-1];//d->v[dim-1] =\r
-       // a->m[(dim-1)*dim+(dim-1)];\r
-       // }\r
-       //\r
-       // for(int k = dim-1; k >= 0; k--) {\r
-       // vec=a.getRowVec(k);//v = a.getPointer(k*dim);//v = &(a->m[k*dim]);\r
-       // if( k < dim-2 ) {\r
-       // for(int i = k+1; i < dim; i++ ){\r
-       // //wv1.clm = wv2.clm = dim-k-1;\r
-       // //wv1.v = &(v[k+1]);\r
-       // //wv2.v = &(a->m[i*dim+k+1]);\r
-       // vec2=a.getRowVec(i);\r
-       //\r
-       // t = vec.vecInnerproduct(vec2,k+1);\r
-       // for(int j = k+1; j < dim; j++ ){\r
-       // NyARException.trap("未チェックパス");\r
-       // a_array[i][j]-=t*vec.v[j];//a.subValue(i*dim+j,t*v.get(j));//a->m[i*dim+j]\r
-       // -= t * v[j];\r
-       // }\r
-       // }\r
-       // }\r
-       // for(int i = 0; i < dim; i++ ){\r
-       // vec.v[i]=0.0;//v.set(i,0.0);//v[i] = 0.0;\r
-       // }\r
-       // vec.v[k]=1;//v.set(k,1);//v[k] = 1;\r
-       // }\r
-       // }\r
        /**\r
         * 現在ラップしている配列を取り外して、新しい配列をラップします。\r
         * \r