Discipline contract

The interface a multidisciplinary framework integrates against is checked on a small inviscid aerofoil case (NACA 0012, M = 0.5, α = 2°, 64 elements of degree 2). (i) The partial-derivative products of residual and outputs (apply_linear) in forward and reverse mode must be transposes of each other, with respect to the state and to all inputs including the surface displacement field, both at the discipline's solution and at a point named by the caller that is not a solution. (ii) The forward products must agree with central finite differences of the residual (apply_nonlinear) and of the outputs, taken along smooth directions (the numerical flux is only piecewise smooth, so a rough direction would test the finite difference and not the derivative). (iii) solve_linear must invert apply_linear in both modes, and the two modes must be transposes of each other. (iv) The reduced (total) products tangent_product and adjoint_product must be transposes of each other. (v) A linear solve that cannot reach its tolerance must raise, or be flagged when the caller asks for that; and adjoint and direct totals, which solve different linear systems of the same linearisation, must approach each other as the tolerance of those solves is tightened. (vi) Wrapped as an OpenMDAO implicit component, the discipline must pass OpenMDAO's own check of partial derivatives, reproduce the discipline's total derivatives, and, coupled to a linear structure in a block Gauss–Seidel cycle, give a coupled total derivative that agrees with a finite difference of re-converged coupled analyses.

• flow solve convergedyesexpected truepass
• apply_linear forward vs reverse: largest relative dot-product defect over both points2.51e-16expected < 1e-11 (round-off)pass
• solve_linear forward vs reverse: largest relative dot-product defect4.05e-12expected < 1e-8 (two Krylov solves to 1e-11)pass
• solve_linear inverts apply_linear in both modes: largest relative residual, recomputed with apply_linear9.55e-12expected < 1e-9 (solves to 1e-11; preconditioned residual)pass
• partials (residual, cl, cd, cm, f_surface) with respect to state and inputs vs central differences: largest relative difference1.15e-08expected < 1e-6 (the bound of the linearization study for JVP vs central differences; smooth directions)pass
• total products, tangent_product vs adjoint_product: relative dot-product defect6.64e-10expected < 1e-7 (two linear solves to 1e-10)pass
• a linear solve that misses its tolerance raises LinearSolveError and reports the reduction it achieved1.65e-14expected raised, with 0 < reduction < 1e-10pass
• with on_failure='flag' the same solve returns with info.convergednoexpected falsepass
• adjoint and direct totals meet as the linear tolerance tightens: largest relative difference at 1e-8, 1e-10, 1e-122.9e-05 / 1.3e-06 / 5.8e-09expected decreasing, below 1e-7 at 1e-12pass
• OpenMDAO is installed (extra 'mdo')yesexpected truepass
• OpenMDAO check_partials of the implicit flow component and the functionals component (forward and reverse vs central differences): largest relative error5.65e-09expected < 1e-6pass
• OpenMDAO totals through the implicit component vs the discipline's adjoint totals (rel., largest)4.14e-16expected < 1e-6 (linear solves to 1e-10)pass
• coupled flow–structure adjoint through the implicit component vs finite difference (rel.)1.23e-07expected < 5e-6 (finite-difference step 1e-4)pass
• coupled totals, forward vs reverse mode (rel.)9.71e-10expected < 1e-6 (linear block Gauss–Seidel to 1e-9)pass
• the coupling is active: largest deflection [chords]0.005555expected > 1e-5pass

Run notes OpenMDAO 3.45.2. The coupled derivative differs from the rigid one by +0.31 %. Measured in this run: largest dot-product defect of apply_linear 2.5e-16, of solve_linear 4.1e-12, of the total products 6.6e-10; largest difference between partials and central differences 1.2e-08; the solve with an unreachable tolerance stopped at a reduction of 1.6e-14.

gyre 0.4.0 · 39979fcf90 · 2026-11-12 11:40 · cpu · quick suite · record (JSON)Same code as this build