diff --git a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlDirectEikonalMethodsBase_impl.h b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlDirectEikonalMethodsBase_impl.h
index c742a11b501c818beccf88b94c1671b133c4b8cb..cd4eb9e24b81373b79335227e67808b6a19fa228 100644
--- a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlDirectEikonalMethodsBase_impl.h
+++ b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlDirectEikonalMethodsBase_impl.h
@@ -155,7 +155,7 @@ initInterface( const MeshFunctionPointer& _input,
     {
         MeshFunctionType input = _input.getData();
         
-        double A[320][320];
+        /*double A[320][320];
         std::ifstream fileInit("/home/maty/Downloads/initData.txt");
 
         for (int i = 0; i < 320; i++)
@@ -164,7 +164,7 @@ initInterface( const MeshFunctionPointer& _input,
         fileInit.close();
         for (int i = 0; i < 320; i++)
             for (int j = 0; j < 320; j++)
-                input[i*320 + j] = A[i][j];
+                input[i*320 + j] = A[i][j];*/
         
         
          MeshFunctionType& output = _output.modifyData();
diff --git a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlFastSweepingMethod2D_impl.h b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlFastSweepingMethod2D_impl.h
index e5a9c39f30ace9f19f3a5d527550fb8d1db9f6f8..83333b190d20c99b4867ab128231f5eaed6d4a7d 100644
--- a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlFastSweepingMethod2D_impl.h
+++ b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/tnlFastSweepingMethod2D_impl.h
@@ -64,7 +64,7 @@ solve( const MeshPointer& mesh,
        const AnisotropyPointer& anisotropy,
        MeshFunctionPointer& u )
 {
-   MeshFunctionType v;
+   /*MeshFunctionType v;
    v.setMesh(mesh);
    double A[320][320];
     for (int i = 0; i < 320; i++)
@@ -80,7 +80,7 @@ solve( const MeshPointer& mesh,
     for (int i = 0; i < 320; i++)
         for (int j = 0; j < 320; j++)
             v[i*320 + j] = A[i][j];
-   v.save("mapa.tnl");
+   v.save("mapa.tnl");*/
    
        
    MeshFunctionPointer auxPtr;
@@ -111,7 +111,7 @@ solve( const MeshPointer& mesh,
                {
                   cell.refresh();
                   if( ! interfaceMap( cell ) )
-                     this->updateCell( aux, cell, v( cell ) );
+                     this->updateCell( aux, cell );
                }
          }
 
@@ -128,7 +128,7 @@ solve( const MeshPointer& mesh,
                   //std::cerr << "2 -> ";
                   cell.refresh();
                   if( ! interfaceMap( cell ) )            
-                     this->updateCell( aux, cell, v( cell ) );
+                     this->updateCell( aux, cell );
                }
          }
 
@@ -145,7 +145,7 @@ solve( const MeshPointer& mesh,
                   //std::cerr << "3 -> ";
                   cell.refresh();
                   if( ! interfaceMap( cell ) )            
-                     this->updateCell( aux, cell, v( cell ) );
+                     this->updateCell( aux, cell );
                }
             }
 
@@ -162,7 +162,7 @@ solve( const MeshPointer& mesh,
                   //std::cerr << "4 -> ";
                   cell.refresh();
                   if( ! interfaceMap( cell ) )            
-                     this->updateCell( aux, cell, v( cell ) );
+                     this->updateCell( aux, cell );
                }
             }