The default dtype of macehamiltonian is float32, this is hardcoded in the parameters method of MACEHamiltonian.
For certain methods, such as optimize, this default is overwritten to float64, but at the moment you cannot change the dtype to compute energy, forces, and/or stresses via the compute method of MACEHamiltonian.
Ideally the user can set dtype when initializing the MACEHamiltonian. However, when set, should we still overwrite the dtype for methods such as optimize?