ptwXY_misc.c File Reference

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
#include "ptwXY.h"

Go to the source code of this file.

Functions
double	ptwXY_limitAccuracy (double accuracy)
void	ptwXY_update_biSectionMax (ptwXYPoints *ptwXY1, double oldLength)
ptwXYPoints *	ptwXY_createFromFunction (statusMessageReporting *smr, int n, double *xs, ptwXY_createFromFunction_callback func, void *argList, double accuracy, int checkForRoots, int biSectionMax)
ptwXYPoints *	ptwXY_createFromFunction2 (statusMessageReporting *smr, ptwXPoints *xs, ptwXY_createFromFunction_callback func, void *argList, double accuracy, int checkForRoots, int biSectionMax)
nfu_status	ptwXY_applyFunction (statusMessageReporting *smr, ptwXYPoints *ptwXY1, ptwXY_applyFunction_callback func, void *argList, int checkForRoots)
ptwXYPoints *	ptwXY_fromString (statusMessageReporting *smr, char const *str, char sep, ptwXY_interpolation interpolation, char const *interpolationString, double biSectionMax, double accuracy, char **endCharacter, int useSystem_strtod)
void	ptwXY_showInteralStructure (ptwXYPoints *ptwXY, FILE *f, int printPointersAsNull)
void	ptwXY_simpleWrite (ptwXYPoints *ptwXY, FILE *f, char const *format)
void	ptwXY_simplePrint (ptwXYPoints *ptwXY, char const *format)

Function Documentation

ptwXY_applyFunction()

nfu_status ptwXY_applyFunction	(	statusMessageReporting *	smr,
		ptwXYPoints *	ptwXY1,
		ptwXY_applyFunction_callback	func,
		void *	argList,
		int	checkForRoots )

Definition at line 165 of file ptwXY_misc.c.

                                           {
 
    int64_t i, originalLength = ptwXY1->length, notFirstPass = 0;
    double y1, y2 = 0;
    nfu_status status;
    ptwXYPoint p1, p2;
 
    checkForRoots = checkForRoots && ptwXY1->biSectionMax;
 
    if( ptwXY1->interpolation == ptwXY_interpolationOther ) {
        smr_setReportError2p( smr, nfu_SMR_libraryID, nfu_otherInterpolation, "Other interpolation not allowed." );
        return( ptwXY1->status = nfu_otherInterpolation );
    }
    if( ptwXY1->interpolation == ptwXY_interpolationFlat ) {
        smr_setReportError2p( smr, nfu_SMR_libraryID, nfu_flatInterpolation, "Flat interpolation not allowed." );
        return( ptwXY1->status = nfu_flatInterpolation );
    }
 
    if( ptwXY_simpleCoalescePoints( smr, ptwXY1 ) != nfu_Okay ) goto Err;
    for( i = originalLength - 1; i >= 0; i-- ) {
        y1 = ptwXY1->points[i].y;
        if( ( status = func( smr, &(ptwXY1->points[i]), argList ) ) != nfu_Okay ) {
            if( ptwXY1->status == nfu_Okay ) ptwXY1->status = status;
            return( status );
        }
        p1 = ptwXY1->points[i];
        if( notFirstPass ) {
            if( ptwXY_applyFunction2( smr, ptwXY1, y1, y2, &p1, &p2, func, argList, 0, checkForRoots ) != nfu_Okay ) goto Err;
        }
        notFirstPass = 1;
        p2 = p1;
        y2 = y1;
    }
    ptwXY_update_biSectionMax( ptwXY1, (double) originalLength );
    return( status );
 
Err:
    smr_setReportError2p( smr, nfu_SMR_libraryID, nfu_Error, "Via." );
    if( ptwXY1->status == nfu_Okay ) ptwXY1->status = nfu_Error;
    return( nfu_Error );
}

Referenced by ptwXY_pow().

ptwXY_createFromFunction()

ptwXYPoints * ptwXY_createFromFunction	(	statusMessageReporting *	smr,
		int	n,
		double *	xs,
		ptwXY_createFromFunction_callback	func,
		void *	argList,
		double	accuracy,
		int	checkForRoots,
		int	biSectionMax )

