Intermediary interface layer between high-level routines for physics-related algorithms and low-level routines dealing with matrix storage and manipulation