#1746: Index error in TwoPunctures function JFD_times_dv() -----------------------------------+---------------------------------------- Reporter: ixr5289@… | Owner: Type: defect | Status: new Priority: minor | Milestone: Component: EinsteinToolkit thorn | Version: development version Keywords: TwoPunctures | -----------------------------------+---------------------------------------- In the TwoPunctures thorn, in the file FuncAndJacobian.c, the function JDF_times_dv() contains an index error. The (A,B,phi) coordinates are defined in terms of indices (i,j,k). The current version has
phi = hp * j;
The j index is used to define the B coordinate. The statement should be
phi = hp * k;