Advanced Computing Platform for Theoretical Physics

commit大文件会使得服务器变得不稳定,请大家尽量只commit代码,不要commit大的文件。

Commit 09e43a90 authored by mikeaclark's avatar mikeaclark
Browse files

Changed quda.h -> quda_internal

git-svn-id: http://lattice.bu.edu/qcdalg/cuda/quda@537 be54200a-260c-0410-bdd7-ce6af2a381ab
parent 14d68f3e
#include <stdlib.h>
#include <stdio.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#define REDUCE_THREADS 128
......
#include <quda.h>
#include <quda_internal.h>
#include <blas_reference.h>
// performs the operation x[i] *= a
......
#ifndef _QUDA_BLAS_REF_H
#define _QUDA_BLAS_REF_H
#include <quda.h>
#include <quda_internal.h>
#ifdef __cplusplus
extern "C" {
......
#include <stdio.h>
#include <stdlib.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <spinor_quda.h>
#include <blas_quda.h>
......
......@@ -2,11 +2,10 @@
#include <stdio.h>
#include <math.h>
#include <quda.h>
#include <quda_internal.h>
#include <spinor_quda.h>
#include <util_quda.h>
void allocateParityClover(ParityClover *ret, int *X, Precision precision)
{
ret->precision = precision;
......
......@@ -3,7 +3,7 @@
#include <cuComplex.h>
#include <quda.h>
#include <quda_internal.h>
#define gaugeSiteSize 18 // real numbers per link
#define spinorSiteSize 24 // real numbers per spinor
......
......@@ -2,7 +2,7 @@
#include <stdlib.h>
#include <math.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <dslash_reference.h>
......
#include <stdio.h>
#include <stdlib.h>
#include <quda.h>
#include <quda_internal.h>
#include <dslash_reference.h>
#include <util_quda.h>
#include <spinor_quda.h>
......@@ -64,7 +64,7 @@ void init() {
inv_param.cpu_prec = QUDA_DOUBLE_PRECISION;
inv_param.cuda_prec = QUDA_SINGLE_PRECISION;
inv_param.sp_pad = 24*24*24;
inv_param.sp_pad = 24*24*12;
if (test_type == 2) inv_param.dirac_order = QUDA_DIRAC_ORDER;
else inv_param.dirac_order = QUDA_DIRAC_ORDER;
......@@ -133,7 +133,7 @@ void init() {
gauge = cudaGaugePrecise;
if (clover_yes) {
loadCloverQuda(hostClover, hostCloverInv, &inv_param);
loadCloverQuda(NULL, hostCloverInv, &inv_param);
clover = cudaCloverPrecise;
cloverInv = cudaCloverInvPrecise;
}
......
#include <stdlib.h>
#include <stdio.h>
#include <quda.h>
#include <quda_internal.h>
#include <gauge_quda.h>
#include <xmmintrin.h>
......
......@@ -3,7 +3,7 @@
#include <math.h>
#include <cuComplex.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <spinor_quda.h>
#include <gauge_quda.h>
......@@ -124,10 +124,12 @@ void invertBiCGstabCuda(ParitySpinor x, ParitySpinor b, ParitySpinor r,
rNorm = sqrt(r2);
if (rNorm > maxrx) maxrx = rNorm;
if (rNorm > maxrr) maxrr = rNorm;
int updateR = (rNorm < delta*maxrr && r0Norm <= maxrr) ? 1 : 0;
int updateX = (rNorm < delta*r0Norm && r0Norm <= maxrx) ? 1 : 0;
//int updateR = (rNorm < delta*maxrr && r0Norm <= maxrr) ? 1 : 0;
//int updateX = (rNorm < delta*r0Norm && r0Norm <= maxrx) ? 1 : 0;
if (updateR || updateX) {
int updateR = (rNorm < delta*maxrr) ? 1 : 0;
if (updateR) {
if (x.precision != x_sloppy.precision) copyCuda(x, x_sloppy);
xpyCuda(x, y);
......
......@@ -2,7 +2,7 @@
#include <stdlib.h>
#include <math.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <spinor_quda.h>
#include <gauge_quda.h>
......
......@@ -4,7 +4,7 @@
#include <cuda_runtime.h>
#include <invert_quda.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <spinor_quda.h>
#include <gauge_quda.h>
......
......@@ -2,7 +2,7 @@
#include <stdio.h>
#include <time.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <dslash_reference.h>
......@@ -41,13 +41,13 @@ int main(int argc, char **argv)
} else {
inv_param.dslash_type = QUDA_WILSON_DSLASH;
}
inv_param.inv_type = QUDA_BICGSTAB_INVERTER;
inv_param.inv_type = QUDA_CG_INVERTER;
double mass = -0.94;
inv_param.kappa = 1.0 / (2.0*(1 + 3/gauge_param->anisotropy + mass));
inv_param.tol = 5e-7;
inv_param.maxiter = 10000;
inv_param.reliable_delta = 1e-1;
inv_param.reliable_delta = 1e-3;
inv_param.matpc_type = QUDA_MATPC_EVEN_EVEN;
inv_param.solution_type = QUDA_MAT_SOLUTION;
......@@ -59,7 +59,7 @@ int main(int argc, char **argv)
inv_param.preserve_source = QUDA_PRESERVE_SOURCE_YES;
inv_param.dirac_order = QUDA_DIRAC_ORDER;
inv_param.sp_pad = 0;
inv_param.sp_pad = 24*24*24;
if (clover_yes) {
inv_param.clover_cpu_prec = QUDA_DOUBLE_PRECISION;
......
#include <stdio.h>
#include <stdlib.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <gauge_quda.h>
......
#include <stdlib.h>
#include <stdio.h>
#include <quda.h>
#include <quda_internal.h>
#include <spinor_quda.h>
#include <util_quda.h>
......
#include <stdio.h>
#include <stdlib.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <gauge_quda.h>
......
......@@ -2,7 +2,7 @@
#include <stdio.h>
#include <sys/time.h>
#include <quda.h>
#include <quda_internal.h>
#include <util_quda.h>
#include <dslash_reference.h>
......
#ifndef _QUDA_UTIL_H
#define _QUDA_UTIL_H
#include <quda.h>
#include <quda_internal.h>
#ifdef __cplusplus
extern "C" {
......
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