Skip to content
Merged
Show file tree
Hide file tree
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
review update
* Rename multiprecision -> mixed precision
* Improve documentation and output
* Improve code style

Co-authored-by: Hartwig Anzt <hanzt@icl.utk.edu>
Co-authored-by: Yuhsiang M. Tsai <yhmtsai@gmail.com>
Co-authored-by: Pratik Nayak <pratikvn@protonmail.com>
Co-authored-by: Terry Cojean <terry.cojean@kit.edu>
  • Loading branch information
5 people committed May 4, 2020
commit c373ed06e203196fc54b06dcaee4c0400b1bd5bc
2 changes: 1 addition & 1 deletion examples/CMakeLists.txt
Original file line number Diff line number Diff line change
Expand Up @@ -12,7 +12,7 @@ add_subdirectory(ilu-preconditioned-solver)
add_subdirectory(inverse-iteration)
add_subdirectory(iterative-refinement)
add_subdirectory(minimal-cuda-solver)
add_subdirectory(multiprecision-ir)
add_subdirectory(mixed-precision-ir)
add_subdirectory(nine-pt-stencil-solver)
add_subdirectory(papi-logging)
add_subdirectory(performance-debugging)
Expand Down
4 changes: 4 additions & 0 deletions examples/mixed-precision-ir/CMakeLists.txt
Original file line number Diff line number Diff line change
@@ -0,0 +1,4 @@
add_executable(mixed-precision-ir mixed-precision-ir.cpp)
target_link_libraries(mixed-precision-ir ginkgo)
target_include_directories(mixed-precision-ir PRIVATE ${PROJECT_SOURCE_DIR})
configure_file(data/A.mtx data/A.mtx COPYONLY)
Original file line number Diff line number Diff line change
Expand Up @@ -32,7 +32,7 @@ if [ -z "${CXX}" ]; then
fi

# build
${CXX} -std=c++11 -o ${THIS_DIR}/multiprecision-ir \
${THIS_DIR}/multiprecision-ir.cpp \
${CXX} -std=c++11 -o ${THIS_DIR}/mixed-precision-ir \
${THIS_DIR}/mixed-precision-ir.cpp \
-I${THIS_DIR}/../../include -I${BUILD_DIR}/include \
-L${THIS_DIR} ${LINK_FLAGS}
1 change: 1 addition & 0 deletions examples/mixed-precision-ir/doc/builds-on
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
iterative-refinement
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
<a name="Multiprecision Iterative Refinement"></a>
<h1>This example manually implements a multiprecision iterative refinement solver.</h1>
<a name="Mixed Precision Iterative Refinement (MPIR)"></a>
<h1>This example manually implements a Mixed Precision Iterative Refinement (MPIR) solver.</h1>

<h3> In this example, we first read in a matrix from file, then generate a
right-hand side and an initial guess. An inaccurate CG solver in single precision
Expand Down
1 change: 1 addition & 0 deletions examples/mixed-precision-ir/doc/kind
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
techniques
Original file line number Diff line number Diff line change
Expand Up @@ -3,16 +3,16 @@ This is the expected output:

@code{.cpp}

Initial residual norm sqrt(r^T r):
Initial residual norm sqrt(r^T r):
%%MatrixMarket matrix array real general
1 1
194.679
Final residual norm sqrt(r^T r):
Final residual norm sqrt(r^T r):
%%MatrixMarket matrix array real general
1 1
4.23821e-11
IR iteration count: 24
IR execution time [ms]: 18.0692
1.22728e-10
MPIR iteration count: 25
MPIR execution time [ms]: 18.0933

@endcode

Expand Down
1 change: 1 addition & 0 deletions examples/mixed-precision-ir/doc/short-intro
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
The Mixed Precision Iterative Refinement (MPIR) solver example.
1 change: 1 addition & 0 deletions examples/mixed-precision-ir/doc/tooltip
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
Manually implement a Mixed Precision Iterative Refinement (MPIR) method in Ginkgo. Solve a linear system.
Original file line number Diff line number Diff line change
Expand Up @@ -81,13 +81,12 @@ int main(int argc, char *argv[])
auto A = share(gko::read<mtx>(std::ifstream("data/A.mtx"), exec));
// Create RHS and initial guess as 1
gko::size_type size = A->get_size()[0];
auto host_x = gko::matrix::Dense<ValueType>::create(exec->get_master(),
gko::dim<2>(size, 1));
auto host_x = vec::create(exec->get_master(), gko::dim<2>(size, 1));
for (auto i = 0; i < size; i++) {
host_x->at(i, 0) = 1.;
}
auto x = gko::matrix::Dense<ValueType>::create(exec);
auto b = gko::matrix::Dense<ValueType>::create(exec);
auto x = vec::create(exec);
auto b = vec::create(exec);
x->copy_from(host_x.get());
b->copy_from(host_x.get());

Expand Down Expand Up @@ -172,7 +171,7 @@ int main(int argc, char *argv[])
write(std::cout, lend(res_vec));

// Print solver statistics
std::cout << "IR iteration count: " << iter << std::endl;
std::cout << "IR execution time [ms]: "
std::cout << "MPIR iteration count: " << iter << std::endl;
std::cout << "MPIR execution time [ms]: "
<< static_cast<double>(time.count()) / 1000000.0 << std::endl;
}
4 changes: 0 additions & 4 deletions examples/multiprecision-ir/CMakeLists.txt

This file was deleted.

1 change: 0 additions & 1 deletion examples/multiprecision-ir/doc/builds-on

This file was deleted.

1 change: 0 additions & 1 deletion examples/multiprecision-ir/doc/kind

This file was deleted.

1 change: 0 additions & 1 deletion examples/multiprecision-ir/doc/short-intro

This file was deleted.

1 change: 0 additions & 1 deletion examples/multiprecision-ir/doc/tooltip

This file was deleted.