[ET Trac] #2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D

[Gabriele Bozzola]

#2741: Some thorns seem to incorrectly use CCTK_GFINDEX4D

 Reporter: Gabriele Bozzola
   Status: open
Milestone: 
  Version: 
     Type: bug
 Priority: major
Component: EinsteinToolkit thorn

Comment (by Gabriele Bozzola):

Yes, we were running an older version and sorry for opening a duplicated ticket.  The problem is now only in `wvuthorns_diagnostics`

```
particle_tracerET/src/compute_particles_four_velocities.C:          const double vU[3] = { vel[CCTK_GFINDEX4D(cctkGH,i,j,k,0)],
particle_tracerET/src/compute_particles_four_velocities.C:                                 vel[CCTK_GFINDEX4D(cctkGH,i,j,k,1)],
particle_tracerET/src/compute_particles_four_velocities.C:                                 vel[CCTK_GFINDEX4D(cctkGH,i,j,k,2)] };
particle_tracerET/src/convert_to_MHD_3velocity.C:       double ETvx = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,0)];
particle_tracerET/src/convert_to_MHD_3velocity.C:       double ETvy = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,1)];
particle_tracerET/src/convert_to_MHD_3velocity.C:       double ETvz = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,2)];
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = -(y[index] + 0.5*dY)*A_b*pow(max(Pressure_at_Ax_stagger-P_cut,0.0),n_s);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          if(!have_two_NSs_along_x_axis) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] =  (xL + 0.5*dX)*A_b*pow(max(Pressure_at_Ay_stagger-P_cut,0.0),n_s);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:            if(r1<=r_NS1) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] =  (x1 + 0.5*dX)*A_b*pow(max(Pressure_at_Ay_stagger-P_cut,0.0),n_s);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:            if(r2<=r_NS2) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] =  (x2 + 0.5*dX)*A_b*pow(max(Pressure_at_Ay_stagger-P_cut,0.0),n_s);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:            if(r1>r_NS1 && r2>r_NS2) Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = 0.0; // No external B-field.
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = - (yL + 0.5 * dY) * (Ap1 + Ap2);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] = (x1 + 0.5 * dX) * Ap1 + (x2 + 0.5 * dX) * Ap2;
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,2)] =  0.0;
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = -y[index]*A_b*pow(max(PL-P_cut,0.0),n_s);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)] =  x[index]*A_b*pow(max(PL-P_cut,0.0),n_s);
Seed_Magnetic_Fields_BNS/src/Seed_Magnetic_Fields.C:          Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,2)] =  0.0;
smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C:   double ETvx = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,0)];
smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C:   double ETvy = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,1)];
smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C:   double ETvz = vel[CCTK_GFINDEX4D(cctkGH,i,j,k,2)];
smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C:   double Bx_center = Bvec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)];
smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C:   double By_center = Bvec[CCTK_GFINDEX4D(cctkGH,i,j,k,1)];
smallbPoynET/src/compute_bi_b2_Poyn_fluxET.C:   double Bz_center = Bvec[CCTK_GFINDEX4D(cctkGH,i,j,k,2)];
```

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2741/some-thorns-seem-to-incorrectly-use
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.einsteintoolkit.org/pipermail/trac/attachments/20230609/3f6235fc/attachment.htm>