use cubes loop and unrolled quadrature

operators.toml

...s/curl_curl_plus_mass/N1E1ElementwiseCurlCurlPlusMass.hpp

...ss/avx/N1E1ElementwiseCurlCurlPlusMass_apply_macro_3D.cpp

...lusMass_computeInverseDiagonalOperatorValues_macro_3D.cpp

...noarch/N1E1ElementwiseCurlCurlPlusMass_apply_macro_3D.cpp

...lusMass_computeInverseDiagonalOperatorValues_macro_3D.cpp

operators/div_k_grad/CMakeLists.txt

operators/div_k_grad/P1ElementwiseDivKGrad.hpp

...s/div_k_grad/avx/P1ElementwiseDivKGrad_apply_macro_2D.cpp

...s/div_k_grad/avx/P1ElementwiseDivKGrad_apply_macro_3D.cpp

...ivKGrad_computeInverseDiagonalOperatorValues_macro_2D.cpp

...ivKGrad_computeInverseDiagonalOperatorValues_macro_3D.cpp

...iv_k_grad/noarch/P1ElementwiseDivKGrad_apply_macro_2D.cpp

...iv_k_grad/noarch/P1ElementwiseDivKGrad_apply_macro_3D.cpp

...ivKGrad_computeInverseDiagonalOperatorValues_macro_2D.cpp

...ivKGrad_computeInverseDiagonalOperatorValues_macro_3D.cpp

...k_grad/noarch/P1ElementwiseDivKGrad_toMatrix_macro_2D.cpp

...k_grad/noarch/P1ElementwiseDivKGrad_toMatrix_macro_3D.cpp