def METHOD_NAME(out, a1, a2, h0, mtx_c, c33, mtx_b, mtx_t, geo,
                  term_mode, fmode):
    if term_mode == 'strain':
        out_qp = membranes.get_green_strain_sym3d(mtx_c, c33)
    elif term_mode == 'stress':
        n_el, n_qp, dm, _ = mtx_c.shape
        dim = dm + 1
        sym = dim2sym(dim)
        out_qp = nm.zeros((n_el, n_qp, sym, 1), dtype=mtx_c.dtype)
        stress = eval_membrane_mooney_rivlin(a1, a2, mtx_c, c33, 0)
        out_qp[..., 0:2, 0] = stress[..., 0:2, 0]
        out_qp[..., 3, 0] = stress[..., 2, 0]
    status = geo.integrate(out, out_qp, fmode)
    out[:, 0, :, 0] = transform_data(out.squeeze(), mtx=mtx_t)
    return status