```diff + iszero(analytical) ? isapprox(analytical, numerical, atol=1e-5) : isapprox(analytical, numerical, rtol=1e-5) - isapprox(analytical, numerical, rtol=1e-5, atol=1e-5) ```