Hi Roland,
After looking at this section of code further, it seems that the termination criterion for the bare mass search could be improved. Attached is an updated patch which implements and improved check. It now continues trying to find the bare masses until the ADM masses are within a specified tolerance of their target. Other than the change in termination criterion, the algorithm is unchanged, although I have simplified the calculation of a new bare mass guess somewhat.
Barry
On 14/05/2010 02:52, Roland Haas wrote:
Hello Barry,
Part of this patch removes the explicit cast to double in the assignment to new_mass as you suggest. I've just noticed that there is a similar cast in the assignment to old_mass which I think should also be removed.
I don't see, to see this second patch. Did you forget to attach it? I'll commit the patch to TwoPunctures once I receive it unless I hear any objections.
Yours, Roland