Add variable coefficients to matrix-free Helmholtz operator.
MR !1284 (merged) adds constant anisotropy to the Helmholtz operator. This should be extended with variable coefficients.
Note that varcoeffs
in ContField.cpp:417 is not fully populated.