reinitialization must be parallelized
refine on phi