@Kim.Liegeois I need your help to adapt build_thermal and build_thermomecanical in Tetra4 and Hex8 classes. I would need at least the formula used to compute the matrices. I could reconstruct it be "de-rolling" the loops, but I ma lazy :). Also, your tests contain no value check, I therefore cannot be sure of what I am doing...