/* remesw.c */

/* Write answers to disc file. */

#include <stdio.h>
#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<n+d+1; 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);
}