Decide on default behaviour for nonstandard input to discretization function
Despite the fact that main input to discretizer will be only the hamiltonian like
hamiltonian = kx * A * kx
and argument passed to discretization function will be multiplied by Psi(x) by us
object_to_discretize = kx * A * kx * Psi
I want to decide what should be output of discretization procedure when passed following input
C = 5.0 # some constant goes here, whatever it is
object_to_discretize = C*kx
To be honest I don't know what that input would even mean here so it's hard to decide on default behavior.