Bridge++  Ver. 1.1.x
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
fprop_4spinor_lex.cpp
Go to the documentation of this file.
1 
14 #include "fprop_4spinor_lex.h"
15 
16 //====================================================================
17 void Fprop_4spinor_lex::calc(Field_F *xq, const Field_F& b)
18 {
19  int Niter = 500;
20  double Stop_cond = 1.0e-24;
21 
22  int Nconv;
23  double diff;
24 
25  // CGNE solver
26  Field_F *b2 = new Field_F;
27 
28  m_fopr->set_mode("Ddag");
29  *b2 = (Field_F)m_fopr->mult(b);
30 
31  Solver *solver = new Solver_CG(m_fopr);
32  solver->set_parameters(Niter, Stop_cond);
33  m_fopr->set_mode("DdagD");
34 
35  // solver->solve((Field)(*xq),(Field)(*b2),Nconv,diff);
36  solver->solve(*xq, *b2, Nconv, diff);
37 
38  vout.general(m_vl, " Nconv = %d\n", Nconv);
39  vout.general(m_vl, " diff = %.8e\n", diff);
40 
41  delete solver;
42  delete b2;
43 }