add check for matrices
add planing for matrices
We do not need fft_3d_r2c_NxN_lower_then_copy.F anymore.
call ifc_matrix(n_atoms, Lx, Ly, Lz, c)
and write a 3x3 matrix at the Gamma point.
done: short-range interaction for sc