Add derivatives of shape function to ShapeFunctionMatrices() #796

Job	Run time
Correctness-MacOS (g++, Release, all)	19m 22s
Correctness-MacOS (g++, Release, none)	23m 13s
Correctness-MacOS (g++, Debug, all)	28m 17s
Correctness-MacOS (g++, Debug, none)	23m 57s
Correctness-MacOS (clang++, Release, all)	28m 0s
Correctness-MacOS (clang++, Release, none)	25m 30s
Correctness-MacOS (clang++, Debug, all)	24m 23s
Correctness-MacOS (clang++, Debug, none)	23m 32s
	3h 16m 14s

Provide feedback