Transforming over 1000 vertices to C_CONTIGUOUS. Transforming over 1000 elements to C_CONTIGUOUS. ================================================================================ OpenSees three-way cross-check: case=gravity, site=baseline ================================================================================ Mesh: reports/full_dome_merged.vtu N_nodes = 20,187 N_tets = 57,725 reorient_tets: flipped 17/57725 tets (0.0%) Clamp NSET (y < 50.0 mm): 163 nodes Material: E=70.80 MPa, nu=0.3, rho=240.0 kg/m^3 Pressure: +0 Pa (sign convention: + = compressive on outer face) skfem-convention value: +0 Pa (+ = outward push) --- scikit-fem solve --- [skfem] 3.0 s |u|max=30.002 |u|p99=11.926 VMmax=4456.2807 VMp99=1.5433 sigma_t_max=+92991.6617 sigma_c_min=-150.7437 --- CalculiX solve --- [ccx] 1.1 s |u|max=6.886 |u|p99=6.886* VMmax=0.1488 VMp99=0.0472 sigma_t_max=+0.1476 sigma_c_min=-0.1809 (* ccx wrapper reports only scalar |u|max; p99 column duplicates max for this solver) --- OpenSees solve --- [opensees] 1.0 s |u|max=6.886 |u|p99=6.138 VMmax=0.1488 VMp99=0.0472 sigma_t_max=+0.1476 sigma_c_min=-0.1809 ============ AGREEMENT ============ metric skfem ccx opensees max diff % ---------------------------------------------------------------------- |u| p99 [mm] 11.926 6.886 6.138 48.5 |u| max [mm] 30.002 6.886 6.886 77.0 VM p99 [MPa] 1.5433 0.0472 0.0472 96.9 sigma_t max [MPa] 92991.6617 0.1476 0.1476 100.0 sigma_c min [MPa] -150.7437 -0.1809 -0.1809 99.9 VERDICT: FAIL (criteria: |u|p99 diff <= 5%, VM p99 diff <= 10%) - |u| p99 max-diff: 48.5 % [FAIL] - VM p99 max-diff: 96.9 % [FAIL] - sigma sign agreement: [yes] (sigma_t: skfem=+, ccx=+, opensees=+; sigma_c: skfem=-, ccx=-, opensees=-) Note: FAIL on this mesh is the documented status quo (NEXT_STEPS.md #1: skfem-vs-ccx p99 VM disagree by ~90% on the merged mesh, likely sliver-tet artefacts at joints). OpenSees acts as tiebreaker. The table above is the value; this tool intentionally does NOT abort on FAIL.