use C++ standard headers
*** 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