QuadraticFitTriggerInterpolantTest.c

Go to the documentation of this file.

/*
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with with program; see the file COPYING. If not, write to the
 * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
 * MA  02110-1301  USA
 *
 * Copyright (C) 2012 Leo Singer
 */
 
 
#include <complex.h>
#include <stdlib.h>
 
#include <lal/TriggerInterpolation.h>
 
int main(__attribute__ ((unused)) int argc, __attribute__ ((unused)) char **argv)
{
    int result;
    double tmax;
 
    QuadraticFitTriggerInterpolant *interp = XLALCreateQuadraticFitTriggerInterpolant(2);
    if (!interp)
        exit(EXIT_FAILURE);
 
    {
        const COMPLEX16 y[] = {0, 3, 4, 3, 0};
        COMPLEX16 ymax;
 
        result = XLALCOMPLEX16ApplyQuadraticFitTriggerInterpolant(interp, &tmax, &ymax, &y[2]);
        if (result)
            exit(EXIT_FAILURE);
 
        if (tmax != 0)
            exit(EXIT_FAILURE);
        if (creal(ymax) != 4)
            exit(EXIT_FAILURE);
        if (cimag(ymax) != 0)
            exit(EXIT_FAILURE);
    }
 
    {
        const COMPLEX8 y[] = {0, 3, 4, 3, 0};
        COMPLEX8 ymax;
 
        result = XLALCOMPLEX8ApplyQuadraticFitTriggerInterpolant(interp, &tmax, &ymax, &y[2]);
        if (result)
            exit(EXIT_FAILURE);
 
        if (tmax != 0)
            exit(EXIT_FAILURE);
        if (crealf(ymax) != 4)
            exit(EXIT_FAILURE);
        if (cimagf(ymax) != 0)
            exit(EXIT_FAILURE);
    }
 
    {
        const REAL8 y[] = {0, 3, 4, 3, 0};
        REAL8 ymax;
 
        result = XLALREAL8ApplyQuadraticFitTriggerInterpolant(interp, &tmax, &ymax, &y[2]);
        if (result)
            exit(EXIT_FAILURE);
 
        if (tmax != 0)
            exit(EXIT_FAILURE);
        if (ymax != 4)
            exit(EXIT_FAILURE);
    }
 
    {
        const REAL4 y[] = {0, 3, 4, 3, 0};
        REAL4 ymax;
 
        result = XLALREAL4ApplyQuadraticFitTriggerInterpolant(interp, &tmax, &ymax, &y[2]);
        if (result)
            exit(EXIT_FAILURE);
 
        if (tmax != 0)
            exit(EXIT_FAILURE);
        if (ymax != 4)
            exit(EXIT_FAILURE);
    }
 
    XLALDestroyQuadraticFitTriggerInterpolant(interp);
    exit(EXIT_SUCCESS);
}