Beispiel #1
0
void random_jacobi_field(su3_vector * const k, const int V)
{
int ix;
su3_vector *s;
double v[6];

 for (ix=0; ix<V ;ix++)
   {
     s=k+ix;
     *s=unif_su3_vector();
   }
#ifdef TM_USE_MPI
 xchange_jacobi(k);
#endif
}
Beispiel #2
0
void random_gauss_jacobi_field(su3_vector * const k, const int V)
{
int ix;
su3_vector *s;
double v[6];

 for (ix=0; ix<V ;ix++) {
     s=k+ix;
     gauss_vector(v,6);
     s->c0 = v[0] + v[1] * I;
     s->c1 = v[2] + v[3] * I;
     s->c2 = v[4] + v[5] * I;
 }
#ifdef TM_USE_MPI
 xchange_jacobi(k);
#endif
}
Beispiel #3
0
void random_gauss_jacobi_field(su3_vector * const k, const int V)
{
int ix;
su3_vector *s;
double v[6];

 for (ix=0; ix<V ;ix++) {
     s=k+ix;
     gauss_vector(v,6);
     (*s).c0.re=v[0];
     (*s).c0.im=v[1];
     (*s).c1.re=v[2];
     (*s).c1.im=v[3];
     (*s).c2.re=v[4];
     (*s).c2.im=v[5];
 }
#ifdef MPI
 xchange_jacobi(k);
#endif
}