Sets the output file or the error message output file.

Synopsis with Optional Arguments

#include <imsl.h>

void imsl_output_file (
IMSL_SET_OUTPUT_FILE, FILE *ofile,
IMSL_GET_OUTPUT_FILE, FILE **pofile,
IMSL_SET_ERROR_FILE, FILE *efile,
IMSL_GET_ERROR_FILE, FILE **pefile,
0)

Optional Arguments

IMSL_SET_OUTPUT_FILE, FILE *ofile   (Input)
Set the output file to ofile.
Default: ofile = stdout

IMSL_GET_OUTPUT_FILE, FILE **pfile   (Output)
Set the FILE pointed to by pfile to the current output file.

IMSL_SET_ERROR_FILE, FILE *efile   (Input)
Set the error message output file to efile.
Default: efile = stderr

IMSL_GET_ERROR_FILE, FILE **pefile   (Output)
Set the FILE pointed to by pefile to the error message output file.

Description

This function allows the file used for printing by IMSL routines to be changed.

If multiple threads are used then default settings are valid for each thread. When using threads it is possible to set different output files for each thread by calling imsl_output_file in
Chapter 15 of the IMSL Stat Numerical Libraries from within each thread.  See Example 2 for details.

Examples

Example 1

This example opens the file myfile and changes the output file to this new file.
The function imsl_f_write_matrix then writes to this file.

#include <stdio.h>
#include <imsl.h>

main()
{
    FILE        *ofile;
    float       x[] = {3.0, 2.0, 1.0};

    imsl_f_write_matrix ("x (default file)", 1, 3, x, 0);

    ofile = fopen("myfile", "w");
    imsl_output_file(IMSL_SET_OUTPUT_FILE, ofile,
                     0);
    imsl_f_write_matrix ("x (myfile)", 1, 3, x, 0);
}

Output

x (default file)
1           2           3
3           2           1

File myfile

x (myfile)
1           2           3
3           2           1

 

Example 2

The following example illustrates how to direct output from IMSL routines that run in separate threads to different files.  First, two threads are created, each calling a different IMSL function, then the results are printed by calling imsl_f_write_matrix from within each thread. Note that imsl_output_file is called from within each thread to change the default output file.

 

#include <pthread.h>

#include <stdio.h>

#include "imsl.h"

 

void *ex1(void* arg);

void *ex2(void* arg);

void main()

{

  pthread_t       thread1;

  pthread_t       thread2;

 

  /* Disable IMSL signal trapping. */

  imsl_error_options(IMSL_SET_SIGNAL_TRAPPING, 0, 0);

 

  /* Create two threads. */

  if (pthread_create(&thread1, NULL ,ex1, (void *)NULL) != 0)

    perror("pthread_create"), exit(1);

  if (pthread_create(&thread2, NULL ,ex2, (void *)NULL) != 0)

    perror("pthread_create"), exit(1);

 

  /* Wait for threads to finish. */

  if (pthread_join(thread1, NULL) != 0)

    perror("pthread_join"),exit(1);

  if (pthread_join(thread2, NULL) != 0)

    perror("pthread_join"),exit(1);

 

}

 

void *ex1(void* arg)

{

  float *rand_nums = NULL;

  FILE  *file_ptr;

  /* Open a file to write the result in. */

  file_ptr = fopen("ex1.out", "w");

  /* Set the output file for this thread. */

  imsl_output_file(IMSL_SET_OUTPUT_FILE, file_ptr, 0);

  /* Compute 5 random numbers. */

  imsl_random_seed_set(12345);

  rand_nums = imsl_f_random_uniform(5, 0);

  /* Output random numbers. */

  imsl_f_write_matrix("Random Numbers", 5, 1, rand_nums, 0);

 

  if (rand_nums) free(rand_nums);

  fclose(file_ptr);

}

void *ex2(void* arg)

{     

  int n =     3;

  float *x;  

  float a[] = {1.0, 3.0, 3.0,

              1.0, 3.0, 4.0,     

              1.0, 4.0, 3.0};    

  float b[] = {1.0, 4.0, -1.0};

  FILE  *file_ptr;

  /* Open a file to write the result in. */

  file_ptr = fopen("ex2.out", "w");

  /* Set the output file for this thread. */

  imsl_output_file(IMSL_SET_OUTPUT_FILE, file_ptr, 0);

  /* Solve Ax = b for x */

  x = imsl_f_lin_sol_gen (n, a, b, 0);

  /* Print x */

  imsl_f_write_matrix ("Solution, x, of Ax = b", 1, 3, x, 0);

 

  if (x) free(x);

  fclose(file_ptr);

}

 

Output

 

ex1.out

 

Random Numbers

 1      0.0966

 2      0.8340

 3      0.9477

 4      0.0359

 5      0.0115

 

 

ex2.out

 

      Solution, x, of Ax = b

         1           2           3

        -2          -2           3

 

Visual Numerics, Inc. Visual Numerics - Developers of IMSL and PV-WAVE http://www.vni.com/ PHONE: 713.784.3131 FAX:713.781.9260