beyn.h

Go to the documentation of this file.

 
#ifndef BEYN_H
#define BEYN_H
 
#include <stddef.h>
#include <complex.h>
 
 
typedef int (*beyn_function_M_t)(complex double *target_M, size_t m, complex double z, void *params);
 
 
typedef int (*beyn_function_M_inv_Vhat_t)(complex double *target_M_inv_Vhat, size_t m, size_t l,
            const complex double *Vhat, complex double z, void *params); 
 
typedef struct beyn_contour_t {
    size_t n; 
 
    complex double centre;
    _Bool (*inside_test)(struct beyn_contour_t *, complex double z);
    complex double z_dz[][2]; 
} beyn_contour_t;
 
 
beyn_contour_t *beyn_contour_ellipse(complex double centre, double halfax_re, double halfax_im, size_t npoints);
 
 
typedef enum {
    BEYN_CONTOUR_HALFELLIPSE_RE_PLUS = 3,
    BEYN_CONTOUR_HALFELLIPSE_RE_MINUS = 1,
    BEYN_CONTOUR_HALFELLIPSE_IM_PLUS = 0,
    BEYN_CONTOUR_HALFELLIPSE_IM_MINUS = 2,
} beyn_contour_halfellipse_orientation;
 
 
 
beyn_contour_t *beyn_contour_halfellipse(complex double centre, double halfax_re, double halfax_im, size_t npoints,
        beyn_contour_halfellipse_orientation or);
 
beyn_contour_t *beyn_contour_kidney(complex double centre, double halfax_re, double halfax_im, 
    double rounding, 
    size_t n, beyn_contour_halfellipse_orientation or);
 
 
typedef struct beyn_result_t {
    size_t neig; 
    size_t vlen; 
    complex double *eigval;
    complex double *eigval_err;
    double *residuals;
    complex double *eigvec; // Rows are the eigenvectors
    double *ranktest_SV;
 
} beyn_result_t;
 
void beyn_result_free(beyn_result_t *result);
 
beyn_result_t *beyn_solve(
        size_t m, 
        size_t l, 
        beyn_function_M_t M, 
        beyn_function_M_inv_Vhat_t M_inv_Vhat, 
        void *params, 
        const beyn_contour_t *contour, 
        double rank_tol, 
        size_t rank_min_sel, 
        double res_tol 
          );
 
#endif // BEYN_H