Implement just rhs assembly
Implement assemble_3d_coo_A()
Precalculate the matrix A in CSR format
Use static arrays
Precalculate the CSR matrix
Move derivatives together with energies
Merge derivative/value free energy calculation