add rokko/utility/magnetization.hpp