hsl_ma77d.h

Go to the documentation of this file.

/* COPYRIGHT (c) 2011, 2013 The Science and Technology Facilities Council (STFC)
 * All Rights Reserved.
 * This code is published under the Eclipse Public License.
 *
 * Authors:  Jonathan Hogg    STFC     2011-05-18
 */
 
#ifndef HSL_MA77D_H
#define HSL_MA77D_H
 
#ifndef ma77_default_control
#define ma77_control ma77_control_d
#define ma77_info ma77_info_d
#define ma77_default_control ma77_default_control_d
#define ma77_open_nelt ma77_open_nelt_d
#define ma77_open ma77_open_d
#define ma77_input_vars ma77_input_vars_d
#define ma77_input_reals ma77_input_reals_d
#define ma77_analyse ma77_analyse_d
#define ma77_factor ma77_factor_d
#define ma77_factor_solve ma77_factor_solve_d
#define ma77_solve ma77_solve_d
#define ma77_resid ma77_resid_d
#define ma77_scale ma77_scale_d
#define ma77_enquire_posdef ma77_enquire_posdef_d
#define ma77_enquire_indef ma77_enquire_indef_d
#define ma77_alter ma77_alter_d
#define ma77_restart ma77_restart_d
#define ma77_finalise ma77_finalise_d
#define ma77_solve_fredholm ma77_solve_fredholm_d
#define ma77_lmultiply ma77_lmultiply_d
#endif
 
typedef double ma77pkgtype_d_;
 
struct ma77_control_d
{
   int f_arrays;
   int print_level;
   int unit_diagnostics;
   int unit_error;
   /* unit for warning messages
    *
    * Printing is suppressed if unit_warning < 0.
    */
   int unit_warning;
   int bits;
   int buffer_lpage[2];
   int buffer_npage[2];
   long int file_size;
   long int maxstore;
   long int storage[3];
   int nemin;
   int maxit;
   int infnorm;
   ma77pkgtype_d_ thresh;
   int nb54;
   int action;
   ma77pkgtype_d_ multiplier;
   int nb64;
   int nbi;
   ma77pkgtype_d_ small;
   ma77pkgtype_d_ static_;
   long int storage_indef;
   ma77pkgtype_d_ u; 
   ma77pkgtype_d_ umin; 
   ma77pkgtype_d_ consist_tol; 
   /* Pad data structure to allow for future growth */
   int ispare[5];
   long int lspare[5];
   ma77pkgtype_d_ rspare[5];
};
 
/***************************************************/
 
struct ma77_info_d
{
   ma77pkgtype_d_ detlog;
   int detsign;
   int flag;
   int iostat;
   int matrix_dup;
   int matrix_rank;
   int matrix_outrange;
   int maxdepth;
   int maxfront;
   long int minstore;
   int ndelay;
   long int nfactor;
   long int nflops;
   int niter;
   int nsup;
   int num_neg;
   int num_nothresh;
   int num_perturbed;
   int ntwo;
   int stat;
   int index[4];
   long int nio_read[2];
   long int nio_write[2];
   long int nwd_read[2];
   long int nwd_write[2];
   int num_file[4];
   long int storage[4];
   int tree_nodes;
   int unit_restart;
   int unused;
   ma77pkgtype_d_ usmall;
 
   /* if we do not have MA77, we assume its is loaded via the linear solver loader, for which we assume HSL 2013 */
   /* Pad data structure to allow for future growth */
   int ispare[5];
   long int lspare[5];
   ma77pkgtype_d_ rspare[5];
};
 
void ma77_default_control_d(
   struct ma77_control_d* control
);
 
void ma77_open_nelt(
   const int                    n,
   const char*                  fname1,
   const char*                  fname2,
   const char*                  fname3,
   const char*                  fname4,
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   const int                    nelt
);
 
void ma77_open_d(
   const int                    n,
   const char*                  fname1,
   const char*                  fname2,
   const char*                  fname3,
   const char*                  fname4,
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_input_vars(
   const int                    idx,
   const int                    nvar,
   const int                    list[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_input_reals_d(
   const int                    idx,
   const int                    length,
   const ma77pkgtype_d_         reals[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_analyse(
   const int                    order[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_factor_d(
   const int                    posdef,
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   const ma77pkgtype_d_*        scale
);
 
void ma77_factor_solve_d(
   const int                    posdef,
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   const ma77pkgtype_d_*        scale,
   const int                    nrhs,
   const int                    lx,
   ma77pkgtype_d_               rhs[]
);
 
void ma77_solve_d(
   const int                    job,
   const int                    nrhs,
   const int                    lx,
   ma77pkgtype_d_               x[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   const ma77pkgtype_d_*        scale
);
 
void ma77_resid_d(
   const int                    nrhs,
   const int                    lx,
   const ma77pkgtype_d_         x[],
   const int                    lresid,
   ma77pkgtype_d_               resid[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   ma77pkgtype_d_*              anorm_bnd
);
 
void ma77_scale_d(
   ma77pkgtype_d_               scale[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   ma77pkgtype_d_*              anorm
);
 
void ma77_enquire_posdef_d(
   ma77pkgtype_d_               d[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_enquire_indef_d(
   int                          piv_order[],
   ma77pkgtype_d_               d[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_alter_d(
   ma77pkgtype_d_               d[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
void ma77_restart_d(
   const char*                  restart_file,
   const char*                  fname1,
   const char*                  fname2,
   const char*                  fname3,
   const char*                  fname4,
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
/* exists only for HSL 2013 */
void ma77_solve_fredholm_d(
   int                          nrhs,
   int                          flag_out[],
   int                          lx,
   ma77pkgtype_d_               x[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   const ma77pkgtype_d_*        scale
);
 
/* exists only for HSL 2013 */
void ma77_lmultiply_d(
   int                          trans,
   int                          k,
   int                          lx,
   ma77pkgtype_d_               x[],
   int                          ly,
   ma77pkgtype_d_               y[],
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info,
   const ma77pkgtype_d_*        scale
);
 
void ma77_finalise_d(
   void**                       keep,
   const struct ma77_control_d* control,
   struct ma77_info_d*          info
);
 
#endif