Function pj_inv

src/inv.cpp:153–186

}

PJ_LP pj_inv(PJ_XY xy, PJ *P) {
    PJ_COORD coo = {{0, 0, 0, 0}};
    coo.xy = xy;

    const int last_errno = P->ctx->last_errno;
    P->ctx->last_errno = 0;

    if (!P->skip_inv_prepare)
        inv_prepare(P, coo);
    if (HUGE_VAL == coo.v[0])
        return proj_coord_error().lp;

    /* Do the transformation, using the lowest dimensional transformer available
     */
    if (P->inv) {
        const auto lp = P->inv(coo.xy, P);
        coo.lp = lp;
    } else if (P->inv3d) {
        const auto lpz = P->inv3d(coo.xyz, P);
        coo.lpz = lpz;
    } else if (P->inv4d)
        P->inv4d(coo, P);
    else {
        proj_errno_set(P, PROJ_ERR_OTHER_NO_INVERSE_OP);
        return proj_coord_error().lp;
    }
    if (HUGE_VAL == coo.v[0])
        return proj_coord_error().lp;

    if (!P->skip_inv_finalize)
        inv_finalize(P, coo);

    return error_or_coord(P, coo, last_errno).lp;
}

PJ_LPZ pj_inv3d(PJ_XYZ xyz, PJ *P) {
    PJ_COORD coo = {{0, 0, 0, 0}};