Essentially fixed this today, at least for hamilt.coupling.n_eqs=1, which is quite nice already. And it fixes the MolVib/H3+ demo.

However, for hamilt.coupling.n_eqs>1 the imaginary time propagator (cheby_imag) is still throwing an error message. But this shouldn't be too difficult to fix ... one day ;-)

On that occasion, the scalar products could be calculated using
coeff = wave.braket ( psi.bound{ii}, psi.dvr{1} );

Perhaps something similar for the normalization ?!?

 

Last edit: Burkhard Schmidt 2019-08-26