... efficient structure-preserving algorithms are developed by combing the Crank–Nicolson (C-N) scheme with the extrapolation method, and the numerical computation is implemented by solving a linear system. Findings The proposed algorithms can enhance the computation efficiency and satisfy some physical...