The work of K. Zhang was supported by the NSF of China under the grant Nos. 11871245, 11771179, 11726102 and by the Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, Jilin University (No. 93K172018Z01). The work of H. Song was supported by the NSF of China under the grant No. 11701210, the Education Department Project of Jilin Province under the grant No. JJKH20180113KJ, the Science and Technology Department Project of Jilin Province under the grant No. 20190103029JH, and the special funds for basic scientific research operating expenses of central universities. They also wish to thank the high performance computing center of Jilin university and computing center of Jilin province for essential computing support.