import { eEstimate } from "big-float-ts";


/**
 * Returns the `p-infinity norm`, i.e. the maximum magnitude absolute value 
 * within the given array of numbers / coefficients (with intermediate 
 * calculations (and the final result) done in double precision).
 * 
 * @param p an array of numbers; can represent an array of polynomial 
 * coefficients
 * 
 * @doc
 */
function ePInfNorm(p: number[][]): number {
    let max = 0;
    for (let i=0; i<p.length; i++) {
        const v = Math.abs(eEstimate(p[i]));
        if (v > max) { 
            max = v; 
        }
    }
    
    return max;
}


export { ePInfNorm }