#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-t...