Open EranOfek opened 4 weeks ago
In some cases celestial.SolarSys.orbitIntegration
takes very long time to finish. This happens when rk1210 is used an the object passes near one of the perturbers. For example: running:
[UpdatedCat, OnlyMP(If), AstrometricCat, Status(If), Path] = VO.ZTF.reanalize_ZTF_psfcat(Files(If).name, 'Path',[], 'OrbEl',Args.OrbEl, 'INPOP',Args.INPOP, 'GeoPos',Args.GeoPos,... 'WriteProd',true,... 'FieldID',Args.FieldID,... 'CCDID',Args.CCDID); on: /raid/eran/projects/telescopes/ZTF/sourceCatalogs/DAO/626/1 ztf_20180806177303_000626_zg_c01_o_q2_psfcat.fits
This does not happen with ode45. Likely due to a problem with rkn1210. If can't be fixed, then, need to modify that will return NaN in finite time and in these cases ode45 will be used
In some cases celestial.SolarSys.orbitIntegration
takes very long time to finish. This happens when rk1210 is used an the object passes near one of the perturbers. For example: running:
[UpdatedCat, OnlyMP(If), AstrometricCat, Status(If), Path] = VO.ZTF.reanalize_ZTF_psfcat(Files(If).name, 'Path',[], 'OrbEl',Args.OrbEl, 'INPOP',Args.INPOP, 'GeoPos',Args.GeoPos,... 'WriteProd',true,... 'FieldID',Args.FieldID,... 'CCDID',Args.CCDID); on: /raid/eran/projects/telescopes/ZTF/sourceCatalogs/DAO/626/1 ztf_20180806177303_000626_zg_c01_o_q2_psfcat.fits
This does not happen with ode45. Likely due to a problem with rkn1210. If can't be fixed, then, need to modify that will return NaN in finite time and in these cases ode45 will be used