Advanced Computing Platform for Theoretical Physics

Commit 0c0d5e54 authored by rbabich's avatar rbabich
Browse files

forgot to change matpc_type for the reference dslash. Asymmetric clover

seems to be working now.


git-svn-id: http://lattice.bu.edu/qcdalg/cuda/quda@469 be54200a-260c-0410-bdd7-ce6af2a381ab
parent 3d8d0d54
...@@ -61,7 +61,7 @@ void init() { ...@@ -61,7 +61,7 @@ void init() {
inv_param.kappa = kappa; inv_param.kappa = kappa;
inv_param.matpc_type = QUDA_MATPC_EVEN_EVEN_ASYMMETRIC; inv_param.matpc_type = QUDA_MATPC_ODD_ODD_ASYMMETRIC;
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_SINGLE_PRECISION;
...@@ -233,7 +233,14 @@ double dslashCUDA() { ...@@ -233,7 +233,14 @@ double dslashCUDA() {
} }
void dslashRef() { void dslashRef() {
// to be removed once reference clover is finished
if (inv_param.matpc_type == QUDA_MATPC_EVEN_EVEN_ASYMMETRIC) {
inv_param.matpc_type = QUDA_MATPC_EVEN_EVEN;
} else if (inv_param.matpc_type == QUDA_MATPC_ODD_ODD_ASYMMETRIC) {
inv_param.matpc_type = QUDA_MATPC_ODD_ODD;
}
// compare to dslash reference implementation // compare to dslash reference implementation
printf("Calculating reference implementation..."); printf("Calculating reference implementation...");
fflush(stdout); fflush(stdout);
...@@ -243,7 +250,7 @@ void dslashRef() { ...@@ -243,7 +250,7 @@ void dslashRef() {
inv_param.cpu_prec, gaugeParam.cpu_prec); inv_param.cpu_prec, gaugeParam.cpu_prec);
break; break;
case 1: case 1:
matpc(spinorRef, hostGauge, spinorEven, kappa, QUDA_MATPC_EVEN_EVEN, DAGGER_BIT, matpc(spinorRef, hostGauge, spinorEven, kappa, inv_param.matpc_type, DAGGER_BIT,
inv_param.cpu_prec, gaugeParam.cpu_prec); inv_param.cpu_prec, gaugeParam.cpu_prec);
break; break;
case 2: case 2:
......
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