{ "schema_version": 1, "reference_code": { "project": "NeoPZ", "repository": "https://github.com/labmec/neopz", "revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "modules": [ "Material/DarcyFlow/TPZMixedDarcyFlow.cpp", "Mesh/pzelchdiv.cpp", "Pre/TPZMHMixedMeshControl.cpp", "Pre/TPZMHMeshControl.cpp" ], "linear_solvers": { "default": "native NeoPZ skyline LDLt", "tetrahedral_well_mhm": "native assembled free matrix solved by SciPy SuperLU with partial pivoting and explicit symmetric equilibration; physical coefficients reimported and original matrix residual and fields evaluated in NeoPZ" }, "native_library_modified": false, "adapter_scope": "Original input/field bridge, boundary quadrature policy and exact cap EquationFilter; all shape functions, element forms, restrictions, assembly and physical field evaluation are NeoPZ; each recorded linear solver is identified separately." }, "data": { "affine_patch": "p=1+x+2y+3z, K=I, f=0, q=(-1,-2,-3), exact pressure on every boundary face.", "well": "Octagonal annular prism ri=0.2m, re=50m, height10m, four geometric radial intervals; K/eta=1e-10 SI, total influx0.01m3/s, p(re)=25MPa; exact radial pressure evaluated on actual planar inner/outer walls, zero normal flux on horizontal caps.", "pressure": "Absolute pressure and increment p-25MPa compared separately for well cases.", "integration": "Both fields evaluated at identical physical interior quadrature points on geometrically identical fine cells; no nodal coefficient comparison or smoothing.", "trace_spaces": "P1 triangular normal traces and Q1 rectangular traces; tetrahedral volume pressure P1 or P2, prism pressure P1(triangle) tensor P1(interval)." }, "classical_cases": [ { "kind": "prism", "degree": 1, "problem": "affine", "fine_factor": 1, "cells": 2, "pressure_l2_difference": 7.302524790888654e-15, "flux_l2_difference": 7.677337242460786e-14, "pressure_relative_difference": 1.762503762797775e-15, "flux_relative_difference": 2.0518546860005778e-14, "native_pressure_l2_error": 1.0314782561325539e-15, "native_flux_l2_error": 2.1072520114612865e-14, "pymhm_pressure_l2_error": 6.753163969314535e-15, "pymhm_flux_l2_error": 7.382368766114369e-14, "native_divergence_l2": 6.462597372198978e-14, "pymhm_divergence_l2": 6.224409805196977e-14, "pymhm_physical_residual": 2.8694650225310296e-15, "native_metadata": { "equations": 62, "active_equations": 62, "relative_residual": 3.8093999874325557e-16 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "driver_sha256": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "binary_sha256": "f051b51711f75646468f527abfa7d452b594f4e7f4fc13e07c520b1cc215818c", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "source_changed": false, "pressure_increment_relative_difference": 1.762503762797775e-15, "fields_sha256": "fb264eb140d80be5a795c566c3b63c22e3d79f796bcc3e5f1e4c06a34b5b4357", "native_pressure_l2_norm": 4.143267631552017, "native_pressure_increment_l2_norm": 4.143267631552017, "native_flux_l2_norm": 3.7416573867739404, "pressure_increment_datum": 0.0, "comparison_method": "classical conforming NeoPZ mixed finite elements; identical complete fine trace in PyMHM", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 27, "pressure_space_dimension_per_fine_cell": 6, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "run_compare.py": "7f5a2fd981dc73549b6dcc42c5ad5fc48e4cf548318d90928e4777d280b2dc45" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "degree": 1, "problem": "affine", "fine_factor": 1, "cells": 6, "pressure_l2_difference": 8.536788064611237e-15, "flux_l2_difference": 6.213770060952358e-14, "pressure_relative_difference": 2.0603998640110695e-15, "flux_relative_difference": 1.6606999034483684e-14, "native_pressure_l2_error": 2.7783485719674606e-15, "native_flux_l2_error": 1.8059862948884552e-14, "pymhm_pressure_l2_error": 1.0179375343870948e-14, "pymhm_flux_l2_error": 6.201617165443315e-14, "native_divergence_l2": 1.427664339007074e-13, "pymhm_divergence_l2": 3.9126192430035195e-14, "pymhm_physical_residual": 1.563186547454891e-15, "native_metadata": { "equations": 114, "active_equations": 114, "relative_residual": 6.747927791219855e-16 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "driver_sha256": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "binary_sha256": "f051b51711f75646468f527abfa7d452b594f4e7f4fc13e07c520b1cc215818c", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "source_changed": false, "pressure_increment_relative_difference": 2.0603998640110695e-15, "fields_sha256": "b0581b36608d0ee0d04be043e15525affd46c94785fff43f4079248befd64ffb", "native_pressure_l2_norm": 4.143267631552015, "native_pressure_increment_l2_norm": 4.143267631552015, "native_flux_l2_norm": 3.7416573867739404, "pressure_increment_datum": 0.0, "comparison_method": "classical conforming NeoPZ mixed finite elements; identical complete fine trace in PyMHM", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 18, "pressure_space_dimension_per_fine_cell": 4, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "run_compare.py": "7f5a2fd981dc73549b6dcc42c5ad5fc48e4cf548318d90928e4777d280b2dc45" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "degree": 2, "problem": "affine", "fine_factor": 1, "cells": 6, "pressure_l2_difference": 2.402426036770042e-14, "flux_l2_difference": 1.4897996703776225e-13, "pressure_relative_difference": 5.798384874959482e-15, "flux_relative_difference": 3.9816571010584386e-14, "native_pressure_l2_error": 1.537984241712687e-14, "native_flux_l2_error": 1.1201883064923927e-13, "pymhm_pressure_l2_error": 1.6440821981912408e-14, "pymhm_flux_l2_error": 9.395863558535954e-14, "native_divergence_l2": 1.0902199250595184e-12, "pymhm_divergence_l2": 1.778770550868112e-13, "pymhm_physical_residual": 5.934614595173473e-15, "native_metadata": { "equations": 234, "active_equations": 234, "relative_residual": 1.792840510970004e-15 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "driver_sha256": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "binary_sha256": "f051b51711f75646468f527abfa7d452b594f4e7f4fc13e07c520b1cc215818c", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "source_changed": false, "pressure_increment_relative_difference": 5.798384874959482e-15, "fields_sha256": "a0a3718bbc207a0ffb5e8f0e43cefd9576513f5e5392fb10227350cfd4f114d2", "native_pressure_l2_norm": 4.1432676315520185, "native_pressure_increment_l2_norm": 4.1432676315520185, "native_flux_l2_norm": 3.7416573867739418, "pressure_increment_datum": 0.0, "comparison_method": "classical conforming NeoPZ mixed finite elements; identical complete fine trace in PyMHM", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 32, "pressure_space_dimension_per_fine_cell": 10, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "run_compare.py": "7f5a2fd981dc73549b6dcc42c5ad5fc48e4cf548318d90928e4777d280b2dc45" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "prism", "degree": 1, "problem": "well", "fine_factor": 1, "cells": 64, "pressure_l2_difference": 1.195496191185978e-05, "flux_l2_difference": 7.964046184081839e-16, "pressure_relative_difference": 1.8606282947807005e-15, "flux_relative_difference": 2.73992991678088e-13, "native_pressure_l2_error": 24298737.04961014, "native_flux_l2_error": 0.0005812552815531045, "pymhm_pressure_l2_error": 24298737.049607143, "pymhm_flux_l2_error": 0.0005812552815531072, "native_divergence_l2": 7.956187680876976e-16, "pymhm_divergence_l2": 1.5431171209936003e-16, "pymhm_physical_residual": 1.57940009484548e-14, "native_metadata": { "equations": 1760, "active_equations": 1376, "relative_residual": 9.751158412266151e-16 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "driver_sha256": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "binary_sha256": "f051b51711f75646468f527abfa7d452b594f4e7f4fc13e07c520b1cc215818c", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "source_changed": false, "pressure_increment_relative_difference": 3.921517569340792e-14, "fields_sha256": "4aec4beeab87db210eb7b7208133bc735d8ee96de7ae8b44251c5f54fd2f8c7d", "native_pressure_l2_norm": 6425228480.828209, "native_pressure_increment_l2_norm": 304855497.91555345, "native_flux_l2_norm": 0.002906660544602078, "pressure_increment_datum": 25000000.0, "comparison_method": "classical conforming NeoPZ mixed finite elements; identical complete fine trace in PyMHM", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 27, "pressure_space_dimension_per_fine_cell": 6, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "run_compare.py": "7f5a2fd981dc73549b6dcc42c5ad5fc48e4cf548318d90928e4777d280b2dc45" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "degree": 1, "problem": "well", "fine_factor": 1, "cells": 192, "pressure_l2_difference": 2.200318424549883e-05, "flux_l2_difference": 3.622538190458245e-15, "pressure_relative_difference": 3.4255956002331134e-15, "flux_relative_difference": 1.2399944788021123e-12, "native_pressure_l2_error": 19007617.6398407, "native_flux_l2_error": 0.0005028746926763006, "pymhm_pressure_l2_error": 19007617.639847875, "pymhm_flux_l2_error": 0.0005028746926762921, "native_divergence_l2": 7.599683134384704e-16, "pymhm_divergence_l2": 5.512292363356272e-17, "pymhm_physical_residual": 5.8848288977957665e-15, "native_metadata": { "equations": 3312, "active_equations": 2928, "relative_residual": 2.9798555168241843e-15 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "driver_sha256": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "binary_sha256": "f051b51711f75646468f527abfa7d452b594f4e7f4fc13e07c520b1cc215818c", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "source_changed": false, "pressure_increment_relative_difference": 7.164767110556006e-14, "fields_sha256": "18afd71f1ac4b885f4fb5c20e0971b13d982ae82722cf456f7bfd8d9b46cac15", "native_pressure_l2_norm": 6423170395.244991, "native_pressure_increment_l2_norm": 307102574.3890693, "native_flux_l2_norm": 0.0029214147743284887, "pressure_increment_datum": 25000000.0, "comparison_method": "classical conforming NeoPZ mixed finite elements; identical complete fine trace in PyMHM", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 18, "pressure_space_dimension_per_fine_cell": 4, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "run_compare.py": "7f5a2fd981dc73549b6dcc42c5ad5fc48e4cf548318d90928e4777d280b2dc45" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "degree": 2, "problem": "well", "fine_factor": 1, "cells": 192, "pressure_l2_difference": 5.796370483608678e-05, "flux_l2_difference": 4.457133950576556e-15, "pressure_relative_difference": 9.024357142446196e-15, "flux_relative_difference": 1.5250405595046387e-12, "native_pressure_l2_error": 12328251.846725753, "native_flux_l2_error": 0.000495757341552343, "pymhm_pressure_l2_error": 12328251.84671865, "pymhm_flux_l2_error": 0.0004957573415522551, "native_divergence_l2": 3.767636834742085e-15, "pymhm_divergence_l2": 5.092223487771713e-16, "pymhm_physical_residual": 1.5318955668063888e-14, "native_metadata": { "equations": 7152, "active_equations": 6768, "relative_residual": 5.46331008575907e-15 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "driver_sha256": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "binary_sha256": "f051b51711f75646468f527abfa7d452b594f4e7f4fc13e07c520b1cc215818c", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "source_changed": false, "pressure_increment_relative_difference": 1.8828394996944194e-13, "fields_sha256": "2b37a934727e38ab33f669599f5e9bf9121473e0a12d6433c35fc707d33b76ed", "native_pressure_l2_norm": 6423028690.149422, "native_pressure_increment_l2_norm": 307852606.9030003, "native_flux_l2_norm": 0.002922633055755786, "pressure_increment_datum": 25000000.0, "comparison_method": "classical conforming NeoPZ mixed finite elements; identical complete fine trace in PyMHM", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 32, "pressure_space_dimension_per_fine_cell": 10, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "run_compare.py": "7f5a2fd981dc73549b6dcc42c5ad5fc48e4cf548318d90928e4777d280b2dc45" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } } ], "native_mhm_cases": [ { "kind": "prism", "pressure_degree": 1, "problem": "affine", "local_refinement": 2, "macro_cells": 2, "fine_cells": 16, "controller": "TPZMHMixedMeshControl", "native_substructuring": false, "trace_degree": 1, "trace_subdivisions": 1, "native_boundary_quadrature_order": 24, "pymhm_boundary_quadrature_points": 12, "boundary_cap_elimination": "exact EquationFilter of prescribed zero normal-flux connects", "maximum_geometry_difference": 0.0, "pressure_l2_difference": 1.202848722157831e-14, "flux_l2_difference": 1.034727532635838e-13, "pressure_relative_difference": 2.9031402967982033e-15, "flux_relative_difference": 2.7654256541323268e-14, "native_pressure_l2_error": 5.010202508919451e-15, "pymhm_pressure_l2_error": 1.0454773225530034e-14, "native_flux_l2_error": 4.596761247203194e-14, "pymhm_flux_l2_error": 9.056790719653394e-14, "native_divergence_l2": 4.916266576195624e-13, "pymhm_divergence_l2": 1.7780909012724905e-13, "pymhm_physical_residual": 4.660036127378543e-15, "native_metadata": { "equations": 344, "active_equations": 344, "relative_residual": 1.1428938752114386e-15 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "native_binary_sha256": "b3375c7a7d0c66ff90044c84f0a39703c933905b523b9815c9a17da3d1f5e845", "source_changed": false, "fields_sha256": "59273140f3183587e43807202da1239fb272de4f31ee1d995665bfd9f24390ca", "native_pressure_l2_norm": 4.143267631552017, "native_pressure_increment_l2_norm": 4.143267631552017, "native_flux_l2_norm": 3.7416573867739418, "pressure_increment_relative_difference": 2.9031402967982033e-15, "pressure_increment_datum": 0.0, "comparison_method": "native TPZMHMixedMeshControl with RestrainSide; BuildComputationalMesh(false), assembled without substructuring", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 27, "pressure_space_dimension_per_fine_cell": 6, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "dc2a70929c1963090bcd646d3ac857ef74fb64684847ec89f5d7c1d06a79c2a7", "mhm_compare.cpp": "d76cfca75fb679105af340b1e9d300fefa5edea1259a8d874a43943c165bbb79", "run_mhm.py": "e765ccb8ceb780919946db5e3467547ede7d75062e65ff1de1143207ebca6d99" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "pressure_degree": 1, "problem": "affine", "local_refinement": 2, "macro_cells": 6, "fine_cells": 48, "controller": "TPZMHMixedMeshControl", "native_substructuring": false, "trace_degree": 1, "trace_subdivisions": 1, "native_boundary_quadrature_order": 24, "pymhm_boundary_quadrature_points": 12, "boundary_cap_elimination": "exact EquationFilter of prescribed zero normal-flux connects", "maximum_geometry_difference": 0.0, "pressure_l2_difference": 7.898658733208856e-15, "flux_l2_difference": 7.126030738830284e-14, "pressure_relative_difference": 1.9063839065231028e-15, "flux_relative_difference": 1.9045118251658935e-14, "native_pressure_l2_error": 3.308419135383407e-15, "pymhm_pressure_l2_error": 7.866938312249503e-15, "native_flux_l2_error": 2.6130241589135248e-14, "pymhm_flux_l2_error": 6.597892570509045e-14, "native_divergence_l2": 2.1928235242977375e-13, "pymhm_divergence_l2": 1.2854217472353918e-13, "pymhm_physical_residual": 2.084181514711298e-15, "native_metadata": { "equations": 678, "active_equations": 678, "relative_residual": 8.391295754584621e-16 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "native_binary_sha256": "d59ced7b754b67a01231e6d5e948ffbe8e8d677ed4c749e7172953b4d7a7e82e", "source_changed": false, "fields_sha256": "4ba68ae715f4c8ea572ab2533e12af68f2f7ba3c361d8f6bda1dab9b290b3be7", "native_pressure_l2_norm": 4.143267631552016, "native_pressure_increment_l2_norm": 4.143267631552016, "native_flux_l2_norm": 3.74165738677394, "pressure_increment_relative_difference": 1.9063839065231028e-15, "pressure_increment_datum": 0.0, "comparison_method": "native TPZMHMixedMeshControl with RestrainSide; BuildComputationalMesh(false), assembled without substructuring", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 18, "pressure_space_dimension_per_fine_cell": 4, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "mhm_compare.cpp": "c6ddc5b33c897fe5557fbe0e9e7ff529514c5dc5830f6538e818fab761a23bbc", "run_mhm.py": "e765ccb8ceb780919946db5e3467547ede7d75062e65ff1de1143207ebca6d99" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "pressure_degree": 2, "problem": "affine", "local_refinement": 2, "macro_cells": 6, "fine_cells": 48, "controller": "TPZMHMixedMeshControl", "native_substructuring": false, "trace_degree": 1, "trace_subdivisions": 1, "native_boundary_quadrature_order": 24, "pymhm_boundary_quadrature_points": 12, "boundary_cap_elimination": "exact EquationFilter of prescribed zero normal-flux connects", "maximum_geometry_difference": 0.0, "pressure_l2_difference": 1.451401829867005e-14, "flux_l2_difference": 1.5140100641565005e-13, "pressure_relative_difference": 3.5030366342116448e-15, "flux_relative_difference": 4.046362100143757e-14, "native_pressure_l2_error": 9.072697085679622e-15, "pymhm_pressure_l2_error": 1.2937855820053978e-14, "native_flux_l2_error": 9.893894702831338e-14, "pymhm_flux_l2_error": 1.202308198812406e-13, "native_divergence_l2": 1.575696815690043e-12, "pymhm_divergence_l2": 8.929091299301189e-13, "pymhm_physical_residual": 8.072018149720215e-15, "native_metadata": { "equations": 1638, "active_equations": 1638, "relative_residual": 1.837667610746022e-15 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "native_binary_sha256": "d59ced7b754b67a01231e6d5e948ffbe8e8d677ed4c749e7172953b4d7a7e82e", "source_changed": false, "fields_sha256": "5992f9dd9b6dd21406d1a1d72b7afb0ccb4553509bb8e3b4fbac046caf1e0e48", "native_pressure_l2_norm": 4.143267631552021, "native_pressure_increment_l2_norm": 4.143267631552021, "native_flux_l2_norm": 3.741657386773941, "pressure_increment_relative_difference": 3.5030366342116456e-15, "pressure_increment_datum": 0.0, "comparison_method": "native TPZMHMixedMeshControl with RestrainSide; BuildComputationalMesh(false), assembled without substructuring", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 32, "pressure_space_dimension_per_fine_cell": 10, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "5bfe924228a696701f89f9e492c7ddf81735c795e164592e1776cda5736263fb", "mhm_compare.cpp": "c6ddc5b33c897fe5557fbe0e9e7ff529514c5dc5830f6538e818fab761a23bbc", "run_mhm.py": "e765ccb8ceb780919946db5e3467547ede7d75062e65ff1de1143207ebca6d99" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "prism", "pressure_degree": 1, "problem": "well", "local_refinement": 2, "macro_cells": 64, "fine_cells": 512, "controller": "TPZMHMixedMeshControl", "native_substructuring": false, "trace_degree": 1, "trace_subdivisions": 1, "native_boundary_quadrature_order": 24, "pymhm_boundary_quadrature_points": 12, "boundary_cap_elimination": "exact EquationFilter of prescribed zero normal-flux connects", "maximum_geometry_difference": 5.0242958677880805e-15, "pressure_l2_difference": 0.0018457862681482192, "flux_l2_difference": 1.16146461242328e-14, "pressure_relative_difference": 2.872741431824413e-13, "flux_relative_difference": 3.994939986129537e-12, "native_pressure_l2_error": 19068854.08996386, "pymhm_pressure_l2_error": 19068854.0909798, "native_flux_l2_error": 0.0005788264422259387, "pymhm_flux_l2_error": 0.0005788264422267114, "native_divergence_l2": 8.90662091254461e-15, "pymhm_divergence_l2": 3.743668607697109e-16, "pymhm_physical_residual": 1.1546086162771478e-13, "native_metadata": { "equations": 10784, "active_equations": 10400, "relative_residual": 1.5769333092662672e-15 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "native_binary_sha256": "b3375c7a7d0c66ff90044c84f0a39703c933905b523b9815c9a17da3d1f5e845", "source_changed": false, "fields_sha256": "c82a190c0316eaa1ca2e7b7e78546b87e523f3aef49a2dece3634c6b935aa1f2", "native_pressure_l2_norm": 6425173695.413312, "native_pressure_increment_l2_norm": 305776423.7294408, "native_flux_l2_norm": 0.002907339325386349, "pressure_increment_relative_difference": 6.0363917061879835e-12, "pressure_increment_datum": 25000000.0, "comparison_method": "native TPZMHMixedMeshControl with RestrainSide; BuildComputationalMesh(false), assembled without substructuring", "native_linear_solver": "skyline-ldlt", "native_equilibration": "none", "flux_space_dimension_per_fine_cell": 27, "pressure_space_dimension_per_fine_cell": 6, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "dc2a70929c1963090bcd646d3ac857ef74fb64684847ec89f5d7c1d06a79c2a7", "mhm_compare.cpp": "d76cfca75fb679105af340b1e9d300fefa5edea1259a8d874a43943c165bbb79", "run_mhm.py": "e765ccb8ceb780919946db5e3467547ede7d75062e65ff1de1143207ebca6d99" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "pressure_degree": 1, "problem": "well", "local_refinement": 2, "macro_cells": 192, "fine_cells": 1536, "controller": "TPZMHMixedMeshControl", "native_linear_solver": "scipy", "external_solver_diagnostics": { "solver": "SciPy SuperLU with COLAMD and partial pivoting", "equilibration": "five symmetric infinity-row Ruiz passes", "original_residuals": [ 1.4020699420971877e-14 ], "original_matrix_shape": [ 20976, 20976 ], "matrix_nnz": 566446, "original_matrix_sha256": "31804353d1978062e400ba6b695021f0a67bd0cf1b02c90097d4c8a40858f2c2", "original_rhs_sha256": "2090f13f99c3d00ae2d1e5f1c16e8315d735df7fb5c5157081f86ee43ec7cca4" }, "native_equilibration": "five symmetric infinity-row Ruiz passes", "native_substructuring": false, "trace_degree": 1, "trace_subdivisions": 1, "native_boundary_quadrature_order": 24, "pymhm_boundary_quadrature_points": 12, "boundary_cap_elimination": "exact EquationFilter of prescribed zero normal-flux connects", "maximum_geometry_difference": 5.0242958677880805e-15, "pressure_l2_difference": 3.0760244226479904e-05, "flux_l2_difference": 6.7623361111394765e-15, "pressure_relative_difference": 4.788975973984831e-15, "flux_relative_difference": 2.3143980283526466e-12, "native_pressure_l2_error": 13423918.621561598, "pymhm_pressure_l2_error": 13423918.621567318, "native_flux_l2_error": 0.0005004325180426161, "pymhm_flux_l2_error": 0.0005004325180426165, "native_divergence_l2": 7.706054381943742e-16, "pymhm_divergence_l2": 1.6404437522839288e-16, "pymhm_physical_residual": 5.06391615430859e-15, "native_metadata": { "equations": 21360, "active_equations": 20976, "relative_residual": 1.441306353325649e-14 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "native_binary_sha256": "50d6a5cf92241223818d2ff718e302a388b37f866562316a49e3c0ce6dab5179", "source_changed": false, "fields_sha256": "ba16e8b08170ee9750ed6e79df83d1a2764aab221517571344de779d82f78514", "native_pressure_l2_norm": 6423136051.126351, "native_pressure_increment_l2_norm": 307673180.8774473, "native_flux_l2_norm": 0.0029218552851744367, "pressure_increment_relative_difference": 9.997700852168964e-14, "pressure_increment_datum": 25000000.0, "comparison_method": "native TPZMHMixedMeshControl with RestrainSide; BuildComputationalMesh(false), assembled without substructuring", "flux_space_dimension_per_fine_cell": 18, "pressure_space_dimension_per_fine_cell": 4, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "83698cb1e403c9ca08520f6a01060ae284b9901309846a4bbaf94ad4259e9534", "mhm_compare.cpp": "92509a842812b533e185c5c5c5cc78071682f5be22b8a5c5b6cdc24c2540c57a", "run_mhm.py": "64cf23800af09853f4fa25a78f9d0ec58fb9040bed6c0b125773e3169e1ecbde" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } }, { "kind": "tetrahedron", "pressure_degree": 2, "problem": "well", "local_refinement": 2, "macro_cells": 192, "fine_cells": 1536, "controller": "TPZMHMixedMeshControl", "native_linear_solver": "scipy", "external_solver_diagnostics": { "solver": "SciPy SuperLU with COLAMD and partial pivoting", "equilibration": "five symmetric infinity-row Ruiz passes", "original_residuals": [ 2.627906091540224e-14 ], "original_matrix_shape": [ 51696, 51696 ], "matrix_nnz": 2231264, "original_matrix_sha256": "61c71a1e132915b97c57456a4a7c4b95a581975007f266e4b75004b0a52019ab", "original_rhs_sha256": "8fcbaa7dc8d793a5e9e756b6483fe0c5057f037970d3938ec1700493414c5962" }, "native_equilibration": "five symmetric infinity-row Ruiz passes", "native_substructuring": false, "trace_degree": 1, "trace_subdivisions": 1, "native_boundary_quadrature_order": 24, "pymhm_boundary_quadrature_points": 12, "boundary_cap_elimination": "exact EquationFilter of prescribed zero normal-flux connects", "maximum_geometry_difference": 5.0242958677880805e-15, "pressure_l2_difference": 4.5607456135242145e-05, "flux_l2_difference": 8.675895962279729e-15, "pressure_relative_difference": 7.100523934977099e-15, "flux_relative_difference": 2.9691047245111038e-12, "native_pressure_l2_error": 12427556.109177675, "pymhm_pressure_l2_error": 12427556.109189041, "native_flux_l2_error": 0.0004992482382842982, "pymhm_flux_l2_error": 0.0004992482382842746, "native_divergence_l2": 3.282164620543693e-15, "pymhm_divergence_l2": 9.101927204705506e-16, "pymhm_physical_residual": 1.4168827335878446e-14, "native_metadata": { "equations": 52080, "active_equations": 51696, "relative_residual": 2.6885578147249686e-14 }, "native_revision": "4c6b6d277ce097b97bfc8dea1b6725860f4fe05a", "source_hashes": { "src/pymhm/hdiv3d_family.py": "29cd2993848bc6ce5bc944a7a732eae3f9f6cdcc957bf787148d92704dc01c47", "src/pymhm/hdiv3d_mesh.py": "ee7968260eadc10ebf664cd79d5a21d501843cebb24723d12d49ba4119799060", "src/pymhm/darcy_hdiv3d.py": "dcbc14d19771bef9d82e6f4661c6839aa92f20b7be9a82422b217cc26da4d619", "src/pymhm/hybrid.py": "288f85a41eeda946fe7b6a2fc9299c2c9dfc9f1d7bd61a4cc1a625708dae998d", "src/pymhm/solvers.py": "197c1fdcd7c36f695a625bbe0ad0a23c2851aab4b6b945bea3f7b6a4204d7fbc" }, "native_binary_sha256": "50d6a5cf92241223818d2ff718e302a388b37f866562316a49e3c0ce6dab5179", "source_changed": false, "fields_sha256": "d0e29c716b6957944de54fc58e74cd4f9a0d887e728a904791b92093cb87cdb7", "native_pressure_l2_norm": 6423111386.271137, "native_pressure_increment_l2_norm": 307756685.4247571, "native_flux_l2_norm": 0.0029220579155247923, "pressure_increment_relative_difference": 1.481932263219433e-13, "pressure_increment_datum": 25000000.0, "comparison_method": "native TPZMHMixedMeshControl with RestrainSide; BuildComputationalMesh(false), assembled without substructuring", "flux_space_dimension_per_fine_cell": 32, "pressure_space_dimension_per_fine_cell": 10, "quadrature": "Positive tensor/Duffy volume rule, five Gauss points per transformed coordinate; exact for squared discrete polynomial differences.", "native_volume_quadrature": "Material forcing-order policy10 for zero body load; native material chooses order>=11, exceeding the polynomial mass/divergence requirement.", "boundary_quadrature": "Native boundary IntegrationRuleOrder24; PyMHM12 Gauss points per transformed face coordinate; both evaluate the same radial pressure on the planar boundary.", "private_adapter_source_sha256": { "compare.cpp": "83698cb1e403c9ca08520f6a01060ae284b9901309846a4bbaf94ad4259e9534", "mhm_compare.cpp": "92509a842812b533e185c5c5c5cc78071682f5be22b8a5c5b6cdc24c2540c57a", "run_mhm.py": "64cf23800af09853f4fa25a78f9d0ec58fb9040bed6c0b125773e3169e1ecbde" }, "acceptance": { "field_relative_tolerance": 1e-09, "native_algebraic_relative_tolerance": 1e-10, "pymhm_physical_relative_tolerance": 1e-10, "passed": true } } ], "interpretation": "Agreement verifies these discrete spaces, signs, planar geometry, boundary data, native MHM restrictions and physical fields. It is not a reproduction of the historical L05 experiment or a certification of discretization accuracy; analytical errors and classical refinement are reported separately in the case campaign." }