/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:30:08 */ /*FOR_C Options SET: ftn=u io=c no=p op=aimnv s=dbov str=l x=f - prototypes */ #include #include "fcrt.h" #include "sconcm.h" #include void /*FUNCTION*/ sconcm( long n, float coeff[]) { long int i, j; float tp; /* Copyright (c) 1996 California Institute of Technology, Pasadena, CA. * ALL RIGHTS RESERVED. * Based on Government Sponsored Research NAS7-03001. *>> 1994-11-11 SCONCM Krogh Declared all vars. *>> 1994-10-20 SCONCM Krogh Changes to use M77CON *>> 1992-06-04 SCONCM Krogh Corrected line below. *>> 1992-01-31 SCONCM Krogh Initial Code. *--S replaces "?": ?CONCM * * Converts COEFF in Chebyshev basis to the monomial basis. * */ if (n <= 0) return; tp = 1.e0; for (j = 0; j <= (n - 2); j++) { for (i = n - 2; i >= j; i--) { coeff[i] -= coeff[i + 2]; } coeff[j + 1] *= .5e0; coeff[j] *= tp; tp *= 2.e0; } coeff[n] *= tp; coeff[n - 1] *= tp; return; } /* end of function */