secular-equation.h

Go to the documentation of this file.

/*
 * This file is part of MPSolve 3.2.2
 *
 * Copyright (C) 2001-2020, Dipartimento di Matematica "L. Tonelli", Pisa.
 * License: http://www.gnu.org/licenses/gpl.html GPL version 3 or higher
 *
 * Authors:
 *   Leonardo Robol <leonardo.robol@unipi.it>
 */

 
#ifndef MPS_SECULAR_H_
#define MPS_SECULAR_H_
 
#include <mps/mps.h>
#include <float.h>
 
MPS_BEGIN_DECLS
 
#define MPS_SECULAR_EQUATION(t) (MPS_POLYNOMIAL_CAST (mps_secular_equation, t))
#define MPS_IS_SECULAR_EQUATION(t) (mps_polynomial_check_type (t, "mps_secular_equation"))
 
#ifdef _MPS_PRIVATE
 
/* CONSTANTS */

#define MPS_SECULAR_STARTING_MP_PRECISION 128

#define MPS_SECULAR_EQUIVALENT_FP_PRECISION (MPS_SECULAR_STARTING_MP_PRECISION / 2)
 
struct mps_secular_equation_double_buffer {
  char active;
  mpc_t *ampc1;
  mpc_t *ampc2;
  mpc_t *bmpc1;
  mpc_t *bmpc2;
};

struct mps_secular_equation {
  struct mps_polynomial __base_class__;
 
  struct mps_secular_equation_double_buffer db;

  cplx_t *afpc;

  cdpe_t *adpc;

  cplx_t *bfpc;

  cdpe_t *bdpc;

  mpc_t *ampc;

  pthread_mutex_t * ampc_mutex;

  mpc_t *bmpc;

  pthread_mutex_t * bmpc_mutex;

  double *aafpc;

  double *abfpc;

  rdpe_t *aadpc;

  rdpe_t *abdpc;

  mpc_t *initial_ampc;

  mpc_t *initial_bmpc;

  mpq_t *initial_ampqrc;

  mpq_t *initial_bmpqrc;

  mpq_t *initial_ampqic;

  mpq_t *initial_bmpqic;

  pthread_mutex_t precision_mutex;
};         /* End of struct mps_secular_equation {... */

struct mps_secular_iteration_data {
  long int k;

  mps_boolean radius_set;

  pthread_mutex_t * gs_mutex;

  mpc_t * local_ampc;

  mpc_t * local_bmpc;

  cplx_t * local_afpc;

  cplx_t * local_bfpc;

  cdpe_t * local_adpc;

  cdpe_t * local_bdpc;
};
 
#endif /* #ifdef _MPS_PRIVATE */
 
 
/* MACROS */
#define mps_secular_equation_from_status(s) (mps_secular_equation*)(s)->secular_equation
 
/* Routines in secular-newton.c */
void mps_secular_fnewton (mps_context * st, mps_polynomial * p, mps_approximation * root, cplx_t corr);
void mps_secular_dnewton (mps_context * st, mps_polynomial * p, mps_approximation * root, cdpe_t corr);
void mps_secular_mnewton (mps_context * st, mps_polynomial * p, mps_approximation * root, mpc_t corr, long int wp);
 
/* Routines in secular-regeneartion.c */
mps_boolean mps_secular_ga_regenerate_coefficients (mps_context * s);
 
/* Routines in secular.c */
void mps_secular_deflate (mps_context * s, mps_secular_equation * sec);
 
void mps_secular_check_data (mps_context * s, char *which_case);
 
void mps_secular_restart (mps_context * s);
 
void mps_secular_switch_phase (mps_context * s, mps_phase phase);
 
long int mps_secular_raise_coefficient_precision (mps_context * s, mps_polynomial * p, long int wp);
 
void mps_secular_raise_precision (mps_context * s, int wp);
 
void mps_secular_raise_root_precision (mps_context * s, int wp);
 
/* Routines in secular-starting.c */
void mps_secular_fstart (mps_context * s, mps_secular_equation * sec, mps_approximation ** approximations);
void mps_secular_dstart (mps_context * s, mps_secular_equation * sec, mps_approximation ** approximations);
void mps_secular_mstart (mps_context * s, mps_secular_equation * sec, mps_approximation ** approximations);
 
/* Routines in secular-iteration.c */
int mps_secular_ga_fiterate (mps_context * s, int maxit, mps_boolean just_regenerated);
 
int mps_secular_ga_diterate (mps_context * s, int maxit, mps_boolean just_regenerated);
 
int mps_secular_ga_miterate (mps_context * s, int maxit, mps_boolean just_regenerated);
 
/* Routines in secular-ga.c */
mps_boolean mps_secular_ga_check_stop (mps_context * s);
 
void mps_secular_ga_update_coefficients (mps_context * s);
 
/* Interface functions in secular.c */
mps_secular_equation *mps_secular_equation_new (mps_context * s,
                                                cplx_t * afpc,
                                                cplx_t * bfpc,
                                                unsigned long int n);
 
mps_secular_equation *mps_secular_equation_new_raw (mps_context * s,
                                                    unsigned long int n);
 
void mps_secular_equation_set_coefficient_f (mps_context *ctx, mps_secular_equation *p, int i, mpc_t a, mpc_t b);
void mps_secular_equation_set_coefficient_q (mps_context *ctx, mps_secular_equation *p, int i, mpq_t ar, mpq_t ai, mpq_t br, mpq_t bi);
 
void mps_secular_equation_free (mps_context * ctx, mps_polynomial * p);
 
void mps_secular_set_radii (mps_context * s);
 
mps_boolean mps_secular_poly_feval_with_error (mps_context * ctx, mps_polynomial * p, cplx_t x, cplx_t value, double * error);
 
mps_boolean mps_secular_poly_deval_with_error (mps_context * ctx, mps_polynomial * p, cdpe_t x, cdpe_t value, rdpe_t error);
 
mps_boolean mps_secular_poly_meval_with_error (mps_context * ctx, mps_polynomial * p, mpc_t x, mpc_t value, rdpe_t error);
 
void mps_secular_poly_fstart (mps_context * ctx, mps_polynomial * p, mps_approximation ** approximations);
 
void mps_secular_poly_dstart (mps_context * ctx, mps_polynomial * p, mps_approximation ** approximations);
 
void mps_secular_poly_mstart (mps_context * ctx, mps_polynomial * p, mps_approximation ** approximations);
 
mps_secular_equation * mps_secular_equation_read_from_stream (mps_context * ctx, mps_input_buffer * buffer,
                                                              mps_structure structure, mps_density density,
                                                              long int precision);
 
 
MPS_END_DECLS
 
#endif                          /* SECULAR_H */