FEM now makes local element stiffness matrix

in the process it constructs all the reference basis functions