Advanced Computing Platform for Theoretical Physics

Commit 8a52d6ab authored by mikeaclark's avatar mikeaclark
Browse files

Small bug fix in invert_test

git-svn-id: http://lattice.bu.edu/qcdalg/cuda/quda@445 be54200a-260c-0410-bdd7-ce6af2a381ab
parent b5b5f4df
...@@ -3,7 +3,7 @@ CCFILES := inv_bicgstab_cuda.cpp inv_cg.cpp util_cuda.cpp gauge_cuda.cpp spinor ...@@ -3,7 +3,7 @@ CCFILES := inv_bicgstab_cuda.cpp inv_cg.cpp util_cuda.cpp gauge_cuda.cpp spinor
CUDA_INSTALL_PATH = /usr/local/cuda CUDA_INSTALL_PATH = /usr/local/cuda
INCLUDES = -I. -I$(CUDA_INSTALL_PATH)/include -I$(CUDA_INSTALL_PATH)/../cuda_sdk/common/inc INCLUDES = -I. -I$(CUDA_INSTALL_PATH)/include -I$(CUDA_INSTALL_PATH)/../cuda_sdk/common/inc
LIB = -L$(CUDA_INSTALL_PATH)/lib -lcudart LIB = -L$(CUDA_INSTALL_PATH)/lib64 -lcudart
DFLAGS = #-D__DEVICE_EMULATION__ DFLAGS = #-D__DEVICE_EMULATION__
......
...@@ -30,10 +30,10 @@ int TRANSFER = 0; // include transfer time in the benchmark? ...@@ -30,10 +30,10 @@ int TRANSFER = 0; // include transfer time in the benchmark?
void init() { void init() {
gaugeParam.X[0] = 64; gaugeParam.X[0] = 24;
gaugeParam.X[1] = 24; gaugeParam.X[1] = 24;
gaugeParam.X[2] = 26; gaugeParam.X[2] = 24;
gaugeParam.X[3] = 26; gaugeParam.X[3] = 32;
setDims(gaugeParam.X); setDims(gaugeParam.X);
......
...@@ -22,34 +22,35 @@ int main(int argc, char **argv) ...@@ -22,34 +22,35 @@ int main(int argc, char **argv)
setDims(Gauge_param.X); setDims(Gauge_param.X);
Gauge_param.blockDim = 64; Gauge_param.blockDim = 64;
Gauge_param.blockDim_sloppy = 64;
Gauge_param.cpu_prec = QUDA_DOUBLE_PRECISION; Gauge_param.cpu_prec = QUDA_DOUBLE_PRECISION;
Gauge_param.cuda_prec = QUDA_SINGLE_PRECISION; Gauge_param.cuda_prec = QUDA_DOUBLE_PRECISION;
Gauge_param.reconstruct = QUDA_RECONSTRUCT_12; Gauge_param.reconstruct = QUDA_RECONSTRUCT_12;
Gauge_param.cuda_prec_sloppy = QUDA_HALF_PRECISION; Gauge_param.cuda_prec_sloppy = QUDA_SINGLE_PRECISION;
Gauge_param.reconstruct_sloppy = QUDA_RECONSTRUCT_8; Gauge_param.reconstruct_sloppy = QUDA_RECONSTRUCT_12;
Gauge_param.gauge_fix = QUDA_GAUGE_FIXED_YES; Gauge_param.gauge_fix = QUDA_GAUGE_FIXED_NO;
Gauge_param.anisotropy = 1.0; Gauge_param.anisotropy = 1.0;
inv_param.inv_type = QUDA_CG_INVERTER; inv_param.inv_type = QUDA_BICGSTAB_INVERTER;
Gauge_param.t_boundary = QUDA_ANTI_PERIODIC_T; Gauge_param.t_boundary = QUDA_ANTI_PERIODIC_T;
Gauge_param.gauge_order = QUDA_QDP_GAUGE_ORDER; Gauge_param.gauge_order = QUDA_QDP_GAUGE_ORDER;
gauge_param = &Gauge_param; gauge_param = &Gauge_param;
double mass = -0.97; double mass = -0.95;
inv_param.kappa = 1.0 / (2.0*(4 + mass)); inv_param.kappa = 1.0 / (2.0*(4 + mass));
inv_param.tol = 1e-12; inv_param.tol = 1e-12;
inv_param.maxiter = 100; inv_param.maxiter = 10000;
inv_param.reliable_delta = 1e-8; inv_param.reliable_delta = 1e-3;
inv_param.mass_normalization = QUDA_KAPPA_NORMALIZATION; inv_param.mass_normalization = QUDA_KAPPA_NORMALIZATION;
inv_param.cpu_prec = QUDA_DOUBLE_PRECISION; inv_param.cpu_prec = QUDA_DOUBLE_PRECISION;
inv_param.cuda_prec = QUDA_SINGLE_PRECISION; inv_param.cuda_prec = QUDA_DOUBLE_PRECISION;
inv_param.cuda_prec_sloppy = QUDA_HALF_PRECISION; inv_param.cuda_prec_sloppy = QUDA_SINGLE_PRECISION;
inv_param.solution_type = QUDA_MAT_SOLUTION; inv_param.solution_type = QUDA_MAT_SOLUTION;
inv_param.matpc_type = QUDA_MATPC_EVEN_EVEN; inv_param.matpc_type = QUDA_MATPC_EVEN_EVEN;
inv_param.preserve_source = QUDA_PRESERVE_SOURCE_YES; // preservation doesn't work with reliable? inv_param.preserve_source = QUDA_PRESERVE_SOURCE_YES; // preservation doesn't work with reliable?
......
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