@@ -54,7 +54,6 @@ RegisterWithCorrespondencesAlgo::RegisterWithCorrespondencesAlgo()
5454
5555AlgorithmOutput RegisterWithCorrespondencesAlgo::run_generic (const AlgorithmInput& input) const
5656{
57- std::cout << " ALGO" << std::endl;
5857 auto input_field = input.get <Field>(Variables::InputField);
5958 auto corres1 = input.get <Field>(Correspondences1);
6059 auto corres2 = input.get <Field>(Correspondences2);
@@ -65,15 +64,12 @@ AlgorithmOutput RegisterWithCorrespondencesAlgo::run_generic(const AlgorithmInpu
6564 switch (op)
6665 {
6766 case 0 :
68- std::cout << " runM" << std::endl;
6967 runM (input_field, corres1, corres2, return_field);
7068 break ;
7169 case 1 :
72- std::cout << " runA" << std::endl;
7370 runA (input_field, corres1, corres2, return_field);
7471 break ;
7572 case 2 :
76- std::cout << " runN" << std::endl;
7773 runN (input_field, corres1, corres2, return_field);
7874 break ;
7975 }
@@ -227,10 +223,10 @@ bool RegisterWithCorrespondencesAlgo::runM(FieldHandle input, FieldHandle Cors1,
227223 Bm (3 , L1) = 1 ;
228224
229225 // vertical x,y,z
230- Bm (L1 + 4 ) = num_cors1, P.x ();
231- Bm (L1 + 4 ) = num_cors1 + 1 , P.y ();
232- Bm (L1 + 4 ) = num_cors1 + 2 , P.z ();
233- Bm (L1 + 4 ) = num_cors1 + 3 , 1 ;
226+ Bm (L1 + 4 , num_cors1 ) = P.x ();
227+ Bm (L1 + 4 , num_cors1 + 1 ) = P.y ();
228+ Bm (L1 + 4 , num_cors1 + 2 ) = P.z ();
229+ Bm (L1 + 4 , num_cors1 + 3 ) = 1 ;
234230 }
235231
236232 for (int L1 = num_cors1; L1 < num_cors1 + 4 ; ++L1)
@@ -244,7 +240,6 @@ bool RegisterWithCorrespondencesAlgo::runM(FieldHandle input, FieldHandle Cors1,
244240 // put in sigmas
245241 DenseMatrixHandle SMat;
246242 radial_basis_func (icors2, icors2, SMat);
247- double temp = 0 ;
248243
249244 for (int i = 0 ; i < num_cors1; ++i)
250245 {
@@ -754,7 +749,6 @@ bool RegisterWithCorrespondencesAlgo::make_new_pointsA(VMesh* points, VMesh* Cor
754749 Point P, Pp;
755750
756751 points->size (num_pts);
757- int msz = coefs.size ();
758752
759753 for (int i = 0 ; i < num_pts; ++i)
760754 {
0 commit comments