move numerical variance computation to header file
compare contributions to numerical and analytical variance