@@ -156,7 +156,7 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
156156 int elmID = elm2ElmConn (iElm,i)-1 ;
157157
158158 // New delta
159- Vec3d center (elmCenter (elmID, 0 ), elmCenter (elmID, 1 ), elmCenter (elmID, 2 ));
159+ Vec3d center (elmCenter (elmID, 0 ), elmCenter (elmID, 1 ), elmCenter (elmID, 2 ));
160160 delta = MPnew - center;
161161
162162 double neighborDistSq = delta[0 ]*delta[0 ] + delta[1 ]*delta[1 ] + delta[2 ]*delta[2 ];
@@ -165,6 +165,7 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
165165 minDistSq = neighborDistSq;
166166 }
167167 }
168+
168169 if (closestElm<0 ){
169170 MPs2Elm (mp) = iElm;
170171 if (elm2Process.size () > 0 )
@@ -206,6 +207,7 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
206207 Kokkos::deep_copy (h_resultRight, resultRight);
207208 Kokkos::deep_copy (h_mpTgtPos, mpTgtPosArray);
208209 // printVTP file
210+
209211 char * fileOutput = (char *)malloc (sizeof (char ) * 256 );
210212 sprintf (fileOutput, " polyMPOCVTTrackingElmCenter_MPtracks_%d.vtp" , printVTPIndex);
211213 FILE * pFile = fopen (fileOutput," w" );
@@ -229,7 +231,6 @@ void MPMesh::CVTTrackingElmCenterBased(const int printVTPIndex){
229231 fprintf (pFile," </DataArray>\n </Lines>\n </Piece>\n </PolyData>\n </VTKFile>\n " );
230232 fclose (pFile);
231233 }
232-
233234 pumipic::RecordTime (" PolyMPO_CVTTrackingElmCenterBased" , timer.seconds ());
234235}
235236
0 commit comments