**Acknowledgements**

This work was financially supported by the Science and Technology Development Plan of Jilin Province (Nos. 20200201057JC and 20190201108JC), and the Technology Research Planning Project of Education Department of Jilin Province (No. JJKH20220690KJ).

The authors would like to thank K. Svanberg for providing the Matlab code of MMA optimizer.