Definition at line 46 of file ptwXY_misc.c.

                                                                              {
 
    int64_t i;
    double x1, y1, x2, y2, eps = ClosestAllowXFactor * DBL_EPSILON;
    ptwXYPoints *ptwXY;
    ptwXYPoint *p1, *p2;
 
    if( n < 2 ) {
        smr_setReportError2( smr, nfu_SMR_libraryID, nfu_tooFewPoints, "Too few point = %d.", (int) n );
        return( NULL );
    }
    for( i = 1; i < n; i++ ) {
        if( xs[i-1] >= xs[i] ) {
            smr_setReportError2( smr, nfu_SMR_libraryID, nfu_XNotAscending, 
                    "Non-ascending domain values: xs[%d] = %.17e >= xs[%d] = %.17e.", 
                    (int) (i-1), xs[i-1], (int) i, xs[i] );
            return( NULL );
        }
    }
 
    x1 = xs[0];
    if( func( smr, x1, &y1, argList ) != nfu_Okay ) {
        return( NULL );
    }
    if( ( ptwXY = ptwXY_new( smr, ptwXY_interpolationLinLin, NULL, biSectionMax, accuracy, 500, 50, 0 ) ) == NULL ) goto Err;
    for( i = 1; i < n; i++ ) {
        if( ptwXY_setValueAtX_overrideIfClose( smr, ptwXY, x1, y1, eps, 0 ) != nfu_Okay ) goto Err;
        x2 = xs[i];
        if( func( smr, x2, &y2, argList ) != nfu_Okay ) goto Err;
        if( ptwXY_createFromFunctionBisect( smr, ptwXY, x1, y1, x2, y2, func, argList, 0, checkForRoots, eps ) != nfu_Okay ) goto Err;
        x1 = x2;
        y1 = y2;
    }
    if( ptwXY_setValueAtX_overrideIfClose( smr, ptwXY, x2, y2, eps, 1 ) != nfu_Okay ) goto Err;
 
    if( checkForRoots ) {
        if( ptwXY_simpleCoalescePoints( NULL, ptwXY ) != nfu_Okay ) goto Err;
        for( i = ptwXY->length - 1, p2 = NULL; i >= 0; i--, p2 = p1 ) { /* Work backward so lower points are still valid if a new point is added. */
            p1 = &(ptwXY->points[i]);
            if( p2 != NULL ) {
                if( ( p1->y * p2->y ) < 0. ) {
                    if( ptwXY_createFromFunctionZeroCrossing( smr, ptwXY, p1->x, p1->y, p2->x, p2->y, func, argList, eps ) != nfu_Okay ) goto Err;
                }
            }
        }
    }
 
    return( ptwXY );
 
Err:
    smr_setReportError2p( smr, nfu_SMR_libraryID, nfu_Error, "Via." );
    if( ptwXY != NULL ) ptwXY_free( ptwXY );
    return( NULL );
}

Referenced by GIDI::Functions::Polynomial1d::asXYs1d(), MCGIDI::Probabilities::NBodyPhaseSpace2d::NBodyPhaseSpace2d(), nf_Legendre_to_ptwXY(), and ptwXY_createFromFunction2().

ptwXY_createFromFunction2()

ptwXYPoints * ptwXY_createFromFunction2	(	statusMessageReporting *	smr,
		ptwXPoints *	xs,
		ptwXY_createFromFunction_callback	func,
		void *	argList,
		double	accuracy,
		int	checkForRoots,
		int	biSectionMax )

Definition at line 104 of file ptwXY_misc.c.

                                                                              {
 
    ptwXYPoints *ptwXY = ptwXY_createFromFunction( smr, (int) xs->length, xs->points, func, argList, accuracy, 
            checkForRoots, biSectionMax );
 
    if( ptwXY == NULL ) smr_setReportError2p( smr, nfu_SMR_libraryID, nfu_Error, "Via." );
    return( ptwXY );
}

ptwXY_fromString()

ptwXYPoints * ptwXY_fromString	(	statusMessageReporting *	smr,
		char const *	str,
		char	sep,
		ptwXY_interpolation	interpolation,
		char const *	interpolationString,
		double	biSectionMax,
		double	accuracy,
		char **	endCharacter,
		int	useSystem_strtod )

Definition at line 275 of file ptwXY_misc.c.

                                                                                                                           {
 
    int64_t numberConverted;
    double  *doublePtr;
    ptwXYPoints *ptwXY = NULL;
 
    if( ( doublePtr = nfu_stringToListOfDoubles( smr, str, sep, &numberConverted, endCharacter, useSystem_strtod ) ) == NULL ) {
        smr_setReportError2p( smr, nfu_SMR_libraryID, nfu_Error, "Via." );
        return( NULL );
    }
    if( ( numberConverted % 2 ) == 0 ) {
        ptwXY = ptwXY_create( smr, interpolation, interpolationString, biSectionMax, accuracy, numberConverted / 2, 10, numberConverted / 2, doublePtr, 0 ); }
    else {
        smr_setReportError2( smr, nfu_SMR_libraryID, nfu_oddNumberOfValues, "Odd number = %d of float for ptwXY.", (int) numberConverted );
    }
    smr_freeMemory2( doublePtr );
    return( ptwXY );
}

ptwXY_limitAccuracy()

