update

Author: RaulPPelaez

README.md

@@ -50,7 +50,7 @@ Provide the solver with an instance of the functor and the target vector:
     thrust::fill(v.begin(), v.end(), 1);
     //A vector to store the result of sqrt(M)*v
     thrust::device_vector<real> result(size);
-    //A functor that multiplies by the identity matrix times two
+    //A functor that multiplies by a diagonal matrix
     MatrixDot dot(size);
     //Call the solver
     real* d_result = thrust::raw_pointer_cast(result.data());

example.cpp

@@ -45,7 +45,7 @@ int main(){
     std::fill(v.begin(), v.end(), 1);
     //A vector to store the result of sqrt(M)*v
     std::vector<real> result(size);
-    //A functor that multiplies by the identity matrix times two
+    //A functor that multiplies by a diagonal matrix
     MatrixDot dot(size);
     //Call the solver
     int numberIterations = lanczos.solve(dot, result.data(), v.data(), size);

example.cu

@@ -45,7 +45,7 @@ int main(){
     lanczos::detail::device_fill(v.begin(), v.end(), 1);
     //A vector to store the result of sqrt(M)*v
     lanczos::device_container<real> result(size);
-    //A functor that multiplies by the identity matrix times two
+    //A functor that multiplies by ta diagonal matrix
     MatrixDot dot(size);
     //Call the solver
     real* d_result = lanczos::detail::getRawPointer(result);