Merge branch 'kohl/split-toml-and-ci-rebased' into 'main'

[Rebased] Splitting toml files and CI for better overview and more parallelism.

See merge request !16

.gitlab-ci.yml

-hyteg-hog-integration-test:
+stages:
+  - pretest
+  - test
+
+compare_number_of_directories_and_toml_files:
+  stage: pretest
+  image: i10git.cs.fau.de:5005/pycodegen/pycodegen/full
+  script:
+    - find operators/* -maxdepth 0 -type d | wc -l > num_directories.txt
+    - find ./*.toml | wc -l > num_toml_files.txt
+    - cat num_directories.txt num_toml_files.txt
+    - diff num_directories.txt num_toml_files.txt
+  tags:
+    - docker
+
+.hyteg-hog-integration-test-template:
   stage: test
   image: i10git.cs.fau.de:5005/pycodegen/pycodegen/full
   script:
@@ -10,10 +25,66 @@ hyteg-hog-integration-test:
     - GIT_CONFIG_COUNT=1 GIT_CONFIG_KEY_0="url.https://gitlab-ci-token:${CI_JOB_TOKEN}@i10git.cs.fau.de.insteadOf" GIT_CONFIG_VALUE_0="ssh://git@i10git.cs.fau.de" python -m pip install -r requirements.txt
     - pip freeze
     - rm -rf ../operators/*
-    - python3 generate.py -o ../operators ../operators.toml --processes 8
+    - python3 generate.py -o ../operators ../operators-$OPERATOR_BASE_NAME.toml --processes 8
     - cd ..
+    - cd operators/$OPERATOR_BASE_NAME
+    - pwd
+    - ls -lhta
     - git status
-    - git diff
-    - git diff --exit-code --stat
+    - git diff .
+    - git diff --exit-code --stat .
+  variables:
+    - OPERATOR_BASE_NAME: "OPERATOR_NOT_SPECIFIED"
   tags:
     - docker
+
+curl_curl:
+  stage: pretest                                       # <=== THIS RUNS IN PRETEST STAGE
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "curl_curl"
+
+diffusion:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "diffusion"
+
+div_k_grad:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "div_k_grad"
+
+divergence:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "divergence"
+
+epsilon:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "epsilon"
+
+full_stokes:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "full_stokes"
+
+gradient:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "gradient"
+
+k_mass:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "k_mass"
+
+mass:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "mass"
+
+terraneo:
+  extends: .hyteg-hog-integration-test-template
+  variables:
+    OPERATOR_BASE_NAME: "terraneo"
\ No newline at end of file

README.md

@@ -5,7 +5,8 @@ The code under [operators](operators) is generated using [HFG](https://i10git.cs
 Each kernel has a platform independent and, where applicable, an AVX-vectorized version.
 The desired variant is configured in CMake.
 
-All operators are described declaratively in [operators.toml](operators.toml).
+All operators are described declaratively in `operators-<operator-group>.toml`.
+They are roughly split into groups for a better overview and faster integration tests.
 
 ## How to generate operators
 
@@ -16,5 +17,7 @@ python -m venv .venv
 source .venv/bin/activate
 python -m pip install -r requirements.txt
 
-python generate.py -o ../operators ../operators.toml
+python generate.py -o ../operators ../operators-<operator-group>.toml
+
+# To generate all operators you have to run the script repeatedly, each time with a different toml file.
 ```

operators-curl_curl.toml

+[[curl_curl]]
+trial-space   = "N1E1"
+test-space    = "N1E1"
+dimensions    = [3]
+quadrature    = 0
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]

operators-diffusion.toml

+[[diffusion]]
+trial-space   = "P1"
+test-space    = "P1"
+dimensions    = [2, 3]
+quadrature    = 0
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[diffusion]]
+trial-space   = "P1"
+test-space    = "P1"
+dimensions    = [2, 3]
+quadrature    = 0
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "ParametricP1Map"
+
+[[diffusion]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[diffusion]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "AnnulusMap"
+
+[[diffusion]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [3]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "IcosahedralShellMap"
+
+[[diffusion]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "ParametricP2Map"

operators-div_k_grad.toml

+[[div_k_grad]]
+trial-space   = "P1"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[div_k_grad]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[div_k_grad]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [2]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "AnnulusMap"
+
+[[div_k_grad]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [3]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "IcosahedralShellMap"
+
+[[div_k_grad]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+operator-name = "DivKGradP1Coefficient"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[div_k_grad]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2]
+quadrature    = 3
+loop-strategy = "sawtooth"
+operator-name = "DivKGradP1Coefficient"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "AnnulusMap"
+
+[[div_k_grad]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [3]
+quadrature    = 3
+loop-strategy = "sawtooth"
+operator-name = "DivKGradP1Coefficient"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "IcosahedralShellMap"
\ No newline at end of file

operators-divergence.toml

+[[divergence]]
+trial-space   = "P2"
+test-space    = "P1"
+components-trial = [0, 1, 2]
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[divergence]]
+trial-space   = "P2"
+test-space    = "P1"
+components-trial = [0, 1]
+dimensions    = [2]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "AnnulusMap"
+
+[[divergence]]
+trial-space   = "P2"
+test-space    = "P1"
+components-trial = [0, 1, 2]
+dimensions    = [3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "IcosahedralShellMap"
+
+[[divergence]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[divergence]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+dimensions    = [2]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "AnnulusMap"
+
+[[divergence]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+dimensions    = [3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "IcosahedralShellMap"

operators-epsilon.toml

+[[epsilon]]
+trial-space               = "P2"
+test-space                = "P2"
+components-trial = [0, 1, 2]
+components-test  = [0, 1, 2]
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2, 3]
+quadrature      = 2
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[epsilon]]
+trial-space               = "P2"
+test-space                = "P2"
+components-trial = [0, 1]
+components-test  = [0, 1]
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "AnnulusMap"
+
+[[epsilon]]
+trial-space               = "P2"
+test-space                = "P2"
+components-trial = [0, 1, 2]
+components-test  = [0, 1, 2]
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [3]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "IcosahedralShellMap"
+
+[[epsilon]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2, 3]
+quadrature      = 2
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[epsilon]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "AnnulusMap"
+
+[[epsilon]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [3]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "IcosahedralShellMap"
+
+[[epsilon]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P1"
+dimensions      = [2]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+operator-name   = "EpsilonP1Viscosity"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "AnnulusMap"
+
+[[epsilon]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P1"
+dimensions      = [3]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+operator-name   = "EpsilonP1Viscosity"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "IcosahedralShellMap"

operators-full_stokes.toml

+[[full_stokes]]
+trial-space               = "P2"
+test-space                = "P2"
+components-trial = [0, 1, 2]
+components-test  = [0, 1, 2]
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2, 3]
+quadrature      = 2
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[full_stokes]]
+trial-space               = "P2"
+test-space                = "P2"
+components-trial = [0, 1]
+components-test  = [0, 1]
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "AnnulusMap"
+
+[[full_stokes]]
+trial-space               = "P2"
+test-space                = "P2"
+components-trial = [0, 1, 2]
+components-test  = [0, 1, 2]
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [3]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "IcosahedralShellMap"
+
+[[full_stokes]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2, 3]
+quadrature      = 2
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[full_stokes]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [2]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "AnnulusMap"
+
+[[full_stokes]]
+trial-space               = "P2Vector"
+test-space                = "P2Vector"
+form-space-args.coefficient_function_space = "P2"
+dimensions      = [3]
+quadrature      = 3
+loop-strategy   = "sawtooth"
+optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending        = "IcosahedralShellMap"

operators-gradient.toml

+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2"
+components-test = [0, 1, 2]
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2"
+components-test = [0, 1]
+dimensions    = [2]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "AnnulusMap"
+
+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2"
+components-test = [0, 1, 2]
+dimensions    = [3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "IcosahedralShellMap"
+
+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2Vector"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "cubes"
+optimizations = ["moveconstants", "vectorize"]
+
+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2Vector"
+dimensions    = [2]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "AnnulusMap"
+
+[[gradient]]
+trial-space   = "P1"
+test-space    = "P2Vector"
+dimensions    = [3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize"]
+blending      = "IcosahedralShellMap"

operators-k_mass.toml

+[[k_mass]]
+trial-space   = "P1"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2, 3]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[k_mass]]
+trial-space   = "P1"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "AnnulusMap"
+
+[[k_mass]]
+trial-space   = "P1"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [3]
+quadrature    = 3
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "IcosahedralShellMap"
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [2, 3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [2]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "AnnulusMap"
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P1"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "IcosahedralShellMap"
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [2, 3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [2]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "AnnulusMap"
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P2"
+dimensions    = [3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "IcosahedralShellMap"
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2, 3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+operator-name = "KMassP1Coefficient"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [2]
+quadrature    = 4
+loop-strategy = "sawtooth"
+operator-name = "KMassP1Coefficient"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "AnnulusMap"
+
+[[k_mass]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.coefficient_function_space = "P1"
+dimensions    = [3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+operator-name = "KMassP1Coefficient"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "IcosahedralShellMap"

operators-mass.toml

+[[mass]]
+trial-space   = "P1"
+test-space    = "P1"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+
+[[mass]]
+trial-space   = "P1"
+test-space    = "P1"
+dimensions    = [2, 3]
+quadrature    = 2
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
+blending      = "ParametricP1Map"
+
+[[mass]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2, 3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[mass]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "AnnulusMap"
+
+[[mass]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "IcosahedralShellMap"
+
+[[mass]]
+trial-space   = "P2"
+test-space    = "P2"
+dimensions    = [2, 3]
+quadrature    = 4
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "ParametricP2Map"

operators-terraneo.toml

+[[shear_heating]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.viscosity_function_space = "P2"
+form-space-args.velocity_function_space = "P2"
+dimensions    = [2, 3]
+quadrature    = 3
+folder-name   = "terraneo"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[shear_heating]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.viscosity_function_space = "P2"
+form-space-args.velocity_function_space = "P2"
+dimensions    = [2]
+quadrature    = 3
+folder-name   = "terraneo"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "AnnulusMap"
+
+[[shear_heating]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.viscosity_function_space = "P2"
+form-space-args.velocity_function_space = "P2"
+dimensions    = [3]
+quadrature    = 3
+folder-name   = "terraneo"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "IcosahedralShellMap"
+
+[[shear_heating]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.viscosity_function_space = "P1"
+form-space-args.velocity_function_space = "P2"
+dimensions    = [2, 3]
+quadrature    = 3
+folder-name   = "terraneo"
+operator-name = "ShearHeatingP1Viscosity"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[shear_heating]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.viscosity_function_space = "P1"
+form-space-args.velocity_function_space = "P2"
+dimensions    = [2]
+quadrature    = 3
+folder-name   = "terraneo"
+operator-name = "ShearHeatingP1Viscosity"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "AnnulusMap"
+
+[[shear_heating]]
+trial-space   = "P2"
+test-space    = "P2"
+form-space-args.viscosity_function_space = "P1"
+form-space-args.velocity_function_space = "P2"
+dimensions    = [3]
+quadrature    = 3
+folder-name   = "terraneo"
+operator-name = "ShearHeatingP1Viscosity"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+blending      = "IcosahedralShellMap"
+
+[[grad_rho_by_rho_dot_u]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+form-space-args.density_function_space = "P2"
+dimensions    = [2, 3]
+quadrature    = 3
+operator-name = "FrozenVelocity"
+folder-name   = "terraneo"
+blending      = "IdentityMap"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[grad_rho_by_rho_dot_u]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+form-space-args.density_function_space = "P2"
+dimensions    = [2]
+quadrature    = 3
+operator-name = "FrozenVelocity"
+folder-name   = "terraneo"
+blending      = "AnnulusMap"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[grad_rho_by_rho_dot_u]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+form-space-args.density_function_space = "P2"
+dimensions    = [3]
+quadrature    = 3
+operator-name = "FrozenVelocity"
+folder-name   = "terraneo"
+blending      = "IcosahedralShellMap"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[grad_rho_by_rho_dot_u]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+form-space-args.density_function_space = "P1"
+dimensions    = [2, 3]
+quadrature    = 3
+operator-name = "FrozenVelocityP1Density"
+folder-name   = "terraneo"
+blending      = "IdentityMap"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[grad_rho_by_rho_dot_u]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+form-space-args.density_function_space = "P1"
+dimensions    = [2]
+quadrature    = 3
+operator-name = "FrozenVelocityP1Density"
+folder-name   = "terraneo"
+blending      = "AnnulusMap"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]
+
+[[grad_rho_by_rho_dot_u]]
+trial-space   = "P2Vector"
+test-space    = "P1"
+form-space-args.density_function_space = "P1"
+dimensions    = [3]
+quadrature    = 3
+operator-name = "FrozenVelocityP1Density"
+folder-name   = "terraneo"
+blending      = "IcosahedralShellMap"
+loop-strategy = "sawtooth"
+optimizations = ["moveconstants", "vectorize", "quadloops"]

operators.toml

-[[curl_curl]]
-trial-space   = "N1E1"
-test-space    = "N1E1"
-dimensions    = [3]
-quadrature    = 0
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-
-[[diffusion]]
-trial-space   = "P1"
-test-space    = "P1"
-dimensions    = [2, 3]
-quadrature    = 0
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-
-[[diffusion]]
-trial-space   = "P1"
-test-space    = "P1"
-dimensions    = [2, 3]
-quadrature    = 0
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "ParametricP1Map"
-
-[[diffusion]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-
-[[diffusion]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [2]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "AnnulusMap"
-
-[[diffusion]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [3]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "IcosahedralShellMap"
-
-[[diffusion]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "ParametricP2Map"
-
-[[div_k_grad]]
-trial-space   = "P1"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[div_k_grad]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[div_k_grad]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [2]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "AnnulusMap"
-
-[[div_k_grad]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [3]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "IcosahedralShellMap"
-
-[[div_k_grad]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-operator-name = "DivKGradP1Coefficient"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[div_k_grad]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2]
-quadrature    = 3
-loop-strategy = "sawtooth"
-operator-name = "DivKGradP1Coefficient"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "AnnulusMap"
-
-[[div_k_grad]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [3]
-quadrature    = 3
-loop-strategy = "sawtooth"
-operator-name = "DivKGradP1Coefficient"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "IcosahedralShellMap"
-
-[[mass]]
-trial-space   = "P1"
-test-space    = "P1"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[mass]]
-trial-space   = "P1"
-test-space    = "P1"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "ParametricP1Map"
-
-[[mass]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [2, 3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[mass]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [2]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "AnnulusMap"
-
-[[mass]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "IcosahedralShellMap"
-
-[[mass]]
-trial-space   = "P2"
-test-space    = "P2"
-dimensions    = [2, 3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "ParametricP2Map"
-
-[[k_mass]]
-trial-space   = "P1"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2, 3]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[k_mass]]
-trial-space   = "P1"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "AnnulusMap"
-
-[[k_mass]]
-trial-space   = "P1"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [3]
-quadrature    = 3
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "IcosahedralShellMap"
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [2, 3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [2]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "AnnulusMap"
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P1"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "IcosahedralShellMap"
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [2, 3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [2]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "AnnulusMap"
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P2"
-dimensions    = [3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "IcosahedralShellMap"
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2, 3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-operator-name = "KMassP1Coefficient"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [2]
-quadrature    = 4
-loop-strategy = "sawtooth"
-operator-name = "KMassP1Coefficient"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "AnnulusMap"
-
-[[k_mass]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.coefficient_function_space = "P1"
-dimensions    = [3]
-quadrature    = 4
-loop-strategy = "sawtooth"
-operator-name = "KMassP1Coefficient"
-optimizations = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending      = "IcosahedralShellMap"
-
-[[divergence]]
-trial-space   = "P2"
-test-space    = "P1"
-components-trial = [0, 1, 2]
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-
-[[divergence]]
-trial-space   = "P2"
-test-space    = "P1"
-components-trial = [0, 1]
-dimensions    = [2]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "AnnulusMap"
-
-[[divergence]]
-trial-space   = "P2"
-test-space    = "P1"
-components-trial = [0, 1, 2]
-dimensions    = [3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "IcosahedralShellMap"
-
-[[divergence]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-
-[[divergence]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-dimensions    = [2]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "AnnulusMap"
-
-[[divergence]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-dimensions    = [3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "IcosahedralShellMap"
-
-[[gradient]]
-trial-space   = "P1"
-test-space    = "P2"
-components-test = [0, 1, 2]
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-
-[[gradient]]
-trial-space   = "P1"
-test-space    = "P2"
-components-test = [0, 1]
-dimensions    = [2]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "AnnulusMap"
-
-[[gradient]]
-trial-space   = "P1"
-test-space    = "P2"
-components-test = [0, 1, 2]
-dimensions    = [3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "IcosahedralShellMap"
-
-[[gradient]]
-trial-space   = "P1"
-test-space    = "P2Vector"
-dimensions    = [2, 3]
-quadrature    = 2
-loop-strategy = "cubes"
-optimizations = ["moveconstants", "vectorize"]
-
-[[gradient]]
-trial-space   = "P1"
-test-space    = "P2Vector"
-dimensions    = [2]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "AnnulusMap"
-
-[[gradient]]
-trial-space   = "P1"
-test-space    = "P2Vector"
-dimensions    = [3]
-quadrature    = 2
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize"]
-blending      = "IcosahedralShellMap"
-
-[[epsilon]]
-trial-space               = "P2"
-test-space                = "P2"
-components-trial = [0, 1, 2]
-components-test  = [0, 1, 2]
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2, 3]
-quadrature      = 2
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[epsilon]]
-trial-space               = "P2"
-test-space                = "P2"
-components-trial = [0, 1]
-components-test  = [0, 1]
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "AnnulusMap"
-
-[[epsilon]]
-trial-space               = "P2"
-test-space                = "P2"
-components-trial = [0, 1, 2]
-components-test  = [0, 1, 2]
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [3]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "IcosahedralShellMap"
-
-[[epsilon]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2, 3]
-quadrature      = 2
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[epsilon]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "AnnulusMap"
-
-[[epsilon]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [3]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "IcosahedralShellMap"
-
-[[epsilon]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P1"
-dimensions      = [2]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-operator-name   = "EpsilonP1Viscosity"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "AnnulusMap"
-
-[[epsilon]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P1"
-dimensions      = [3]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-operator-name   = "EpsilonP1Viscosity"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "IcosahedralShellMap"
-
-[[full_stokes]]
-trial-space               = "P2"
-test-space                = "P2"
-components-trial = [0, 1, 2]
-components-test  = [0, 1, 2]
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2, 3]
-quadrature      = 2
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[full_stokes]]
-trial-space               = "P2"
-test-space                = "P2"
-components-trial = [0, 1]
-components-test  = [0, 1]
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "AnnulusMap"
-
-[[full_stokes]]
-trial-space               = "P2"
-test-space                = "P2"
-components-trial = [0, 1, 2]
-components-test  = [0, 1, 2]
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [3]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "IcosahedralShellMap"
-
-[[full_stokes]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2, 3]
-quadrature      = 2
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-
-[[full_stokes]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [2]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "AnnulusMap"
-
-[[full_stokes]]
-trial-space               = "P2Vector"
-test-space                = "P2Vector"
-form-space-args.coefficient_function_space = "P2"
-dimensions      = [3]
-quadrature      = 3
-loop-strategy   = "sawtooth"
-optimizations   = ["moveconstants", "vectorize", "quadloops", "tabulate"]
-blending        = "IcosahedralShellMap"
-
-[[shear_heating]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.viscosity_function_space = "P2"
-form-space-args.velocity_function_space = "P2"
-dimensions    = [2, 3]
-quadrature    = 3
-folder-name   = "terraneo"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[shear_heating]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.viscosity_function_space = "P2"
-form-space-args.velocity_function_space = "P2"
-dimensions    = [2]
-quadrature    = 3
-folder-name   = "terraneo"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "AnnulusMap"
-
-[[shear_heating]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.viscosity_function_space = "P2"
-form-space-args.velocity_function_space = "P2"
-dimensions    = [3]
-quadrature    = 3
-folder-name   = "terraneo"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "IcosahedralShellMap"
-
-[[shear_heating]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.viscosity_function_space = "P1"
-form-space-args.velocity_function_space = "P2"
-dimensions    = [2, 3]
-quadrature    = 3
-folder-name   = "terraneo"
-operator-name = "ShearHeatingP1Viscosity"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[shear_heating]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.viscosity_function_space = "P1"
-form-space-args.velocity_function_space = "P2"
-dimensions    = [2]
-quadrature    = 3
-folder-name   = "terraneo"
-operator-name = "ShearHeatingP1Viscosity"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "AnnulusMap"
-
-[[shear_heating]]
-trial-space   = "P2"
-test-space    = "P2"
-form-space-args.viscosity_function_space = "P1"
-form-space-args.velocity_function_space = "P2"
-dimensions    = [3]
-quadrature    = 3
-folder-name   = "terraneo"
-operator-name = "ShearHeatingP1Viscosity"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-blending      = "IcosahedralShellMap"
-
-[[grad_rho_by_rho_dot_u]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-form-space-args.density_function_space = "P2"
-dimensions    = [2, 3]
-quadrature    = 3
-operator-name = "FrozenVelocity"
-folder-name   = "terraneo"
-blending      = "IdentityMap"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[grad_rho_by_rho_dot_u]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-form-space-args.density_function_space = "P2"
-dimensions    = [2]
-quadrature    = 3
-operator-name = "FrozenVelocity"
-folder-name   = "terraneo"
-blending      = "AnnulusMap"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[grad_rho_by_rho_dot_u]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-form-space-args.density_function_space = "P2"
-dimensions    = [3]
-quadrature    = 3
-operator-name = "FrozenVelocity"
-folder-name   = "terraneo"
-blending      = "IcosahedralShellMap"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[grad_rho_by_rho_dot_u]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-form-space-args.density_function_space = "P1"
-dimensions    = [2, 3]
-quadrature    = 3
-operator-name = "FrozenVelocityP1Density"
-folder-name   = "terraneo"
-blending      = "IdentityMap"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[grad_rho_by_rho_dot_u]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-form-space-args.density_function_space = "P1"
-dimensions    = [2]
-quadrature    = 3
-operator-name = "FrozenVelocityP1Density"
-folder-name   = "terraneo"
-blending      = "AnnulusMap"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]
-
-[[grad_rho_by_rho_dot_u]]
-trial-space   = "P2Vector"
-test-space    = "P1"
-form-space-args.density_function_space = "P1"
-dimensions    = [3]
-quadrature    = 3
-operator-name = "FrozenVelocityP1Density"
-folder-name   = "terraneo"
-blending      = "IcosahedralShellMap"
-loop-strategy = "sawtooth"
-optimizations = ["moveconstants", "vectorize", "quadloops"]