double ptwXY_limitAccuracy

(

double

accuracy

)

Definition at line 28 of file ptwXY_misc.c.

                                              {
 
    if( accuracy < ptwXY_minAccuracy ) accuracy = ptwXY_minAccuracy;
    if( accuracy > 1 ) accuracy = 1.;
    return( accuracy );
}

Referenced by ptwXY_setAccuracy(), ptwXY_thin(), and ptwXY_toOtherInterpolation().

ptwXY_showInteralStructure()

void ptwXY_showInteralStructure	(	ptwXYPoints *	ptwXY,
		FILE *	f,
		int	printPointersAsNull )

Definition at line 297 of file ptwXY_misc.c.

                                                                                        {
 
    int64_t i, n1;
    ptwXYPoint *point = ptwXY->points;
    ptwXYOverflowPoint *overflowPoint;
 
    n1 = ptwXY_getNonOverflowLength( NULL, ptwXY );
 
    fprintf( f, "status = %d  interpolation = %d  length = %d  allocatedSize = %d\n", 
        (int) ptwXY->status, (int) ptwXY->interpolation, (int) ptwXY->length, (int) ptwXY->allocatedSize );
    fprintf( f, "userFlag = %d  biSectionMax = %.8e  accuracy = %.2e  minFractional_dx = %.6e\n", 
        ptwXY->userFlag, ptwXY->biSectionMax, ptwXY->accuracy, ptwXY->minFractional_dx );
    fprintf( f, "interpolationString = %s\n", ptwXY->interpolationString );
    fprintf( f, "  overflowLength = %d  overflowAllocatedSize = %d  mallocFailedSize = %d\n", 
        (int) ptwXY->overflowLength, (int) ptwXY->overflowAllocatedSize, (int) ptwXY->mallocFailedSize );
    fprintf( f, "  Points data, points = %20p\n", ( printPointersAsNull ? NULL : ptwXY->points ) );
    for( i = 0; i < n1; i++,  point++ ) fprintf( f, "    %14.7e %14.7e\n", point->x, point->y );
    fprintf( f, "  Overflow points data; %20p\n", ( printPointersAsNull ? NULL : &(ptwXY->overflowHeader) ) );
    for( overflowPoint = ptwXY->overflowHeader.next; overflowPoint != &(ptwXY->overflowHeader); overflowPoint = overflowPoint->next ) {
        fprintf( f, "    %14.7e %14.7e %8d %20p %20p %20p\n", overflowPoint->point.x, overflowPoint->point.y, (int) overflowPoint->index, 
            ( printPointersAsNull ? NULL : overflowPoint ), ( printPointersAsNull ? NULL : overflowPoint->prior ), 
            ( printPointersAsNull ? NULL : overflowPoint->next ) );
    }
    fprintf( f, "  Points in order\n" );
    for( i = 0; i < ptwXY->length; i++ ) {
        point = ptwXY_getPointAtIndex_Unsafely( ptwXY, i );
        fprintf( f, "    %14.7e %14.7e\n", point->x, point->y );
    }
}

ptwXY_simplePrint()

void ptwXY_simplePrint	(	ptwXYPoints *	ptwXY,
		char const *	format )

Definition at line 342 of file ptwXY_misc.c.

                                                                 {
 
    ptwXY_simpleWrite( ptwXY, stdout, format );
}

ptwXY_simpleWrite()

void ptwXY_simpleWrite	(	ptwXYPoints *	ptwXY,
		FILE *	f,
		char const *	format )

Definition at line 329 of file ptwXY_misc.c.

                                                                          {
 
    int64_t i;
    ptwXYPoint *point;
 
    for( i = 0; i < ptwXY->length; i++ ) {
        point = ptwXY_getPointAtIndex_Unsafely( ptwXY, i );
        fprintf( f, format, point->x, point->y );
    }
}

Referenced by ptwXY_simplePrint(), and GIDI::Functions::XYs1d::write().

ptwXY_update_biSectionMax()

void ptwXY_update_biSectionMax	(	ptwXYPoints *	ptwXY1,
		double	oldLength )

Definition at line 37 of file ptwXY_misc.c.

                                                                        {
 
    ptwXY1->biSectionMax = ptwXY1->biSectionMax - 1.442695 * log( ptwXY1->length / oldLength ); /* 1.442695 = 1 / log( 2. ) */
    if( ptwXY1->biSectionMax < 0 ) ptwXY1->biSectionMax = 0;
    if( ptwXY1->biSectionMax > ptwXY_maxBiSectionMax ) ptwXY1->biSectionMax = ptwXY_maxBiSectionMax;
}

Referenced by ptwXY_applyFunction(), ptwXY_div_ptwXY(), and ptwXY_mul2_ptwXY().