added a part to compute error
move basic parameter reading to problem.h
control of computing error or not
init with initFEMProblem()