*** empty log message ***
options and defaults
QR maximum iteration: 30 for double, 50 for quad
cout replaces printf for quadmath overloading
#define in namespace nominally
Use JX Zhou 2D formulation
Elec_PK2.cpp was missing for some reason....