smooth initial density and magnetization
new Matching_coefficients class