fiqus 2025.2.0__py3-none-any.whl → 2025.10.0__py3-none-any.whl

fiqus/pro_material_functions/ironBHcurves.pro

@@ -1,246 +1,246 @@
-/* *****************************************************************************
-* CERN
-* BH curves for iron
-*******************************************************************************/
-
-Function{
-
-    magneticFluxBHiron1() = { // B(T)
-        0.000000, 0.800000, 1.200000, 1.400000, 1.500000, 1.550000, 1.600000, 1.650000, 1.700000, 1.750000, 1.800000,
-        1.850000, 1.900000, 2.000000, 2.050000, 2.100000, 2.150000, 2.200000, 2.249000, 2.275000, 2.300000, 2.344300,
-        2.399600, 2.490500, 2.562700, 2.670600, 2.849800, 3.207400, 4.278200, 4.813400, 5.705200, 6.418600, 7.488700
-	};
-    magneticFieldBHiron1() = { // H(A/m)
-        0.0, 159.2, 294.4, 501.3, 795.8, 1154.0, 1795.0, 2862.0, 4383.0, 6044.0, 8122.0, 10590.0, 13160.0, 21170.0,
-        26750.0, 33760.0, 43800.0, 66000.0, 99740.0, 120960.0, 141210.0, 169600.0, 212170.0, 283130.0, 339890.0,
-        425040.0, 566950.0, 850760.0, 1702300.0, 2128000.0, 2837700.0, 3405100.0, 4256700.0
-	};
-    iron_nuBHiron1()  = magneticFieldBHiron1() / magneticFluxBHiron1(); 	// nu = H/B
-    iron_nuBHiron1(0) = iron_nuBHiron1(1); // Avoid 0/0
-    BHiron1()   = ListAlt[magneticFluxBHiron1()^2, iron_nuBHiron1()];
-
-    magneticFluxBHiron2() = { // B(T)
-        0.000000, 0.800000, 1.200000, 1.400000, 1.500000, 1.550000, 1.600000, 1.650000, 1.700000, 1.750000, 1.800000,
-        1.850000, 1.900000, 2.000000, 2.050000, 2.100000, 2.150000, 2.200000, 2.249000, 2.275000, 2.300000, 2.344300,
-        2.399600, 2.490500, 2.562700, 2.670600, 2.849800, 3.207400, 4.278200, 4.813400, 5.705200, 6.418600, 7.488700, 10.0
-	};
-    magneticFieldBHiron2() = { // H(A/m)
-        0.0, 159.2, 294.4, 501.3, 795.8, 1154.0, 1795.0, 2862.0, 4383.0, 6044.0, 8122.0, 10590.0, 13160.0, 21170.0,
-        26750.0, 33760.0, 43800.0, 66000.0, 99740.0, 120960.0, 141210.0, 169600.0, 212170.0, 283130.0, 339890.0,
-        425040.0, 566950.0, 850760.0, 1702300.0, 2128000.0, 2837700.0, 3405100.0, 4256700.0, 6255129.043
-	};
-    /*magneticFluxBHiron2() = { // B(T)
-	    0.000, 0.012, 0.021, 0.038, 0.056, 0.097, 0.518, 0.683, 0.752, 0.899, 0.988, 1.082, 1.161, 1.237, 1.317, 1.374,
-	    1.432, 1.495, 1.505, 1.541, 1.572, 1.594, 1.614, 1.632, 1.650, 1.655, 1.666, 1.678, 1.682, 1.687, 1.698, 1.710,
-	    1.714, 1.722, 1.730, 1.748, 1.749, 1.766, 1.779, 1.784, 1.787, 1.805, 1.828, 1.850, 1.853, 1.875, 1.902, 1.931,
-	    1.963, 1.989, 1.994, 2.026, 2.080, 2.142, 2.209, 2.257, 2.285, 2.333, 2.386, 2.457, 2.497, 2.527, 2.557, 2.588,
-	    2.619, 2.650, 2.680, 4.278, 4.813, 5.705, 6.419, 7.489
-	};
-    magneticFieldBHiron2() = { // H(A/m)
-        0.00, 17.90, 24.60, 32.80, 39.60, 61.49, 122.97, 147.56, 159.86, 196.75, 221.35, 258.24, 295.13, 344.32, 405.80,
-        467.29, 541.07, 639.45, 651.74, 737.82, 860.79, 996.06, 1155.92, 1340.38, 1561.73, 1647.81, 1807.67, 2029.01,
-        2102.80, 2213.47, 2447.11, 2742.24, 2852.92, 3049.67, 3295.61, 3836.68, 3885.87, 4476.13, 4968.01, 5164.76,
-        5263.14, 6000.96, 7033.92, 8066.87, 8239.03, 9394.95, 10870.60, 12641.38, 14658.09, 16478.06, 17019.13,
-        19626.10, 26660.02, 35907.41, 48401.22, 65912.23, 82390.29, 118838.77, 159763.38, 215641.21, 247170.87,
-        270535.28, 295129.39, 319723.51, 344317.62, 368911.74, 393505.86, 1702300.00, 2128000.00, 2837700.00,
-        3405100.00, 4256700.00
-	}; */
-    iron_nuBHiron2()  = magneticFieldBHiron2() / magneticFluxBHiron2(); 	// nu = H/B
-    iron_nuBHiron2(0) = iron_nuBHiron2(1); // Avoid 0/0
-    BHiron2()   = ListAlt[magneticFluxBHiron2()^2, iron_nuBHiron2()];
-
-    magneticFluxBHiron3() = { // B(T)
-        0, 0.8, 1.2, 1.4, 1.5, 1.55, 1.6, 1.65, 1.7, 1.75, 1.8, 1.85, 1.9, 2, 2.05, 2.1, 2.15, 2.2, 2.25, 2.28, 2.308,
-        2.3447, 2.3996, 2.4905, 2.5627, 2.6706, 2.8498, 3.2074, 4.2782, 4.8134, 5.7055, 6.4187, 7.4888
-	};
-    magneticFieldBHiron3() = { // H(A/m)
-        0, 159.2, 294.4, 501.3, 795.8, 1154, 1795, 2862, 4383, 6044, 8122, 10590, 13160, 21170, 26750, 33760, 43800,
-        66000, 99740, 120960, 141210, 169600, 212170, 283130, 339890, 425040, 566950, 850760, 1702300, 2128000, 2837700,
-        3405100, 4256700
-	};
-	iron_nuBHiron3()  = magneticFieldBHiron3() / magneticFluxBHiron3(); 	// nu = H/B
-    iron_nuBHiron3(0) = iron_nuBHiron3(1); // Avoid 0/0
-    BHiron3()   = ListAlt[magneticFluxBHiron3()^2, iron_nuBHiron3()];
-
-	magneticFluxBHiron4() = { // B(T)
-        0.000000, 5.000000, 10.00000
-	};
-    magneticFieldBHiron4() = { // H(A/m)
-        0.000000, 3959078.186, 7918156.372
-	};
-	iron_nuBHiron4()  = magneticFieldBHiron4() / magneticFluxBHiron4(); 	// nu = H/B
-    iron_nuBHiron4(0) = iron_nuBHiron4(1); // Avoid 0/0
-    BHiron4()   = ListAlt[magneticFluxBHiron4()^2, iron_nuBHiron4()];
-
-	magneticFluxBHiron5() = { // B(T)
-        0.000000, 5.000000, 10.00000
-	};
-	magneticFieldBHiron5() = { // H(A/m)
-        0.000000, 3968951.199, 7937902.399
-	};
-	iron_nuBHiron5()  = magneticFieldBHiron5() / magneticFluxBHiron5(); 	// nu = H/B
-    iron_nuBHiron5(0) = iron_nuBHiron5(1); // Avoid 0/0
-    BHiron5()   = ListAlt[magneticFluxBHiron5()^2, iron_nuBHiron5()];
-
-	magneticFluxBHiron6() = { // B(T)
-        0.000000, 5.000000, 10.00000
-	};
-	magneticFieldBHiron6() = { // H(A/m)
-        0.000000, 3939478.789, 7878957.579
-	};
-	iron_nuBHiron6()  = magneticFieldBHiron6() / magneticFluxBHiron6(); 	// nu = H/B
-    iron_nuBHiron6(0) = iron_nuBHiron6(1); // Avoid 0/0
-    BHiron6()   = ListAlt[magneticFluxBHiron6()^2, iron_nuBHiron6()];
-
-    magneticFluxBHiron7() = { // B(T)
-        0.000000, 0.100000, 0.200000, 0.500000, 1.000000, 2.000000, 4.000000, 8.000000
-	};
-    magneticFieldBHiron7() = { // H(A/m)
-        0.000000, 76590.44422, 154744.7186, 391814.2371, 788285.9985, 1582057.088, 3168839.086, 6342097.752
-	};
-	iron_nuBHiron7()  = magneticFieldBHiron7() / magneticFluxBHiron7(); 	// nu = H/B
-    iron_nuBHiron7(0) = iron_nuBHiron7(1); // Avoid 0/0
-    BHiron7()   = ListAlt[magneticFluxBHiron7()^2, iron_nuBHiron7()];
-
-    magneticFluxBHiron8() = { // B(T)
-        0.000000, 0.900000, 1.800000
-	};
-    magneticFieldBHiron8() = { // H(A/m)
-        0.000000, 7.0E+05, 1.4E+06
-	};
-	iron_nuBHiron8()  = magneticFieldBHiron8() / magneticFluxBHiron8(); 	// nu = H/B
-    iron_nuBHiron8(0) = iron_nuBHiron8(1); // Avoid 0/0
-    BHiron8()   = ListAlt[magneticFluxBHiron8()^2, iron_nuBHiron8()];
-
-    magneticFluxBHiron9() = { // B(T)
-        0.000000, 0.900000, 1.800000
-	};
-    magneticFieldBHiron9() = { // H(A/m)
-        0.000000, 7.0E+05, 1.4E+06
-	};
-	iron_nuBHiron9()  = magneticFieldBHiron9() / magneticFluxBHiron9(); 	// nu = H/B
-    iron_nuBHiron9(0) = iron_nuBHiron9(1); // Avoid 0/0
-    BHiron9()   = ListAlt[magneticFluxBHiron9()^2, iron_nuBHiron9()];
-
-    magneticFluxBHiron10() = { // B(T)
-        0.0, 0.0066, 0.011, 0.0153, 0.0196, 0.0237, 0.0282, 0.0323, 0.0366, 0.0411, 0.0455, 0.0501, 0.0544, 0.0587,
-        0.0632, 0.0672, 0.0712, 0.0752, 0.0795, 0.0839, 0.088, 0.0921, 0.0964, 0.101, 0.1054, 0.1094, 0.1135, 0.1175,
-        0.1217, 0.1263, 0.131, 0.1356, 0.1401, 0.1441, 0.1484, 0.1524, 0.157, 0.1614, 0.1657, 0.1702, 0.1744, 0.1788,
-        0.183, 0.187, 0.1917, 0.1957, 0.2001, 0.2044, 0.2088, 0.2129, 0.217, 0.2214, 0.2259, 0.2304, 0.235, 0.2391,
-        0.2435, 0.2477, 0.2518, 0.2564, 0.2605, 0.2646, 0.2691, 0.2731, 0.2775, 0.2817, 0.2864, 0.2908, 0.2951, 0.2994,
-        0.3034, 0.3075, 0.3122, 0.3165, 0.3205, 0.3247, 0.329, 0.3337, 0.338, 0.3425, 0.3469, 0.3512, 0.3556, 0.3598,
-        0.3643, 0.3683, 0.3723, 0.3764, 0.3812, 0.3853, 0.3895, 0.3936, 0.3976, 0.4018, 0.4058, 0.4098, 0.4144, 0.4188,
-        0.4234, 0.4277, 0.432, 0.4361, 0.4402, 0.4446, 0.4488, 0.4533, 0.4574, 0.4616, 0.4657, 0.4698, 0.4739, 0.4783,
-        0.4829, 0.4871, 0.4911, 0.4954, 0.4998, 0.5038, 0.5083, 0.5124, 0.5168, 0.5212, 0.5252, 0.5295, 0.5329, 0.5357,
-        0.54, 0.5441, 0.5482, 0.5526, 0.5568, 0.5608, 0.5652, 0.5696, 0.5745, 0.5788, 0.5831, 0.5876, 0.592, 0.5963,
-        0.6007, 0.6048, 0.6088, 0.6133, 0.6175, 0.6221, 0.6267, 0.6308, 0.635, 0.6393, 0.6433, 0.6474, 0.6514, 0.6554,
-        0.6596, 0.664, 0.6683, 0.6727, 0.6771, 0.6818, 0.6861, 0.6905, 0.6945, 0.6988, 0.7029, 0.7076, 0.7119, 0.7163,
-        0.7204, 0.7248, 0.7289, 0.7333, 0.7375, 0.7419, 0.7462, 0.7503, 0.7527, 0.7546, 0.7593, 0.7633, 0.7678, 0.7718,
-        0.7761, 0.7805, 0.7849, 0.7892, 0.7936, 0.7977, 0.8021, 0.8063, 0.8109, 0.8156, 0.82, 0.8241, 0.8282, 0.8324,
-        0.8364, 0.8405, 0.8446, 0.8488, 0.8533, 0.8577, 0.8619, 0.8663, 0.8705, 0.875, 0.8792, 0.8832, 0.8874, 0.8915,
-        0.8957, 0.9001, 0.9042, 0.9084, 0.9125, 0.9165, 0.9209, 0.9249, 0.929, 0.9336, 0.9376, 0.9419, 0.9461, 0.9503,
-        0.9505, 0.9544, 0.9573, 0.9591, 0.963, 0.9675, 0.9719, 0.9759, 0.9799, 0.9844, 0.9884, 0.9926, 0.9969, 1.0009,
-        1.0049, 1.0092, 1.0136, 1.0177, 1.0221, 1.0268, 1.0312, 1.0353, 1.0394, 1.0434, 1.0479, 1.0524, 1.0569, 1.0613,
-        1.0657, 1.0703, 1.0744, 1.0789, 1.083, 1.0878, 1.0919, 1.0963, 1.1004, 1.1046, 1.1092, 1.1139, 1.1179, 1.1224,
-        1.1268, 1.131, 1.1353, 1.1393, 1.1438, 1.1478, 1.1518, 1.1562, 1.1606, 1.1648, 1.1692, 1.1738, 1.1782, 1.1829,
-        1.1871, 1.1914, 1.1954, 1.1997, 1.2038, 1.2084, 1.2124, 1.2167, 1.2211, 1.2252, 1.2295, 1.2336, 1.2382, 1.2426,
-        1.247, 1.2511, 1.2553, 1.2599, 1.2647, 1.269, 1.273, 1.2772, 1.2815, 1.2861, 1.2903, 1.2946, 1.299, 1.3037,
-        1.308, 1.3121, 1.3164, 1.3205, 1.3245, 1.3287, 1.3327, 1.3372, 1.342, 1.3464, 1.3506, 1.3549, 1.3592, 1.3635,
-        1.3677, 1.3717, 1.3761, 1.3803, 1.3845, 1.3889, 1.3929, 1.3975, 1.4017, 1.4062, 1.4105, 1.415, 1.4194, 1.4235,
-        1.4277, 1.4318, 1.4361, 1.4405, 1.4445, 1.4486, 1.4526, 1.4566, 1.461, 1.4651, 1.4694, 1.4735, 1.4777, 1.4823,
-        1.4866, 1.4907, 1.495, 1.4992, 1.5034, 1.508, 10.0
-	};
-    magneticFieldBHiron10() = { // H(A/m)
-        0.0, 10024.92038, 16280.65287, 21866.24204, 27652.86624, 32949.04459, 38722.92994, 43773.88535, 49027.86624,
-        54224.52229, 58894.90446, 63624.20382, 67796.97452, 71825.63694, 75840.76433, 79227.70701, 82551.75159,
-        85755.57325, 89135.35032, 92494.42675, 95514.33121, 98519.10828, 101640.12739, 104985.66879, 108128.18471,
-        110996.01911, 113856.6879, 116681.52866, 119693.47134, 123010.35032, 126435.50955, 129794.58599, 133046.17834,
-        135971.33758, 139121.01911, 142046.97452, 145521.49682, 148722.13376, 151929.93631, 155224.52229, 158380.57325,
-        161682.32484, 164853.50318, 167865.44586, 171419.58599, 174417.19745, 177804.93631, 180997.61146, 184342.35669,
-        187722.13376, 190468.15287, 193898.88535, 197351.91083, 200790.6051, 204265.92357, 207386.1465, 210745.22293,
-        213973.72611, 217166.40127, 220655.25478, 223804.93631, 226925.15924, 230392.51592, 233504.77707, 236921.97452,
-        240121.81529, 243777.07006, 247136.1465, 250458.59873, 253796.17834, 256865.44586, 260101.11465, 263699.04459,
-        266934.71338, 270105.09554, 273333.59873, 276699.84076, 280355.09554, 283620.22293, 287188.69427, 290576.43312,
-        293942.67516, 297344.74522, 300588.3758, 304098.72611, 307211.78344, 310338.3758, 313588.3758, 317280.25478,
-        320487.26115, 323614.64968, 326937.10191, 330121.81529, 333321.65605, 336448.24841, 339546.97452, 343151.27389,
-        346640.12739, 350230.89172, 353596.33758, 356948.24841, 360168.78981, 363375.79618, 366821.65605, 370064.49045,
-        373604.29936, 376839.17197, 380082.80255, 383332.80255, 386524.68153, 389760.35032, 393191.0828, 396839.17197,
-        400068.47134, 403238.05732, 406575.63694, 410078.82166, 413176.75159, 416745.22293, 419879.77707, 423354.29936,
-        426871.81529, 429984.87261, 433364.64968, 436014.33121, 438186.30573, 441566.0828, 444787.42038, 448051.75159,
-        451569.26752, 454769.10828, 457896.49682, 461341.56051, 464859.07643, 468717.35669, 472083.59873, 475492.03822,
-        478960.19108, 482463.3758, 485857.48408, 489317.67516, 492567.67516, 495622.61146, 499234.07643, 502520.70064,
-        506154.4586, 509802.54777, 512994.42675, 516281.05096, 519675.95541, 522845.5414, 526095.5414, 529165.6051,
-        532364.64968, 535621.81529, 539139.33121, 542534.23567, 546037.42038, 549483.28025, 553195.85987, 556562.10191,
-        560022.29299, 563177.54777, 566558.12102, 569793.78981, 573535.82803, 576858.28025, 580361.46497, 583625.79618,
-        587071.65605, 590335.98726, 593766.71975, 597110.66879, 600570.85987, 603951.43312, 607215.76433, 609097.92994,
-        610566.87898, 614294.58599, 617450.63694, 621033.43949, 624145.70064, 627569.26752, 631029.4586, 634518.3121,
-        637927.54777, 641387.73885, 644637.73885, 648111.46497, 651426.75159, 655060.50955, 658788.21656, 662262.73885,
-        665368.63057, 668704.61783, 672048.56688, 675204.61783, 678454.61783, 681675.15924, 684990.44586, 688544.58599,
-        692019.10828, 695319.26752, 698837.57962, 702123.40764, 705691.87898, 709035.82803, 712133.75796, 715535.82803,
-        718757.16561, 722042.99363, 725560.50955, 728767.51592, 732111.46497, 735361.46497, 738502.38854, 741991.24204,
-        745132.96178, 748397.29299, 752031.05096, 755201.43312, 758581.21019, 761910.82803, 765254.77707, 765385.35032,
-        768476.11465, 770698.24841, 772218.15287, 775323.24841, 778848.72611, 782337.57962, 785522.29299, 788677.54777,
-        792246.01911, 795402.07006, 798703.02548, 802126.59236, 805282.64331, 808452.2293, 811847.13376, 815335.98726,
-        818600.31847, 822089.17197, 825816.87898, 829277.07006, 832541.40127, 835820.06369, 838946.65605, 842500.79618,
-        846026.27389, 849609.07643, 853112.26115, 856571.65605, 860234.07643, 863469.74522, 867023.88535, 870288.21656,
-        874081.21019, 877345.5414, 880834.3949, 884012.02229, 887384.79299, 891061.6242, 894803.66242, 897945.14331,
-        901513.53503, 904987.89809, 908317.43631, 911712.18153, 914896.8949, 918450.8758, 921577.94586, 924777.07006,
-        928251.43312, 931783.43949, 935098.56688, 938601.75159, 942264.17197, 945666.48089, 949437.34076, 952766.87898,
-        956161.6242, 959317.51592, 962712.26115, 965991.00318, 969667.83439, 972823.72611, 976204.06051, 979707.24522,
-        982942.75478, 986380.73248, 989601.83121, 993235.42994, 996781.84713, 1000227.38854, 1003462.89809,
-        1006806.84713, 1010483.67834, 1014247.69108, 1017642.43631, 1020870.38217, 1024185.50955, 1027623.48726,
-        1031271.49682, 1034557.80255, 1038010.19108, 1041484.55414, 1045226.59236, 1048578.1051, 1051813.61465,
-        1055294.82484, 1058487.10191, 1061686.22611, 1065015.76433, 1068171.65605, 1071754.4586, 1075590.52548,
-        1079064.88854, 1082351.19427, 1085774.76115, 1089241.56051, 1092621.8949, 1095980.25478, 1099107.32484,
-        1102624.92038, 1105925.63694, 1109255.17516, 1112729.53822, 1115899.84076, 1119533.43949, 1122862.97771,
-        1126460.19108, 1129898.16879, 1133480.97134, 1136969.74522, 1140205.25478, 1143549.20382, 1146827.94586,
-        1150237.10191, 1153711.46497, 1156881.76752, 1160131.6879, 1163287.57962, 1166371.4172, 1169932.24522,
-        1173210.98726, 1176605.73248, 1179870.06369, 1183199.60191, 1186847.61146, 1190199.1242, 1193434.63376,
-        1196786.1465, 1200101.27389, 1203459.63376, 1207150.8758, 7880160.66488
-	};
-    iron_nuBHiron10()  = magneticFieldBHiron10() / magneticFluxBHiron10(); 	// nu = H/B
-    iron_nuBHiron10(0) = iron_nuBHiron10(1); // Avoid 0/0
-    BHiron10()   = ListAlt[magneticFluxBHiron10()^2, iron_nuBHiron10()];
-
-//
-    nuBHiron1[]        = InterpolationLinear[SquNorm[$1]]{ BHiron1() };
-    dnuHelperBHiron1[] = dInterpolationLinear[SquNorm[$1]]{ BHiron1() };
-    dnuBHiron1[] = 2*dnuHelperBHiron1[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron2[]        = InterpolationLinear[SquNorm[$1]]{ BHiron2() };
-    dnuHelperBHiron2[] = dInterpolationLinear[SquNorm[$1]]{ BHiron2() };
-    dnuBHiron2[] = 2*dnuHelperBHiron2[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron3[]        = InterpolationLinear[SquNorm[$1]]{ BHiron3() };
-    dnuHelperBHiron3[] = dInterpolationLinear[SquNorm[$1]]{ BHiron3() };
-    dnuBHiron3[] = 2*dnuHelperBHiron3[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron4[]        = InterpolationLinear[SquNorm[$1]]{ BHiron4() };
-    dnuHelperBHiron4[] = dInterpolationLinear[SquNorm[$1]]{ BHiron4() };
-    dnuBHiron4[] = 2*dnuHelperBHiron4[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron5[]        = InterpolationLinear[SquNorm[$1]]{ BHiron5() };
-    dnuHelperBHiron5[] = dInterpolationLinear[SquNorm[$1]]{ BHiron5() };
-    dnuBHiron5[] = 2*dnuHelperBHiron5[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron6[]        = InterpolationLinear[SquNorm[$1]]{ BHiron6() };
-    dnuHelperBHiron6[] = dInterpolationLinear[SquNorm[$1]]{ BHiron6() };
-    dnuBHiron6[] = 2*dnuHelperBHiron6[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron7[]        = InterpolationLinear[SquNorm[$1]]{ BHiron7() };
-    dnuHelperBHiron7[] = dInterpolationLinear[SquNorm[$1]]{ BHiron7() };
-    dnuBHiron7[] = 2*dnuHelperBHiron7[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron8[]        = InterpolationLinear[SquNorm[$1]]{ BHiron8() };
-    dnuHelperBHiron8[] = dInterpolationLinear[SquNorm[$1]]{ BHiron8() };
-    dnuBHiron8[] = 2*dnuHelperBHiron8[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron9[]        = InterpolationLinear[SquNorm[$1]]{ BHiron9() };
-    dnuHelperBHiron9[] = dInterpolationLinear[SquNorm[$1]]{ BHiron9() };
-    dnuBHiron9[] = 2*dnuHelperBHiron9[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-    nuBHiron10[]        = InterpolationLinear[SquNorm[$1]]{ BHiron10() };
-    dnuHelperBHiron10[] = dInterpolationLinear[SquNorm[$1]]{ BHiron10() };
-    dnuBHiron10[] = 2*dnuHelperBHiron10[$1] * SquDyadicProduct[$1] ; // Jacobian
-
-}
+/* *****************************************************************************
+* CERN
+* BH curves for iron
+*******************************************************************************/
+
+Function{
+
+    magneticFluxBHiron1() = { // B(T)
+        0.000000, 0.800000, 1.200000, 1.400000, 1.500000, 1.550000, 1.600000, 1.650000, 1.700000, 1.750000, 1.800000,
+        1.850000, 1.900000, 2.000000, 2.050000, 2.100000, 2.150000, 2.200000, 2.249000, 2.275000, 2.300000, 2.344300,
+        2.399600, 2.490500, 2.562700, 2.670600, 2.849800, 3.207400, 4.278200, 4.813400, 5.705200, 6.418600, 7.488700
+	};
+    magneticFieldBHiron1() = { // H(A/m)
+        0.0, 159.2, 294.4, 501.3, 795.8, 1154.0, 1795.0, 2862.0, 4383.0, 6044.0, 8122.0, 10590.0, 13160.0, 21170.0,
+        26750.0, 33760.0, 43800.0, 66000.0, 99740.0, 120960.0, 141210.0, 169600.0, 212170.0, 283130.0, 339890.0,
+        425040.0, 566950.0, 850760.0, 1702300.0, 2128000.0, 2837700.0, 3405100.0, 4256700.0
+	};
+    iron_nuBHiron1()  = magneticFieldBHiron1() / magneticFluxBHiron1(); 	// nu = H/B
+    iron_nuBHiron1(0) = iron_nuBHiron1(1); // Avoid 0/0
+    BHiron1()   = ListAlt[magneticFluxBHiron1()^2, iron_nuBHiron1()];
+
+    magneticFluxBHiron2() = { // B(T)
+        0.000000, 0.800000, 1.200000, 1.400000, 1.500000, 1.550000, 1.600000, 1.650000, 1.700000, 1.750000, 1.800000,
+        1.850000, 1.900000, 2.000000, 2.050000, 2.100000, 2.150000, 2.200000, 2.249000, 2.275000, 2.300000, 2.344300,
+        2.399600, 2.490500, 2.562700, 2.670600, 2.849800, 3.207400, 4.278200, 4.813400, 5.705200, 6.418600, 7.488700, 10.0
+	};
+    magneticFieldBHiron2() = { // H(A/m)
+        0.0, 159.2, 294.4, 501.3, 795.8, 1154.0, 1795.0, 2862.0, 4383.0, 6044.0, 8122.0, 10590.0, 13160.0, 21170.0,
+        26750.0, 33760.0, 43800.0, 66000.0, 99740.0, 120960.0, 141210.0, 169600.0, 212170.0, 283130.0, 339890.0,
+        425040.0, 566950.0, 850760.0, 1702300.0, 2128000.0, 2837700.0, 3405100.0, 4256700.0, 6255129.043
+	};
+    /*magneticFluxBHiron2() = { // B(T)
+	    0.000, 0.012, 0.021, 0.038, 0.056, 0.097, 0.518, 0.683, 0.752, 0.899, 0.988, 1.082, 1.161, 1.237, 1.317, 1.374,
+	    1.432, 1.495, 1.505, 1.541, 1.572, 1.594, 1.614, 1.632, 1.650, 1.655, 1.666, 1.678, 1.682, 1.687, 1.698, 1.710,
+	    1.714, 1.722, 1.730, 1.748, 1.749, 1.766, 1.779, 1.784, 1.787, 1.805, 1.828, 1.850, 1.853, 1.875, 1.902, 1.931,
+	    1.963, 1.989, 1.994, 2.026, 2.080, 2.142, 2.209, 2.257, 2.285, 2.333, 2.386, 2.457, 2.497, 2.527, 2.557, 2.588,
+	    2.619, 2.650, 2.680, 4.278, 4.813, 5.705, 6.419, 7.489
+	};
+    magneticFieldBHiron2() = { // H(A/m)
+        0.00, 17.90, 24.60, 32.80, 39.60, 61.49, 122.97, 147.56, 159.86, 196.75, 221.35, 258.24, 295.13, 344.32, 405.80,
+        467.29, 541.07, 639.45, 651.74, 737.82, 860.79, 996.06, 1155.92, 1340.38, 1561.73, 1647.81, 1807.67, 2029.01,
+        2102.80, 2213.47, 2447.11, 2742.24, 2852.92, 3049.67, 3295.61, 3836.68, 3885.87, 4476.13, 4968.01, 5164.76,
+        5263.14, 6000.96, 7033.92, 8066.87, 8239.03, 9394.95, 10870.60, 12641.38, 14658.09, 16478.06, 17019.13,
+        19626.10, 26660.02, 35907.41, 48401.22, 65912.23, 82390.29, 118838.77, 159763.38, 215641.21, 247170.87,
+        270535.28, 295129.39, 319723.51, 344317.62, 368911.74, 393505.86, 1702300.00, 2128000.00, 2837700.00,
+        3405100.00, 4256700.00
+	}; */
+    iron_nuBHiron2()  = magneticFieldBHiron2() / magneticFluxBHiron2(); 	// nu = H/B
+    iron_nuBHiron2(0) = iron_nuBHiron2(1); // Avoid 0/0
+    BHiron2()   = ListAlt[magneticFluxBHiron2()^2, iron_nuBHiron2()];
+
+    magneticFluxBHiron3() = { // B(T)
+        0, 0.8, 1.2, 1.4, 1.5, 1.55, 1.6, 1.65, 1.7, 1.75, 1.8, 1.85, 1.9, 2, 2.05, 2.1, 2.15, 2.2, 2.25, 2.28, 2.308,
+        2.3447, 2.3996, 2.4905, 2.5627, 2.6706, 2.8498, 3.2074, 4.2782, 4.8134, 5.7055, 6.4187, 7.4888
+	};
+    magneticFieldBHiron3() = { // H(A/m)
+        0, 159.2, 294.4, 501.3, 795.8, 1154, 1795, 2862, 4383, 6044, 8122, 10590, 13160, 21170, 26750, 33760, 43800,
+        66000, 99740, 120960, 141210, 169600, 212170, 283130, 339890, 425040, 566950, 850760, 1702300, 2128000, 2837700,
+        3405100, 4256700
+	};
+	iron_nuBHiron3()  = magneticFieldBHiron3() / magneticFluxBHiron3(); 	// nu = H/B
+    iron_nuBHiron3(0) = iron_nuBHiron3(1); // Avoid 0/0
+    BHiron3()   = ListAlt[magneticFluxBHiron3()^2, iron_nuBHiron3()];
+
+	magneticFluxBHiron4() = { // B(T)
+        0.000000, 5.000000, 10.00000
+	};
+    magneticFieldBHiron4() = { // H(A/m)
+        0.000000, 3959078.186, 7918156.372
+	};
+	iron_nuBHiron4()  = magneticFieldBHiron4() / magneticFluxBHiron4(); 	// nu = H/B
+    iron_nuBHiron4(0) = iron_nuBHiron4(1); // Avoid 0/0
+    BHiron4()   = ListAlt[magneticFluxBHiron4()^2, iron_nuBHiron4()];
+
+	magneticFluxBHiron5() = { // B(T)
+        0.000000, 5.000000, 10.00000
+	};
+	magneticFieldBHiron5() = { // H(A/m)
+        0.000000, 3968951.199, 7937902.399
+	};
+	iron_nuBHiron5()  = magneticFieldBHiron5() / magneticFluxBHiron5(); 	// nu = H/B
+    iron_nuBHiron5(0) = iron_nuBHiron5(1); // Avoid 0/0
+    BHiron5()   = ListAlt[magneticFluxBHiron5()^2, iron_nuBHiron5()];
+
+	magneticFluxBHiron6() = { // B(T)
+        0.000000, 5.000000, 10.00000
+	};
+	magneticFieldBHiron6() = { // H(A/m)
+        0.000000, 3939478.789, 7878957.579
+	};
+	iron_nuBHiron6()  = magneticFieldBHiron6() / magneticFluxBHiron6(); 	// nu = H/B
+    iron_nuBHiron6(0) = iron_nuBHiron6(1); // Avoid 0/0
+    BHiron6()   = ListAlt[magneticFluxBHiron6()^2, iron_nuBHiron6()];
+
+    magneticFluxBHiron7() = { // B(T)
+        0.000000, 0.100000, 0.200000, 0.500000, 1.000000, 2.000000, 4.000000, 8.000000
+	};
+    magneticFieldBHiron7() = { // H(A/m)
+        0.000000, 76590.44422, 154744.7186, 391814.2371, 788285.9985, 1582057.088, 3168839.086, 6342097.752
+	};
+	iron_nuBHiron7()  = magneticFieldBHiron7() / magneticFluxBHiron7(); 	// nu = H/B
+    iron_nuBHiron7(0) = iron_nuBHiron7(1); // Avoid 0/0
+    BHiron7()   = ListAlt[magneticFluxBHiron7()^2, iron_nuBHiron7()];
+
+    magneticFluxBHiron8() = { // B(T)
+        0.000000, 0.900000, 1.800000
+	};
+    magneticFieldBHiron8() = { // H(A/m)
+        0.000000, 7.0E+05, 1.4E+06
+	};
+	iron_nuBHiron8()  = magneticFieldBHiron8() / magneticFluxBHiron8(); 	// nu = H/B
+    iron_nuBHiron8(0) = iron_nuBHiron8(1); // Avoid 0/0
+    BHiron8()   = ListAlt[magneticFluxBHiron8()^2, iron_nuBHiron8()];
+
+    magneticFluxBHiron9() = { // B(T)
+        0.000000, 0.900000, 1.800000
+	};
+    magneticFieldBHiron9() = { // H(A/m)
+        0.000000, 7.0E+05, 1.4E+06
+	};
+	iron_nuBHiron9()  = magneticFieldBHiron9() / magneticFluxBHiron9(); 	// nu = H/B
+    iron_nuBHiron9(0) = iron_nuBHiron9(1); // Avoid 0/0
+    BHiron9()   = ListAlt[magneticFluxBHiron9()^2, iron_nuBHiron9()];
+
+    magneticFluxBHiron10() = { // B(T)
+        0.0, 0.0066, 0.011, 0.0153, 0.0196, 0.0237, 0.0282, 0.0323, 0.0366, 0.0411, 0.0455, 0.0501, 0.0544, 0.0587,
+        0.0632, 0.0672, 0.0712, 0.0752, 0.0795, 0.0839, 0.088, 0.0921, 0.0964, 0.101, 0.1054, 0.1094, 0.1135, 0.1175,
+        0.1217, 0.1263, 0.131, 0.1356, 0.1401, 0.1441, 0.1484, 0.1524, 0.157, 0.1614, 0.1657, 0.1702, 0.1744, 0.1788,
+        0.183, 0.187, 0.1917, 0.1957, 0.2001, 0.2044, 0.2088, 0.2129, 0.217, 0.2214, 0.2259, 0.2304, 0.235, 0.2391,
+        0.2435, 0.2477, 0.2518, 0.2564, 0.2605, 0.2646, 0.2691, 0.2731, 0.2775, 0.2817, 0.2864, 0.2908, 0.2951, 0.2994,
+        0.3034, 0.3075, 0.3122, 0.3165, 0.3205, 0.3247, 0.329, 0.3337, 0.338, 0.3425, 0.3469, 0.3512, 0.3556, 0.3598,
+        0.3643, 0.3683, 0.3723, 0.3764, 0.3812, 0.3853, 0.3895, 0.3936, 0.3976, 0.4018, 0.4058, 0.4098, 0.4144, 0.4188,
+        0.4234, 0.4277, 0.432, 0.4361, 0.4402, 0.4446, 0.4488, 0.4533, 0.4574, 0.4616, 0.4657, 0.4698, 0.4739, 0.4783,
+        0.4829, 0.4871, 0.4911, 0.4954, 0.4998, 0.5038, 0.5083, 0.5124, 0.5168, 0.5212, 0.5252, 0.5295, 0.5329, 0.5357,
+        0.54, 0.5441, 0.5482, 0.5526, 0.5568, 0.5608, 0.5652, 0.5696, 0.5745, 0.5788, 0.5831, 0.5876, 0.592, 0.5963,
+        0.6007, 0.6048, 0.6088, 0.6133, 0.6175, 0.6221, 0.6267, 0.6308, 0.635, 0.6393, 0.6433, 0.6474, 0.6514, 0.6554,
+        0.6596, 0.664, 0.6683, 0.6727, 0.6771, 0.6818, 0.6861, 0.6905, 0.6945, 0.6988, 0.7029, 0.7076, 0.7119, 0.7163,
+        0.7204, 0.7248, 0.7289, 0.7333, 0.7375, 0.7419, 0.7462, 0.7503, 0.7527, 0.7546, 0.7593, 0.7633, 0.7678, 0.7718,
+        0.7761, 0.7805, 0.7849, 0.7892, 0.7936, 0.7977, 0.8021, 0.8063, 0.8109, 0.8156, 0.82, 0.8241, 0.8282, 0.8324,
+        0.8364, 0.8405, 0.8446, 0.8488, 0.8533, 0.8577, 0.8619, 0.8663, 0.8705, 0.875, 0.8792, 0.8832, 0.8874, 0.8915,
+        0.8957, 0.9001, 0.9042, 0.9084, 0.9125, 0.9165, 0.9209, 0.9249, 0.929, 0.9336, 0.9376, 0.9419, 0.9461, 0.9503,
+        0.9505, 0.9544, 0.9573, 0.9591, 0.963, 0.9675, 0.9719, 0.9759, 0.9799, 0.9844, 0.9884, 0.9926, 0.9969, 1.0009,
+        1.0049, 1.0092, 1.0136, 1.0177, 1.0221, 1.0268, 1.0312, 1.0353, 1.0394, 1.0434, 1.0479, 1.0524, 1.0569, 1.0613,
+        1.0657, 1.0703, 1.0744, 1.0789, 1.083, 1.0878, 1.0919, 1.0963, 1.1004, 1.1046, 1.1092, 1.1139, 1.1179, 1.1224,
+        1.1268, 1.131, 1.1353, 1.1393, 1.1438, 1.1478, 1.1518, 1.1562, 1.1606, 1.1648, 1.1692, 1.1738, 1.1782, 1.1829,
+        1.1871, 1.1914, 1.1954, 1.1997, 1.2038, 1.2084, 1.2124, 1.2167, 1.2211, 1.2252, 1.2295, 1.2336, 1.2382, 1.2426,
+        1.247, 1.2511, 1.2553, 1.2599, 1.2647, 1.269, 1.273, 1.2772, 1.2815, 1.2861, 1.2903, 1.2946, 1.299, 1.3037,
+        1.308, 1.3121, 1.3164, 1.3205, 1.3245, 1.3287, 1.3327, 1.3372, 1.342, 1.3464, 1.3506, 1.3549, 1.3592, 1.3635,
+        1.3677, 1.3717, 1.3761, 1.3803, 1.3845, 1.3889, 1.3929, 1.3975, 1.4017, 1.4062, 1.4105, 1.415, 1.4194, 1.4235,
+        1.4277, 1.4318, 1.4361, 1.4405, 1.4445, 1.4486, 1.4526, 1.4566, 1.461, 1.4651, 1.4694, 1.4735, 1.4777, 1.4823,
+        1.4866, 1.4907, 1.495, 1.4992, 1.5034, 1.508, 10.0
+	};
+    magneticFieldBHiron10() = { // H(A/m)
+        0.0, 10024.92038, 16280.65287, 21866.24204, 27652.86624, 32949.04459, 38722.92994, 43773.88535, 49027.86624,
+        54224.52229, 58894.90446, 63624.20382, 67796.97452, 71825.63694, 75840.76433, 79227.70701, 82551.75159,
+        85755.57325, 89135.35032, 92494.42675, 95514.33121, 98519.10828, 101640.12739, 104985.66879, 108128.18471,
+        110996.01911, 113856.6879, 116681.52866, 119693.47134, 123010.35032, 126435.50955, 129794.58599, 133046.17834,
+        135971.33758, 139121.01911, 142046.97452, 145521.49682, 148722.13376, 151929.93631, 155224.52229, 158380.57325,
+        161682.32484, 164853.50318, 167865.44586, 171419.58599, 174417.19745, 177804.93631, 180997.61146, 184342.35669,
+        187722.13376, 190468.15287, 193898.88535, 197351.91083, 200790.6051, 204265.92357, 207386.1465, 210745.22293,
+        213973.72611, 217166.40127, 220655.25478, 223804.93631, 226925.15924, 230392.51592, 233504.77707, 236921.97452,
+        240121.81529, 243777.07006, 247136.1465, 250458.59873, 253796.17834, 256865.44586, 260101.11465, 263699.04459,
+        266934.71338, 270105.09554, 273333.59873, 276699.84076, 280355.09554, 283620.22293, 287188.69427, 290576.43312,
+        293942.67516, 297344.74522, 300588.3758, 304098.72611, 307211.78344, 310338.3758, 313588.3758, 317280.25478,
+        320487.26115, 323614.64968, 326937.10191, 330121.81529, 333321.65605, 336448.24841, 339546.97452, 343151.27389,
+        346640.12739, 350230.89172, 353596.33758, 356948.24841, 360168.78981, 363375.79618, 366821.65605, 370064.49045,
+        373604.29936, 376839.17197, 380082.80255, 383332.80255, 386524.68153, 389760.35032, 393191.0828, 396839.17197,
+        400068.47134, 403238.05732, 406575.63694, 410078.82166, 413176.75159, 416745.22293, 419879.77707, 423354.29936,
+        426871.81529, 429984.87261, 433364.64968, 436014.33121, 438186.30573, 441566.0828, 444787.42038, 448051.75159,
+        451569.26752, 454769.10828, 457896.49682, 461341.56051, 464859.07643, 468717.35669, 472083.59873, 475492.03822,
+        478960.19108, 482463.3758, 485857.48408, 489317.67516, 492567.67516, 495622.61146, 499234.07643, 502520.70064,
+        506154.4586, 509802.54777, 512994.42675, 516281.05096, 519675.95541, 522845.5414, 526095.5414, 529165.6051,
+        532364.64968, 535621.81529, 539139.33121, 542534.23567, 546037.42038, 549483.28025, 553195.85987, 556562.10191,
+        560022.29299, 563177.54777, 566558.12102, 569793.78981, 573535.82803, 576858.28025, 580361.46497, 583625.79618,
+        587071.65605, 590335.98726, 593766.71975, 597110.66879, 600570.85987, 603951.43312, 607215.76433, 609097.92994,
+        610566.87898, 614294.58599, 617450.63694, 621033.43949, 624145.70064, 627569.26752, 631029.4586, 634518.3121,
+        637927.54777, 641387.73885, 644637.73885, 648111.46497, 651426.75159, 655060.50955, 658788.21656, 662262.73885,
+        665368.63057, 668704.61783, 672048.56688, 675204.61783, 678454.61783, 681675.15924, 684990.44586, 688544.58599,
+        692019.10828, 695319.26752, 698837.57962, 702123.40764, 705691.87898, 709035.82803, 712133.75796, 715535.82803,
+        718757.16561, 722042.99363, 725560.50955, 728767.51592, 732111.46497, 735361.46497, 738502.38854, 741991.24204,
+        745132.96178, 748397.29299, 752031.05096, 755201.43312, 758581.21019, 761910.82803, 765254.77707, 765385.35032,
+        768476.11465, 770698.24841, 772218.15287, 775323.24841, 778848.72611, 782337.57962, 785522.29299, 788677.54777,
+        792246.01911, 795402.07006, 798703.02548, 802126.59236, 805282.64331, 808452.2293, 811847.13376, 815335.98726,
+        818600.31847, 822089.17197, 825816.87898, 829277.07006, 832541.40127, 835820.06369, 838946.65605, 842500.79618,
+        846026.27389, 849609.07643, 853112.26115, 856571.65605, 860234.07643, 863469.74522, 867023.88535, 870288.21656,
+        874081.21019, 877345.5414, 880834.3949, 884012.02229, 887384.79299, 891061.6242, 894803.66242, 897945.14331,
+        901513.53503, 904987.89809, 908317.43631, 911712.18153, 914896.8949, 918450.8758, 921577.94586, 924777.07006,
+        928251.43312, 931783.43949, 935098.56688, 938601.75159, 942264.17197, 945666.48089, 949437.34076, 952766.87898,
+        956161.6242, 959317.51592, 962712.26115, 965991.00318, 969667.83439, 972823.72611, 976204.06051, 979707.24522,
+        982942.75478, 986380.73248, 989601.83121, 993235.42994, 996781.84713, 1000227.38854, 1003462.89809,
+        1006806.84713, 1010483.67834, 1014247.69108, 1017642.43631, 1020870.38217, 1024185.50955, 1027623.48726,
+        1031271.49682, 1034557.80255, 1038010.19108, 1041484.55414, 1045226.59236, 1048578.1051, 1051813.61465,
+        1055294.82484, 1058487.10191, 1061686.22611, 1065015.76433, 1068171.65605, 1071754.4586, 1075590.52548,
+        1079064.88854, 1082351.19427, 1085774.76115, 1089241.56051, 1092621.8949, 1095980.25478, 1099107.32484,
+        1102624.92038, 1105925.63694, 1109255.17516, 1112729.53822, 1115899.84076, 1119533.43949, 1122862.97771,
+        1126460.19108, 1129898.16879, 1133480.97134, 1136969.74522, 1140205.25478, 1143549.20382, 1146827.94586,
+        1150237.10191, 1153711.46497, 1156881.76752, 1160131.6879, 1163287.57962, 1166371.4172, 1169932.24522,
+        1173210.98726, 1176605.73248, 1179870.06369, 1183199.60191, 1186847.61146, 1190199.1242, 1193434.63376,
+        1196786.1465, 1200101.27389, 1203459.63376, 1207150.8758, 7880160.66488
+	};
+    iron_nuBHiron10()  = magneticFieldBHiron10() / magneticFluxBHiron10(); 	// nu = H/B
+    iron_nuBHiron10(0) = iron_nuBHiron10(1); // Avoid 0/0
+    BHiron10()   = ListAlt[magneticFluxBHiron10()^2, iron_nuBHiron10()];
+
+//
+    nuBHiron1[]        = InterpolationLinear[SquNorm[$1]]{ BHiron1() };
+    dnuHelperBHiron1[] = dInterpolationLinear[SquNorm[$1]]{ BHiron1() };
+    dnuBHiron1[] = 2*dnuHelperBHiron1[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron2[]        = InterpolationLinear[SquNorm[$1]]{ BHiron2() };
+    dnuHelperBHiron2[] = dInterpolationLinear[SquNorm[$1]]{ BHiron2() };
+    dnuBHiron2[] = 2*dnuHelperBHiron2[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron3[]        = InterpolationLinear[SquNorm[$1]]{ BHiron3() };
+    dnuHelperBHiron3[] = dInterpolationLinear[SquNorm[$1]]{ BHiron3() };
+    dnuBHiron3[] = 2*dnuHelperBHiron3[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron4[]        = InterpolationLinear[SquNorm[$1]]{ BHiron4() };
+    dnuHelperBHiron4[] = dInterpolationLinear[SquNorm[$1]]{ BHiron4() };
+    dnuBHiron4[] = 2*dnuHelperBHiron4[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron5[]        = InterpolationLinear[SquNorm[$1]]{ BHiron5() };
+    dnuHelperBHiron5[] = dInterpolationLinear[SquNorm[$1]]{ BHiron5() };
+    dnuBHiron5[] = 2*dnuHelperBHiron5[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron6[]        = InterpolationLinear[SquNorm[$1]]{ BHiron6() };
+    dnuHelperBHiron6[] = dInterpolationLinear[SquNorm[$1]]{ BHiron6() };
+    dnuBHiron6[] = 2*dnuHelperBHiron6[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron7[]        = InterpolationLinear[SquNorm[$1]]{ BHiron7() };
+    dnuHelperBHiron7[] = dInterpolationLinear[SquNorm[$1]]{ BHiron7() };
+    dnuBHiron7[] = 2*dnuHelperBHiron7[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron8[]        = InterpolationLinear[SquNorm[$1]]{ BHiron8() };
+    dnuHelperBHiron8[] = dInterpolationLinear[SquNorm[$1]]{ BHiron8() };
+    dnuBHiron8[] = 2*dnuHelperBHiron8[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron9[]        = InterpolationLinear[SquNorm[$1]]{ BHiron9() };
+    dnuHelperBHiron9[] = dInterpolationLinear[SquNorm[$1]]{ BHiron9() };
+    dnuBHiron9[] = 2*dnuHelperBHiron9[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+    nuBHiron10[]        = InterpolationLinear[SquNorm[$1]]{ BHiron10() };
+    dnuHelperBHiron10[] = dInterpolationLinear[SquNorm[$1]]{ BHiron10() };
+    dnuBHiron10[] = 2*dnuHelperBHiron10[$1] * SquDyadicProduct[$1] ; // Jacobian
+
+}