added :mag to default_dec_generate_matrix

docs/src/cism/cism.md

@@ -202,12 +202,8 @@ We provide here the mapping from symbols to differential operators. As more of t
 function generate(sd, my_symbol; hodge=GeometricHodge())
   # We pre-allocate matrices that encode differential operators.
   op = @match my_symbol begin
-    :mag => x -> norm.(x)
-    :♯  => begin
-      sharp_mat = ♯_mat(sd, AltPPSharp())
-      x -> sharp_mat * x
-    end
-    x => error("Unmatched operator $my_symbol")
+    x => default_dec_matrix_generate(sd, x, hodge) 
+    _ => error("Unmatched operator $my_symbol")
   end
   return (args...) -> op(args...)
 end

docs/src/ebm_melt/ebm_melt.md

@@ -273,14 +273,8 @@ constants_and_parameters = (
 # Define how symbols map to Julia functions
 function generate(sd, my_symbol; hodge=GeometricHodge())
   op = @match my_symbol begin
-    :♯ => begin
-      sharp_mat = ♯_mat(sd, AltPPSharp())
-      x -> sharp_mat * x
-    end
-    :mag => x -> begin
-      norm.(x)
-    end
-    x => error("Unmatched operator $my_symbol")
+    x => default_dec_matrix_generate(sd, x, hodge)
+    _ => error("Unmatched operator $my_symbol")
   end
   return (args...) -> op(args...)
 end

docs/src/grigoriev/grigoriev.md

@@ -150,12 +150,8 @@ to_graphviz(ice_dynamics)
 ``` @example DEC
 function generate(sd, my_symbol; hodge=GeometricHodge())
   op = @match my_symbol begin
-    :mag => x -> norm.(x)
-    :♯ => begin
-      sharp_mat = ♯_mat(sd, AltPPSharp())
-      x -> sharp_mat * x
-    end
-    x => error("Unmatched operator $my_symbol")
+    x => default_dec_matrix_generate(sd, x, hodge)
+    _ => error("Unmatched operator $my_symbol")
   end
   return op
 end

src/operators.jl

@@ -80,6 +80,7 @@ function default_dec_matrix_generate(sd::HasDeltaSet, my_symbol::Symbol, hodge::
     :avg₀₁ => dec_avg₀₁(sd)
 
     :neg => x -> -1 .* x
+    :mag => x -> norm.(x)
      _ => error("Unmatched operator $my_symbol")
   end