add some omp directives in the cal_phys.c and odefun.c