fnft_kdvv.h

Go to the documentation of this file.

/*
* This file is part of FNFT.
*
* FNFT is free software; you can redistribute it and/or
* modify it under the terms of the version 2 of the GNU General
* Public License as published by the Free Software Foundation.
*
* FNFT is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*
* Contributors:
* Sander Wahls (TU Delft) 2017-2018, 2023.
* Shrinivas Chimmalgi (TU Delft) 2019-2020.
* Peter J Prins (TU Delft) 2020-2021.
*/
 
#ifndef FNFT_KDVV_H
#define FNFT_KDVV_H
 
#include "fnft__kdv_fscatter.h"
#include "fnft__kdv_scatter.h"
#include "fnft__misc.h" // for l2norm
#include <string.h> // for memcpy
#include <stdio.h>
#include "fnft__errwarn.h"
#include "fnft__poly_roots_fasteigen.h"
#include "fnft__poly_chirpz.h"
 
typedef enum {
    fnft_kdvv_bsloc_NEWTON,
    fnft_kdvv_bsloc_GRIDSEARCH_AND_REFINE
} fnft_kdvv_bsloc_t;
 
typedef enum {
    fnft_kdvv_dstype_NORMING_CONSTANTS,
    fnft_kdvv_dstype_RESIDUES,
    fnft_kdvv_dstype_BOTH
} fnft_kdvv_dstype_t;
 
typedef enum {
    fnft_kdvv_cstype_REFLECTION_COEFFICIENT,
    fnft_kdvv_cstype_AB,
    fnft_kdvv_cstype_BOTH
} fnft_kdvv_cstype_t;
 
typedef struct {
    fnft_kdvv_bsloc_t bound_state_localization;
    FNFT_UINT niter;
    fnft_kdvv_dstype_t discspec_type;
    fnft_kdvv_cstype_t contspec_type;
    FNFT_INT normalization_flag;
    fnft_kdv_discretization_t discretization;
    FNFT_UINT richardson_extrapolation_flag;
    FNFT_REAL grid_spacing;
} fnft_kdvv_opts_t;
 
fnft_kdvv_opts_t fnft_kdvv_default_opts();
 
FNFT_INT fnft_kdvv(const FNFT_UINT D, FNFT_COMPLEX const * const q,
    FNFT_REAL const * const T, const FNFT_UINT M,
    FNFT_COMPLEX * const contspec, FNFT_REAL const * const XI,
    FNFT_UINT * const K_ptr, FNFT_COMPLEX * const bound_states,
    FNFT_COMPLEX * const normconsts_or_residues, //const FNFT_INT kappa,
    fnft_kdvv_opts_t *opts);
 
 
#ifdef FNFT_ENABLE_SHORT_NAMES
#define kdvv_bsloc_NEWTON fnft_kdvv_bsloc_NEWTON
#define kdvv_bsloc_GRIDSEARCH_AND_REFINE fnft_kdvv_bsloc_GRIDSEARCH_AND_REFINE
#define kdvv_dstype_NORMING_CONSTANTS fnft_kdvv_dstype_NORMING_CONSTANTS
#define kdvv_dstype_RESIDUES fnft_kdvv_dstype_RESIDUES
#define kdvv_dstype_BOTH fnft_kdvv_dstype_BOTH
#define kdvv_cstype_REFLECTION_COEFFICIENT fnft_kdvv_cstype_REFLECTION_COEFFICIENT
#define kdvv_cstype_AB fnft_kdvv_cstype_AB
#define kdvv_cstype_BOTH fnft_kdvv_cstype_BOTH
#endif
 
#endif