diff --git a/examples/ir-ilu-preconditioned-solver/CMakeLists.txt b/examples/ir-ilu-preconditioned-solver/CMakeLists.txt index 78da3312512..dd77e163e59 100644 --- a/examples/ir-ilu-preconditioned-solver/CMakeLists.txt +++ b/examples/ir-ilu-preconditioned-solver/CMakeLists.txt @@ -2,5 +2,3 @@ add_executable(ir-ilu-preconditioned-solver ir-ilu-preconditioned-solver.cpp) target_link_libraries(ir-ilu-preconditioned-solver ginkgo) target_include_directories(ir-ilu-preconditioned-solver PRIVATE ${PROJECT_SOURCE_DIR}) configure_file(data/A.mtx data/A.mtx COPYONLY) -configure_file(data/b.mtx data/b.mtx COPYONLY) -configure_file(data/x0.mtx data/x0.mtx COPYONLY) diff --git a/examples/ir-ilu-preconditioned-solver/data/b.mtx b/examples/ir-ilu-preconditioned-solver/data/b.mtx deleted file mode 100644 index 05d92ecc6f7..00000000000 --- a/examples/ir-ilu-preconditioned-solver/data/b.mtx +++ /dev/null @@ -1,21 +0,0 @@ -%%MatrixMarket matrix array real general -19 1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 diff --git a/examples/ir-ilu-preconditioned-solver/data/x0.mtx b/examples/ir-ilu-preconditioned-solver/data/x0.mtx deleted file mode 100644 index 91d470cdbcd..00000000000 --- a/examples/ir-ilu-preconditioned-solver/data/x0.mtx +++ /dev/null @@ -1,21 +0,0 @@ -%%MatrixMarket matrix array real general -19 1 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 diff --git a/examples/ir-ilu-preconditioned-solver/doc/results.dox b/examples/ir-ilu-preconditioned-solver/doc/results.dox index 5de3e3b01f5..eaaaa5758cd 100644 --- a/examples/ir-ilu-preconditioned-solver/doc/results.dox +++ b/examples/ir-ilu-preconditioned-solver/doc/results.dox @@ -26,11 +26,11 @@ Solution (x): 0.0121141 0.0123025 GMRES iteration count: 7 -GMRES execution time [ms]: 5.49966 +GMRES execution time [ms]: 2.64993 Residual norm sqrt(r^T r): %%MatrixMarket matrix array real general 1 1 -7.83607e-13 +2.23805e-10 @endcode diff --git a/examples/ir-ilu-preconditioned-solver/ir-ilu-preconditioned-solver.cpp b/examples/ir-ilu-preconditioned-solver/ir-ilu-preconditioned-solver.cpp index e688e8573d9..14cadf8c912 100644 --- a/examples/ir-ilu-preconditioned-solver/ir-ilu-preconditioned-solver.cpp +++ b/examples/ir-ilu-preconditioned-solver/ir-ilu-preconditioned-solver.cpp @@ -76,8 +76,16 @@ int main(int argc, char *argv[]) // Read data auto A = gko::share(gko::read(std::ifstream("data/A.mtx"), exec)); - auto b = gko::read(std::ifstream("data/b.mtx"), exec); - auto x = gko::read(std::ifstream("data/x0.mtx"), exec); + // Create RHS and initial guess as 1 + gko::size_type num_rows = A->get_size()[0]; + auto host_x = vec::create(exec->get_master(), gko::dim<2>(num_rows, 1)); + for (gko::size_type i = 0; i < num_rows; i++) { + host_x->at(i, 0) = 1.; + } + auto x = vec::create(exec); + auto b = vec::create(exec); + x->copy_from(host_x.get()); + b->copy_from(host_x.get()); auto clone_x = vec::create(exec); clone_x->copy_from(lend(x));