39 set_threadtask(ith, nth, is, ns,
m_Nvol);
41 for (
int mu = 0; mu <
m_Nex; ++mu) {
42 for (
int site = is; site < ns; ++site) {
57 set_threadtask(ith, nth, is, ns,
m_Nvol);
59 for (
int mu = 0; mu <
m_Nex; ++mu) {
60 for (
int site = is; site < ns; ++site) {
61 this->
mat(ut, site, mu);
75 set_threadtask(ith, nth, is, ns,
m_Nvol);
77 for (
int mu = 0; mu <
m_Nex; ++mu) {
78 for (
int site = is; site < ns; ++site) {
89 const Field_G& U1,
const int ex1,
90 const Field_G& U2,
const int ex2)
93 assert(ex1 < U1.
nex());
94 assert(ex2 < U2.
nex());
98 const int Nvol = W.
nvol();
100 int ith, nth, is, ns;
101 set_threadtask(ith, nth, is, ns, Nvol);
103 for (
int site = is; site < ns; ++site) {
111 const Field_G& U1,
const int ex1,
112 const Field_G& U2,
const int ex2)
114 assert(ex < W.
nex());
115 assert(ex1 < U1.
nex());
116 assert(ex2 < U2.
nex());
120 const int Nvol = W.
nvol();
122 int ith, nth, is, ns;
123 set_threadtask(ith, nth, is, ns, Nvol);
125 for (
int site = is; site < ns; ++site) {
133 const Field_G& U1,
const int ex1,
134 const Field_G& U2,
const int ex2)
136 assert(ex < W.
nex());
137 assert(ex1 < U1.
nex());
138 assert(ex2 < U2.
nex());
142 const int Nvol = W.
nvol();
144 int ith, nth, is, ns;
145 set_threadtask(ith, nth, is, ns, Nvol);
147 for (
int site = is; site < ns; ++site) {
155 const Field_G& U1,
const int ex1,
156 const Field_G& U2,
const int ex2)
158 assert(ex < W.
nex());
159 assert(ex1 < U1.
nex());
160 assert(ex2 < U2.
nex());
164 const int Nvol = W.
nvol();
166 int ith, nth, is, ns;
167 set_threadtask(ith, nth, is, ns, Nvol);
169 for (
int site = is; site < ns; ++site) {
177 const Field_G& U1,
const int ex1,
178 const Field_G& U2,
const int ex2,
181 assert(ex < W.
nex());
182 assert(ex1 < U1.
nex());
183 assert(ex2 < U2.
nex());
187 const int Nvol = W.
nvol();
189 int ith, nth, is, ns;
190 set_threadtask(ith, nth, is, ns, Nvol);
192 for (
int site = is; site < ns; ++site) {
193 W.
add_mat(site, ex, U1.
mat(site, ex1) * U2.
mat(site, ex2) * ff);
200 const Field_G& U1,
const int ex1,
201 const Field_G& U2,
const int ex2,
204 assert(ex < W.
nex());
205 assert(ex1 < U1.
nex());
206 assert(ex2 < U2.
nex());
210 const int Nvol = W.
nvol();
212 int ith, nth, is, ns;
213 set_threadtask(ith, nth, is, ns, Nvol);
215 for (
int site = is; site < ns; ++site) {
223 const Field_G& U1,
const int ex1,
224 const Field_G& U2,
const int ex2,
227 assert(ex < W.
nex());
228 assert(ex1 < U1.
nex());
229 assert(ex2 < U2.
nex());
233 const int Nvol = W.
nvol();
235 int ith, nth, is, ns;
236 set_threadtask(ith, nth, is, ns, Nvol);
238 for (
int site = is; site < ns; ++site) {
246 const Field_G& U1,
const int ex1,
247 const Field_G& U2,
const int ex2,
250 assert(ex < W.
nex());
251 assert(ex1 < U1.
nex());
252 assert(ex2 < U2.
nex());
256 const int Nvol = W.
nvol();
258 int ith, nth, is, ns;
259 set_threadtask(ith, nth, is, ns, Nvol);
261 for (
int site = is; site < ns; ++site) {
270 assert(ex < W.
nex());
272 const int Nvol = W.
nvol();
274 int ith, nth, is, ns;
275 set_threadtask(ith, nth, is, ns, Nvol);
277 for (
int site = is; site < ns; ++site) {
286 assert(ex < W.
nex());
288 const int Nvol = W.
nvol();
290 int ith, nth, is, ns;
291 set_threadtask(ith, nth, is, ns, Nvol);
293 for (
int site = is; site < ns; ++site) {
301 const double alpha,
const Field_G& iP,
const Field_G& U,
const int Nprec)
305 const int Nvol = U.
nvol();
306 const int Nex = U.
nex();
311 int ith, nth, is, ns;
312 set_threadtask(ith, nth, is, ns, Nvol);
314 for (
int ex = 0; ex < Nex; ++ex) {
315 for (
int site = is; site < ns; ++site) {
316 u0 = U.
mat(site, ex);
317 v0 = iP.
mat(site, ex);