Function denominator

src/api/c/deconvolution.cpp:241–263

template<typename CT>
Array<CT> denominator(const Array<CT>& I, const Array<CT>& P, const float gamma,
                      const af_inverse_deconv_algo algo) {
    using T = typename af::dtype_traits<CT>::base_type;

    auto RCNST = createValueArray(I.dims(), scalar<T>(gamma));

    if (algo == AF_INVERSE_DECONV_TIKHONOV) {
        auto normP = complexNorm<T, CT>(P);
        auto denom = arithOp<T, af_add_t>(normP, RCNST, normP.dims());

        return cast<CT, T>(denom);
    } else {
        // TODO(pradeep) Wiener Filter code path is disabled.
        // This code path doesn't is not exposed using current API
        auto normI = complexNorm<T, CT>(I);
        auto sRes  = arithOp<T, af_sub_t>(normI, RCNST, normI.dims());
        auto dRes  = arithOp<T, af_div_t>(RCNST, sRes, RCNST.dims());
        auto normP = complexNorm<T, CT>(P);
        auto denom = arithOp<T, af_add_t>(normP, dRes, normP.dims());

        return cast<CT, T>(denom);
    }
}

template<typename InputType, typename RealType = float>
af_array invDeconv(const af_array in, const af_array ker, const float gamma,