Function af_pinverse

src/api/c/pinverse.cpp:170–201

}

af_err af_pinverse(af_array *out, const af_array in, const double tol,
                   const af_mat_prop options) {
    try {
        const ArrayInfo &i_info = getInfo(in);

        af_dtype type = i_info.getType();

        if (options != AF_MAT_NONE) {
            AF_ERROR("Using this property is not yet supported in inverse",
                     AF_ERR_NOT_SUPPORTED);
        }

        ARG_ASSERT(1, i_info.isFloating());  // Only floating and complex types
        ARG_ASSERT(2, tol >= 0.);            // Ensure tolerance is not negative

        af_array output;

        if (i_info.ndims() == 0) { return af_retain_array(out, in); }

        switch (type) {
            case f32: output = pinverse<float>(in, tol); break;
            case f64: output = pinverse<double>(in, tol); break;
            case c32: output = pinverse<cfloat>(in, tol); break;
            case c64: output = pinverse<cdouble>(in, tol); break;
            default: TYPE_ERROR(1, type);
        }
        swap(*out, output);
    }
    CATCHALL;

    return AF_SUCCESS;
}