Ticket #254: polynomialBug.c

#include <stdio.h>
#include "pslib.h"

#define ORDER 3
#define BAD 1

int main(int argc, char **argv)
{
    psPolynomial1D *poly = psPolynomial1DAlloc(ORDER, PS_POLYNOMIAL_ORD);

    for (int i = 0; i < poly->n; i++) {
        poly->coeff[i] = (float)i;
        if (i == BAD) {
            poly->mask[i] = 1;  // Don't use this coefficient!
        } else {
            poly->mask[i] = 0;
        }
    }

    fprintf(stderr,"Value using mask: %f\n", psPolynomial1DEval(0.5, poly));

    // A different method of doing the same thing
    poly->coeff[BAD] = 0.0;             // Set coefficient to zero instead of masking
    poly->mask[BAD] = 0;                // Use this coefficient

    fprintf(stderr,"Value using coefficient zeroing: %f\n", psPolynomial1DEval(0.5, poly));

    psFree(poly);

    exit(EXIT_SUCCESS);
}