Category : C Source Code
Archive   : CEPHES22.ZIP
Filename : REMESW.C

 
Output of file : REMESW.C contained in archive : CEPHES22.ZIP
/* remesw.c */

/* Write answers to disc file. */

#include
#include "remes.h"

remesw()
{
char s[40];
int i;
FILE *f, *fopen();
double log10();

reopn:

printf( "Log file name ? " );
gets(s);
f = fopen( s, "w" );
if( f == 0 )
{
printf( "Can't open %s\n", s );
goto reopn;
}

fprintf( f, "\n%s = ", funnam );
if( config & CW )
fprintf( f, "x + x**2 * " );
if( config & XPX )
fprintf( f, "z " );
if( config & X2PX )
fprintf( f, "z**2 " );
if( config & PXCU )
fprintf( f, "P(z**3)" );
else if( config & PXSQ )
fprintf( f, "P(z**2)" );
else
fprintf( f, "P(z)" );
if( d > 0 )
{
if( config & PXCU )
fprintf( f, "/Q(z**2)" );
else if( config & PXSQ )
fprintf( f, "/Q(z**2)" );
else
fprintf( f, "/Q(z)" );
}
fprintf( f, "\nz(x) = %s\n", znam );
if( relerr )
fprintf( f, "Relative error\n" );
else
fprintf( f, "Absolute error\n" );

fprintf( f, "n = %d, d = %d\n", n, d );
fprintf( f, "precision = %23.15E ", -log10(farther) );
fprintf( f, "error=%23.15E\n", farther );
fprintf( f, "leveled=%23.4E\n", spread );
fprintf( f, "\nNumerator:\n" );
for( i=0; i {
if( i == n+1 )
fprintf( f, "\nDenominator\n" );
fprintf( f, "%23.15E\n", param[i] );
}

if( d > 0 )
fprintf( f, "%23.15E\n", 1.0 );

if( (d > 0) || ((config & ZER) != 0) )
{
fprintf( f, "\n Locations of zero error:\n" );
for( i=0; i<=n+d; i++ )
fprintf( f, "%23.15E\n", xx[i] );
}

fprintf( f, "\n Locations of peak error:\n" );
for( i=0; i<=n+d+1; i++ )
fprintf( f, "%23.15E\n", mm[i] );
fclose(f);
}


  3 Responses to “Category : C Source Code
Archive   : CEPHES22.ZIP
Filename : REMESW.C
”

  1. Daniel says:
    January 27, 2013 at 3:59 pm

    Very nice! Thank you for this wonderful archive. I wonder why I found it only now. Long live the BBS file archives!

  2. Joshie says:
    March 18, 2014 at 4:57 pm

    This is so awesome! 😀 I’d be cool if you could download an entire archive of this at once, though.

  3. DiskingRound says:
    January 14, 2015 at 10:57 pm

    But one thing that puzzles me is the “mtswslnkmcjklsdlsbdmMICROSOFT” string. There is an article about it here. It is definitely worth a read: http://www.os2museum.com/wp/mtswslnk/