def METHOD_NAME():
    thr = 1.0e-8
    return_result = np.array(
        [
            [-1.0184389290110640e-05, -1.2902131714706869e-05, -4.7991389860088184e-05],
            [+3.9008038585852290e-05, -8.6543337237836941e-05, +1.0086268708372256e-04],
            [-2.1792111331028049e-05, +2.8414264003578678e-05, +1.8916279332120321e-05],
            [-3.8847190487495055e-06, +6.4069296782535270e-06, -1.7362611928711208e-05],
            [-1.5358378150451566e-05, -5.7611513737809471e-06, -1.9871764856476196e-05],
            [-3.7225957329162548e-05, -1.5687552833166010e-06, -3.4717349142275990e-05],
            [+3.2408291983641028e-07, -6.7683941908276048e-05, -6.1602680430568183e-05],
            [-8.1819149671629329e-06, -9.5470867672313039e-06, -1.7007294079721094e-05],
            [-1.7245516321582448e-05, +1.6416848978506148e-05, +1.6005589814916458e-05],
            [-3.0116095523698349e-06, +7.3382310356380993e-07, -1.3652036487357768e-05],
            [+1.0378663330288130e-04, +4.4180555120489236e-05, +3.2576724474326750e-05],
            [-2.8021453980080289e-05, +3.9008636389849655e-05, -1.2083469461741150e-04],
            [-3.0452136766248191e-05, +7.0712047105794630e-05, +8.3682948356528971e-05],
            [+1.6902554427499383e-06, -2.9147110970219634e-05, +3.6402009185785855e-05],
            [+4.6966237842630462e-05, -1.2554523547742754e-05, +3.0600183501647347e-05],
            [-1.6417061357004406e-05, +1.9834934423075429e-05, +1.3993399653561848e-05],
        ]
    )
    atomic_input = qcel.models.AtomicInput(
        molecule=qcng.get_molecule("mindless-03"),
        keywords={"level_hint": "D4"},
        driver="gradient",
        model={"method": "r2scan"},
    )
    atomic_result = qcng.compute(atomic_input, "dftd4")
    print(atomic_result.return_result)
    assert atomic_result.success
    assert pytest.approx(atomic_result.return_result, abs=thr) == return_result