def METHOD_NAME():
    thr = 1.0e-8
    return_result = np.array(
        [
            [-9.15721221e-06, -8.18139252e-06, -6.04628002e-05],
            [+2.63353857e-05, -5.74347326e-05, +6.43925910e-05],
            [-2.99422439e-05, +1.88531187e-05, +3.80203831e-05],
            [-1.05590494e-05, -6.00459729e-05, +3.93481945e-05],
            [+1.59600548e-05, +3.66166973e-05, -2.69939628e-05],
            [-7.21060928e-05, +1.35991320e-05, -3.71000739e-05],
            [-9.01933781e-06, -3.41989101e-05, -4.92317946e-05],
            [-2.38625512e-06, -4.42678339e-05, -1.95513968e-05],
            [-3.09663159e-05, +3.41418638e-05, +2.51926884e-05],
            [+5.60572318e-06, -1.06356845e-05, -3.91159008e-05],
            [+7.73090102e-05, +2.75681060e-05, +2.02933984e-05],
            [-3.23109274e-05, +5.39319105e-05, -1.44309772e-04],
            [-4.05785623e-05, +5.03251549e-05, +3.92193348e-05],
            [+8.46365844e-06, -3.35282588e-05, +3.76937643e-05],
            [+1.04861677e-04, -2.99444252e-05, +8.20822571e-05],
            [-1.50951292e-06, +4.32012272e-05, +3.05230899e-05],
        ]
    )
    atomic_input = qcel.models.AtomicInput(
        molecule=qcng.get_molecule("mindless-03"),
        keywords={"level_hint": "d3bj"},
        driver="gradient",
        model={"method": "r2scan"},
    )
    atomic_result = qcng.compute(atomic_input, "s-dftd3")
    assert atomic_result.success
    assert pytest.approx(return_result, abs=thr) == atomic_result.return_result
    assert pytest.approx(atomic_result.return_result, abs=thr) == return_result