Commit 72d3c5e9 authored by Andreas Marek's avatar Andreas Marek

Call proper double or single precision routines

parent 2cb6e81a
......@@ -242,7 +242,7 @@ program test_complex_double_precision
allocate(xr(na_rows,na_cols))
call prepare_matrix(na, myid, sc_desc, iseed, xr, a, z, as)
call prepare_matrix_double(na, myid, sc_desc, iseed, xr, a, z, as)
deallocate(xr)
#ifdef HAVE_DETAILED_TIMINGS
......
......@@ -256,7 +256,7 @@ program test_complex2_double_precision
allocate(ev(na))
allocate(xr(na_rows,na_cols))
call prepare_matrix(na, myid, sc_desc, iseed, xr, a, z, as)
call prepare_matrix_double(na, myid, sc_desc, iseed, xr, a, z, as)
deallocate(xr)
......
......@@ -233,7 +233,7 @@ program test_real_double_precision
allocate(ev(na))
call prepare_matrix(na, myid, sc_desc, iseed, a, z, as)
call prepare_matrix_double(na, myid, sc_desc, iseed, a, z, as)
#ifdef HAVE_DETAILED_TIMINGS
call timer%stop("set up matrix")
......
......@@ -246,7 +246,7 @@ program test_real2_double_precision
allocate(ev(na))
call prepare_matrix(na, myid, sc_desc, iseed, a, z, as)
call prepare_matrix_double(na, myid, sc_desc, iseed, a, z, as)
#ifdef HAVE_DETAILED_TIMINGS
call timer%stop("set up matrix")
......
......@@ -207,10 +207,10 @@ program test_transpose_multiply
allocate(ev(na))
allocate(xr(na_rows,na_cols))
call prepare_matrix(na, myid, sc_desc, iseed, xr, a, z, as)
call prepare_matrix_double(na, myid, sc_desc, iseed, xr, a, z, as)
deallocate(xr)
b(:,:) = 2.0 * a(:,:)
c(:,:) = 0.0
b(:,:) = 2.0_ck8 * a(:,:)
c(:,:) = 0.0_ck8
#ifdef HAVE_DETAILED_TIMINGS
call timer%stop("set up matrix")
......
......@@ -209,8 +209,8 @@ program test_transpose_multiply
allocate(xr(na_rows,na_cols))
call prepare_matrix_single(na, myid, sc_desc, iseed, xr, a, z, as)
deallocate(xr)
b(:,:) = 2.0 * a(:,:)
c(:,:) = 0.0
b(:,:) = 2.0_ck4 * a(:,:)
c(:,:) = 0.0_ck4
#ifdef HAVE_DETAILED_TIMINGS
call timer%stop("set up matrix")
......
......@@ -206,11 +206,11 @@ program test_transpose_multiply
allocate(ev(na))
call prepare_matrix(na, myid, sc_desc, iseed, a, z, as)
call prepare_matrix_double(na, myid, sc_desc, iseed, a, z, as)
b(:,:) = 2.0 * a(:,:)
c(:,:) = 0.0
b(:,:) = 2.0_rk8 * a(:,:)
c(:,:) = 0.0_rk8
#ifdef HAVE_DETAILED_TIMINGS
call timer%stop("set up matrix")
......
......@@ -209,8 +209,8 @@ program test_transpose_multiply
call prepare_matrix_single(na, myid, sc_desc, iseed, a, z, as)
b(:,:) = 2.0 * a(:,:)
c(:,:) = 0.0
b(:,:) = 2.0_rk4 * a(:,:)
c(:,:) = 0.0_rk4
#ifdef HAVE_DETAILED_TIMINGS
call timer%stop("set up matrix")
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment