I take exception to what you wrote about SPMCHF.
The test example showed that, if a user seletect the 'all' parameters for both level 1 and level 2 and the default accuracy for each level as specified, that the Virial theorem changed as shown:
EIV : Ratio -2.000000311943682
NR : Ratio -1.999999999999929
which is essential machine precision for NR.
This code is written more in Oleg's style with default parameters. The user should be able to select only level 1 (I will check that) The parameters are always printed for each level. But you are right -- the "intelligent control' by the code needs to be tested further. But actually, I think the user should specify the defaults at compile time and then modify with command-line arguments the way Oleg always did.
I take exception to what you wrote about SPMCHF.
The test example showed that, if a user seletect the 'all' parameters for both level 1 and level 2 and the default accuracy for each level as specified, that the Virial theorem changed as shown:
EIV : Ratio -2.000000311943682
NR : Ratio -1.999999999999929
which is essential machine precision for NR.
This code is written more in Oleg's style with default parameters. The user should be able to select only level 1 (I will check that) The parameters are always printed for each level. But you are right -- the "intelligent control' by the code needs to be tested further. But actually, I think the user should specify the defaults at compile time and then modify with command-line arguments the way Oleg always did.