This work was supported by the U.S. Department of Energy contract DE-AC05-06OR23177, under which Jefferson Science Associates, LLC operates Jefferson Lab, and also by the U.S. Department of Energy Grant Nos. DE-FG02-87ER40365 and DE-FG02-92ER40735, by U.S. National Science Foundation Grant No. PHY-2209183, by the Spanish Ministerio de Ciencia e Innovacion (MICINN) Grant Nos. PID2020-112777GBI00, PID2019–106080GB-C21, and PID2020-118758GB-I00. DW is supported by National Natural Science Foundation of China Grant No. 12035007 and the NSFC and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through the funds provided to the Sino-German Collaborative Research Center TRR110 ``Symmetries and the Emergence of Structure in QCD'' (NSFC Grant No. 12070131001, DFG ProjectID 196253076-TRR 110). CFR is supported by Spanish Ministerio de Educacion y Formacion Profesional under Grant No. BG20/00133. ANHB is supported by the DFG through the Research Unit FOR 2926 (project number 409651613). MA is supported by Generalitat Valenciana under Grant No. CIDEGENT/2020/002. NH is supported by a Polish research project with no. 2018/29/B/ST2/02576 (National Science Center). VM is a Serra Hunter fellow. VS acknowledges the support of the USDOE ExoHad Topical Collaboration, contract DE-SC0023598. This research was supported in part by Lilly Endowment, Inc., through its support for the Indiana University Pervasive Technology Institute. We acknowledge the computational resources and assistance provided by the Centro de Computacion de Alto Rendimiento CCAR-UNED. This work contributes to the aims of the U.S. Department of Energy ExoHad Topical Collaboration, contract DE-SC0023598.

This version of the code corresponds to the published version of the paper arXiv:2305.01449