def save_int(self):
    self.fha.write(('\n  // PF integral, m=%d  \n' % 0))
    self.fha.write('  if(I == 1 && J == 3){ \n')
    self.fha.write('    PFint_0 pf(PAx, PAy, PAz, PBx, PBy, PBz, PCx, PCy, PCz, TwoZetaInv, store, YVerticalTemp); \n')
    self.fha.write('#ifdef REG_PF \n')
    for i in range(0, 10):
        for j in range(0, 3):
            self.fha.write(('    LOCSTORE(store, %d, %d, STOREDIM, STOREDIM) = pf.x_%d_%d;\n' % ((j + 1), (i + 10), (j + 1), (i + 10))))
    self.fha.write('#endif \n')
    if (OEint.debug == 1):
        self.fha.write('\n#ifdef DEBUG_OEI \n')
        for i in range(0, 10):
            for j in range(0, 3):
                self.fha.write(('    printf("II %%d JJ %%d %s store[%d,%d] = %%f \\n", II, JJ, LOCSTORE(store, %d, %d, STOREDIM, STOREDIM)); \n' % ('PF', (j + 1), (i + 10), (j + 1), (i + 10))))
        self.fha.write('#endif \n\n')
    self.fha.write('  } \n')
def save_int(self):
    self.fha.write(('\n  // PF integral, m=%d  \n' % 0))
    self.fha.write('  if(I == 1 && J == 3){ \n')
    self.fha.write('    PFint_0 pf(PAx, PAy, PAz, PBx, PBy, PBz, PCx, PCy, PCz, TwoZetaInv, store, YVerticalTemp); \n')
    self.fha.write('#ifdef REG_PF \n')
    for i in range(0, 10):
        for j in range(0, 3):
            self.fha.write(('    LOCSTORE(store, %d, %d, STOREDIM, STOREDIM) = pf.x_%d_%d;\n' % ((j + 1), (i + 10), (j + 1), (i + 10))))
    self.fha.write('#endif \n')
    if (OEint.debug == 1):
        self.fha.write('\n#ifdef DEBUG_OEI \n')
        for i in range(0, 10):
            for j in range(0, 3):
                self.fha.write(('    printf("II %%d JJ %%d %s store[%d,%d] = %%f \\n", II, JJ, LOCSTORE(store, %d, %d, STOREDIM, STOREDIM)); \n' % ('PF', (j + 1), (i + 10), (j + 1), (i + 10))))
        self.fha.write('#endif \n\n')
    self.fha.write('  } \n')