added check for gcr inverter

checks/check.cpp

@@ -15,12 +15,8 @@ void run_checks()
 	{
 		lprintf("CHECK", INFO, "Status: PASSED (all checks succeded)");
 	}
-	else if(failed == 1)
-	{
-		lprintf("CHECK", INFO, "Status: FAILED (%d check reported an error)", failed);
-	}
 	else
 	{
-		lprintf("CHECK", INFO, "Status: FAILED (%d checks reported an error)", failed);
+		lprintf("CHECK", INFO, "Status: FAILED (%d check(s) reported an error)", failed);
 	}
 }

checks/inverters.cpp

@@ -20,6 +20,7 @@ static int check_bicgstab()
 	spinor_allocate(Mu);
 
 	// Random spinor fields
+	spinor_zero(Mu);
 	spinor_random(u);
 	norm = spinor_sqnorm(u);
 
@@ -61,6 +62,7 @@ static int check_minres()
 	spinor_allocate(Mu);
 
 	// Random spinor fields
+	spinor_zero(Mu);
 	spinor_random(u);
 	norm = spinor_sqnorm(u);
 
@@ -89,6 +91,48 @@ static int check_minres()
 	}
 }
 
+static int check_gcr()
+{
+	SpinorField u, Mu, v;
+	double r1, norm;
+	inv_par par;
+	int it;
+
+	// Allocate memory
+	spinor_allocate(u);
+	spinor_allocate(v);
+	spinor_allocate(Mu);
+
+	// Random spinor fields
+	spinor_zero(Mu);
+	spinor_random(u);
+	norm = spinor_sqnorm(u);
+
+	// Apply dirac operator
+	par.prec = prec;
+	par.mass = mass;
+	par.log = 0;
+	par.mvm = Dphi;
+	it = gcr_inverter(par, Mu, u);
+
+	// Check result
+	Dphi(mass, v, Mu);
+	spinor_sub_assign(v,u);
+	r1 = spinor_sqnorm(v) / norm;
+
+	// Evaluate result
+	if(r1 < prec)
+	{
+		lprintf("CHECK", INFO, "GCR inverter: PASSED [diff = %1.6e, iter = %d]", r1, it);
+		return 0;
+	}
+	else
+	{
+		lprintf("CHECK", INFO, "GCR inverter: FAILED [diff = %1.6e, iter = %d]", r1, it);
+		return 1;
+	}
+}
+
 static int check_cg()
 {
 	SpinorField u, Mu, v;
@@ -102,6 +146,7 @@ static int check_cg()
 	spinor_allocate(Mu);
 
 	// Random spinor fields
+	spinor_zero(Mu);
 	spinor_random(u);
 	norm = spinor_sqnorm(u);
 
@@ -195,6 +240,7 @@ int check_inverters()
 	int ret = 0;
 	ret += check_bicgstab();
 	ret += check_minres();
+	ret += check_gcr();
 	ret += check_cg();
 	ret += check_cg_mshift();
 	return ret;