RunComputeProfiles and RunVIP concurrence

[stephane-caron]

In the current version,

ret = x.RunComputeProfiles(0, 1e-4)
print "RunComputeProfiles:", ret
ret = x.RunVIP(sdbeg_min, sdbeg_max)
print "RunVIP:", ret

Yields incorrect VIP results (because RunComputeProfiles affects x in a way that RunVIP does not detect). We should make the two independent in the future.

[quangounet]

Could you provide an example?

[stephane-caron]

Sure:

discrtimestep = 1e-4
sdbeg_min, sdbeg_max = 0., 1e-4
trajectorystring = """1.000000
2
1.02380743853 -0.621211430375 6.66133814775e-16 -4.4408920985e-16
-0.294007923442 -0.228798693785 1.66533453694e-16 0.0"""

When I run it with only RunVIP, I get:

sdendmin = 4.79738632613
sdendmax = 9.42902801911

If there is a call to RunComputeProfiles before, it becomes:

sdendmin = 0.0095855458967
sdendmax = 0.306737468694

[quangounet]

I reset resprofileslist at the beginning of Preprocess and it solved the bug on this instance.