Procedures

ProcedureLocationProcedure TypeDescription
alps_erroralps_ioSubroutine

Error catching subroutine.

alps_error_initalps_ioSubroutine

Open a file for the error log.

besselIalps_nhdsFunction

Calculates the modified Bessel function of argument x and order n.

BESSIalps_nhdsFunction

Function to calculate the first kind modified Bessel function of integer order N for any real X.

BESSI0generate_distribution.f90Function

This function calculates the modified Bessel function of the first kind of order zero for any positive real argument x.

BESSI0alps_nhdsFunction

Auxiliary Bessel functions for N=0, N=1

BESSI1generate_distribution.f90Function

This function calculates the modified Bessel function of the first kind of order one for any positive real argument x.

BESSI1alps_nhdsFunction

Modified Bessel function of order 1.

BESSJalps_fns_relFunction

This function calculates the first kind Bessel function of integer order N, for any REAL X. We use here the classical recursion formula, when X > N. For X < N, Miller's algorithm is used to avoid overflows. Reference: C.W.Clenshaw, Chebyshev Series for Mathematical Functions, Mathematical Tables, Vol. 5, 1962.

BESSJ0alps_fns_relFunction

This function calculates the first kind Bessel function of order 0, for any REAL X. The polynomial approximation by series of Chebyshev polynomials is used for 0<X<8 and 0<8/X<1. References: M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, 1965. C.W.Clenshaw, Chebyshev Series for Mathematical Functions, Mathematical Tables, Vol. 5, 1962.

BESSJ1alps_fns_relFunction

This subroutine calculates the First Kind Bessel Function of order 1, for any real number X. The polynomial approximation by series of Chebyshev polynomials is used for 0<X<8 and 0<8/X<1. References: M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, 1965. C.W.Clenshaw, Chebyshev Series for Mathematical Functions, Mathematical Tables, Vol. 5, 1962.

BESSKgenerate_distribution.f90Function

This function calculates the modified Bessel function of the third kind of integer order N, for any REAL X. The classical recursion formula is used. Reference: C.W.Clenshaw, Chebyshev Series for Mathematical Functions, Mathematical Tables, Vol. 5, 1962.

BESSK0generate_distribution.f90Function

This function calculates the modified Bessel function of the third kind of order zero for any positive real argument x.

BESSK1generate_distribution.f90Function

This function calculates the modified Bessel function of the third kind of order zero for any positive real argument x.

bM_readalps_ioSubroutine

Reads in bi-Maxwellian parameters.

calc_chialps_nhdsSubroutine

Subroutine that calculates the susceptibility of species j based on NHDS.

calc_eigenalps_fnsSubroutine

This subroutine calculates the relative electric and magnetic field amplitudes, the relative fluctuations in the density and velocity of all species, and the heating rates of the given solution. It is based on the calc_eigen routine by Greg Howes and Kris Klein.

calc_ypsilonalps_nhdsSubroutine

Calculates the Y-tensor according to Stix for a bi-Maxwelling, using the NHDS calculation.

CBESSJalps_fns_relSubroutine

This subroutine calculates the complex Bessel function. It is based on the CBESSJ function release 1.1 by J-P Moreau, Paris (www.jpmoreau.fr).

derivative_f0alps_fnsSubroutine

This subroutine calculates the perpendicular and parallel derivatives of the background velocity distribution function f0.

derivative_f0_relalps_fns_relSubroutine

This subroutine calculates the derivatives of the background velocity distribution function f0 for the relativistic calculation.

determine_bessel_arrayalps_fnsSubroutine

This subroutine determines the array of Bessel functions that is used in the T-tensor of Eq. (2.10) of the code paper.

determine_JTalps_analytSubroutine

This subroutine calculates the transposed Jacobian matrix of the fit function with respect to the fit parameter array.

determine_nmaxalps_fnsSubroutine

This subroutine determines the maximum required order of the Bessel functions in Eq. (2.9) of the code paper.

determine_param_fitalps_analytSubroutine

This is the fitting routine for the hybrid analytic continuation. It determines the full field param_fit.

determine_resonancesalps_fnsSubroutine

This subroutine determines whether any kinetic resonances are located in the integration domain.

determine_sproc_relalps_fns_relSubroutine

This subroutine determines sproc_rel for the given process.

dispalps_fnsFunction

This function returns the determinant of the dispersion tensor for a given frequency om.

dispfunctalps_nhdsFunction

Calculates the dispersion function based on the complex error function.

display_creditsalps_ioSubroutine

Writes the opening credits.

distribution_analytalps_distribution_analytFunction

This function returns the pre-defined function as f0.

eval_fitalps_analytFunction

This function evaluates the fit to f0 at and the complex parallel momentum ppar_valC. It requires the fit parameters that will be determined by the subroutine determine_param_fit.

Factalps_fns_relFunction

This function returns the factorial k! of its argument k.

find_minimaalps_fnsSubroutine

This subroutine identifies the minima of the coarse map grid. It is called by map_search. The code is based on a routine by Greg Howes, 2006.

fit_functionalps_analytFunction

This function evaluates the fit to f0 at real pperp_val and complex ppar_val, provided that the one-dimensional fit-parameter array params is fed into the function. This is only used during the fitting. For the evaluation in ALPS, use eval_fit.

fit_readalps_ioSubroutine

Reads in fit parameters for component is.

full_integratealps_fnsFunction

This function returns the full integral expression according to Eq. (2.9) in the code paper.

funct_galps_fnsFunction

This function returns the function from Eq. (3.2) of the code paper.

funct_g_relalps_fns_relFunction

This function returns the function from Eq. (3.2) of the code paper for a relativistic calculation.

Gammaalps_fns_relFunction

This function returns the Gamma-function.

get_indexed_double_namelist_unitalps_ioSubroutine

A version of get_indexed_namelist_unit, extended to allow for double indexing in order to read in multiple fits for a single species.

get_indexed_namelist_unitalps_ioSubroutine

Determines unused I/O unit.

get_indexed_namelist_unitgenerate_distributionSubroutine

Determines unused I/O unit.

get_runnamealps_ioSubroutine

Get runname for output files from input argument.

get_runnameinterpolateSubroutine

Get runname for output files from input argument.

get_runnamegenerate_distributionSubroutine

Get runname for output files from input argument.

get_unused_unitalps_ioSubroutine

Determine unused number for I/O index.

get_unused_unitinterpolateSubroutine

This subroutine returns an available unit number. It is taken from the AstroGK code.

get_unused_unitgenerate_distributionSubroutine

Determine unused number for I/O index.

init_paramalps_ioSubroutine

Read in system parameters from *.in file. Only processor 0 calls this routine:

input_unitalps_ioFunction

Assigns input unit for namelist opening.

input_unitinterpolateFunction

This function returns a unit number for a namelist. It is taken from the AstroGK code.

input_unitgenerate_distributionFunction

Assigns input unit for namelist opening.

input_unit_existalps_ioFunction

Determine if a particular namelist already opened.

input_unit_existinterpolateFunction

This function checks whether a unit number exists. It is taken from the AstroGK code.

input_unit_existgenerate_distributionFunction

Determine if a particular namelist already opened.

int_eealps_fnsFunction

This function returns the ee term in Eq. (2.9).

int_ee_relalps_fns_relFunction

This function returns the ee term in Eq. (2.9) for the relativistic calculation.

int_Talps_fnsFunction

This function returns the T-tensor according to Eq. (2.10) of the code paper.

int_T_relalps_fns_relFunction

This function returns the T-tensor according to Eq. (2.10) of the code paper for the relativistic calculation.

int_T_resalps_fnsFunction

This function returns the T-tensor according to Eq. (2.10) of the code paper for the case in which it is evaluated at the complex resonance momentum.

int_T_res_relalps_fns_relFunction

This function returns the T-tensor according to Eq. (2.10) of the code paper for the case in which it is evaluated at the complex resonance momentum for the relativistic calculation.

integratealps_fnsFunction

This function performs the integral in Eq. (2.9) of the code paper, but without accounting for the Landau contour integral. It is called by full_integrate.

integrate_resalps_fnsFunction

This function performs the integration near resonances as described in Section 3.1 of the code paper. It is only called if resonances are present in or near the integration domain.

integrate_res_relalps_fns_relFunction

This function performs the integration near resonances as described in Section 3.1 of the code paper for a relativistic calculation. It is only called if resonances are present in or near the integration domain.

integrate_resU_relalps_fns_relFunction

This function evaluates the integral with the integrand proportional to in Eq. (2.9) of the code paper for a relativistic calculation.

isnancheckalps_ioFunction

Checks if double precision number input is NaN.

landau_integratealps_fnsFunction

This function evaluates the Landau contour according to Eqs. (3.8) and (3.9) of the code paper.

landau_integrate_relalps_fns_relFunction

This function evaluates the Landau contour according to Eqs. (3.8) and (3.9) of the code paper for a relativistic calculation.

LM_nonlinear_fitalps_analytSubroutine

This subroutine processes the nonlinear Levenberg-Marquart algorithm and returns the one-dimensional array params at a given iperp. The variable quality is the sum of the squares of all residuals.

map_readalps_ioSubroutine

Reads in complex frequency map parameters.

map_searchalps_fnsSubroutine

This subroutine calculates the map of the determinant of the dispersion tensor in complex frequency space.

normalizeinterpolation.f90Subroutine

This subroutine normalises the fine interpolation grid.

om_double_scanalps_fnsSubroutine

This subroutine scans along a prescribed plane in wavevector space to map out in this space. It is required that n_scan=2.

om_scanalps_fnsSubroutine

This subroutine scans solutions along a single prescribed path in wavevector space.

output_fitalps_analytSubroutine

This subroutine outputs the fit parameters for iperp=0 to stdout to monitor the fit.

output_timealps_ioSubroutine

Outputs the date and time in a given format using intrinsic FORTRAN function.

pass_distributionalps_comSubroutine

Passes distribution functions and associated parameters.

pass_instructionsalps_comSubroutine

Passes information between processes.

polyharmonic_splineinterpolation.f90Subroutine

This soubroutine interpolates the grid with a polyharmonic thin-plate spline. This subroutine needs the LUPACK and BLAS libraries to evoke the dgesv subroutine. The method uses the Thin Plate Spline. We use these resources: http://cseweb.ucsd.edu/~sjb/eccv_tps.pdf http://www.univie.ac.at/nuhag-php/bibtex/open_files/po94_M%20J%20D%20Powell%2003%2093.pdf http://vision.ucsd.edu/sites/default/files/fulltext(4).pdf

polyharmonic_splinealps_fns_relSubroutine

This soubroutine interpolates the grid with a polyharmonic thin-plate spline. This subroutine needs the LUPACK and BLAS libraries to evoke the dgesv subroutine. The method uses the Thin Plate Spline. We use these resources: http://cseweb.ucsd.edu/~sjb/eccv_tps.pdf http://www.univie.ac.at/nuhag-php/bibtex/open_files/po94_M%20J%20D%20Powell%2003%2093.pdf http://vision.ucsd.edu/sites/default/files/fulltext(4).pdf

principal_integral_relalps_fns_relFunction

This function performs the integration near resonances as described in Section 3.1 of the code paper for a relativistic calculation. It is only called if resonances are present in or near the integration domain.

read_f0alps_ioSubroutine

Subroutine for reading in background distribution function

read_in_paramsinterpolateSubroutine

This subroutine reads in system parameters input file (namelist) as argument after executable: ./interpolate input.in

read_in_paramsgenerate_distributionSubroutine

This subroutine reads in the parameters for the generation of distribution functions for ALPS.

refine_guessalps_fnsSubroutine

This subroutine refines the guess at the starting point of the search for solutions to the dispersion relation when scanning. It is also used by map_search to identify the roots on the map.

resUalps_fnsFunction

This function evaluates the term proportional to in Eq. (2.9) of the code paper.

resU_relalps_fns_relFunction

This function evaluates the term proportional to in Eq. (2.9) of the code paper for the relativistic calculation.

scan_readalps_ioSubroutine

The most important subroutine. Reads in wavevector scan parameters. Defines scanner, which controls the behavior of the wavevector scan.

secantalps_fnsSubroutine

This subroutine applies the secant method to find the roots of the dispersion tensor.

solution_readalps_ioSubroutine

Reads in initial guesses for dispersion solutions.

spec_readalps_ioSubroutine

Subroutine for reading in species parameters

spec_readgenerate_distributionSubroutine

Subroutine for reading in species parameters.

split_processesalps_fnsSubroutine

This subroutine defines the tasks for the individual processes. It uses the number of species and the required orders of the Bessel functions to define the splitting across the MPI processes.

WOFZalps_nhdsSubroutine

Given a complex number Z = (XI,YI), this subroutine computes the value of the Faddeeva function W(Z) = exp(-Z2)erfc(-IZ), where erfc is the complex complementary error function and I is the imaginary unit.

call~~graph~~CallGraph proc~alps_error alps_io::alps_error proc~output_time alps_io::output_time proc~alps_error->proc~output_time mpi_abort mpi_abort proc~alps_error->mpi_abort proc~alps_error_init alps_io::alps_error_init proc~get_runname alps_io::get_runname proc~alps_error_init->proc~get_runname proc~get_unused_unit alps_io::get_unused_unit proc~alps_error_init->proc~get_unused_unit proc~besseli alps_nhds::besselI proc~bessi alps_nhds::BESSI proc~besseli->proc~bessi proc~bessi0 alps_nhds::BESSI0 proc~bessi->proc~bessi0 proc~bessi1 alps_nhds::BESSI1 proc~bessi->proc~bessi1 dfloat dfloat proc~bessi->dfloat float float proc~bessi->float proc~bessi0~2 generate_distribution.f90::BESSI0 proc~bessi1~2 generate_distribution.f90::BESSI1 proc~bessj alps_fns_rel::BESSJ proc~bessj0 alps_fns_rel::BESSJ0 proc~bessj->proc~bessj0 proc~bessj1 alps_fns_rel::BESSJ1 proc~bessj->proc~bessj1 proc~bessj->float proc~bessk generate_distribution.f90::BESSK proc~bessk->dfloat proc~bessk0 generate_distribution.f90::BESSK0 proc~bessk1 generate_distribution.f90::BESSK1 proc~bm_read alps_io::bM_read proc~calc_chi alps_nhds::calc_chi proc~calc_chi->proc~besseli proc~calc_ypsilon alps_nhds::calc_ypsilon proc~calc_chi->proc~calc_ypsilon proc~calc_eigen alps_fns::calc_eigen proc~disp alps_fns::disp proc~calc_eigen->proc~disp proc~calc_ypsilon->proc~besseli proc~dispfunct alps_nhds::dispfunct proc~calc_ypsilon->proc~dispfunct proc~cbessj alps_fns_rel::CBESSJ proc~fact alps_fns_rel::Fact proc~cbessj->proc~fact proc~derivative_f0 alps_fns::derivative_f0 proc~derivative_f0_rel alps_fns_rel::derivative_f0_rel proc~derivative_f0->proc~derivative_f0_rel proc~derivative_f0->proc~get_unused_unit proc~derivative_f0_rel->proc~get_unused_unit proc~polyharmonic_spline alps_fns_rel::polyharmonic_spline proc~derivative_f0_rel->proc~polyharmonic_spline proc~determine_bessel_array alps_fns::determine_bessel_array proc~determine_bessel_array->proc~bessj proc~determine_jt alps_analyt::determine_JT proc~determine_nmax alps_fns::determine_nmax proc~determine_nmax->proc~bessj mpi_barrier mpi_barrier proc~determine_nmax->mpi_barrier proc~determine_param_fit alps_analyt::determine_param_fit proc~lm_nonlinear_fit alps_analyt::LM_nonlinear_fit proc~determine_param_fit->proc~lm_nonlinear_fit proc~output_fit alps_analyt::output_fit proc~determine_param_fit->proc~output_fit proc~determine_resonances alps_fns::determine_resonances proc~determine_sproc_rel alps_fns_rel::determine_sproc_rel proc~disp->proc~calc_chi proc~disp->proc~determine_resonances proc~full_integrate alps_fns::full_integrate proc~disp->proc~full_integrate proc~int_ee alps_fns::int_ee proc~disp->proc~int_ee proc~int_ee_rel alps_fns_rel::int_ee_rel proc~disp->proc~int_ee_rel proc~disp->mpi_barrier mpi_bcast mpi_bcast proc~disp->mpi_bcast mpi_reduce mpi_reduce proc~disp->mpi_reduce proc~wofz alps_nhds::WOFZ proc~dispfunct->proc~wofz proc~display_credits alps_io::display_credits proc~distribution_analyt alps_distribution_analyt::distribution_analyt proc~eval_fit alps_analyt::eval_fit proc~eval_fit->proc~distribution_analyt proc~fit_function alps_analyt::fit_function proc~eval_fit->proc~fit_function proc~find_minima alps_fns::find_minima proc~fit_read alps_io::fit_read proc~integrate alps_fns::integrate proc~full_integrate->proc~integrate proc~integrate_res alps_fns::integrate_res proc~full_integrate->proc~integrate_res proc~integrate_res_rel alps_fns_rel::integrate_res_rel proc~full_integrate->proc~integrate_res_rel proc~landau_integrate alps_fns::landau_integrate proc~full_integrate->proc~landau_integrate proc~landau_integrate_rel alps_fns_rel::landau_integrate_rel proc~full_integrate->proc~landau_integrate_rel proc~funct_g alps_fns::funct_g proc~int_t alps_fns::int_T proc~funct_g->proc~int_t proc~funct_g_rel alps_fns_rel::funct_g_rel proc~int_t_rel alps_fns_rel::int_T_rel proc~funct_g_rel->proc~int_t_rel proc~gamma alps_fns_rel::Gamma proc~get_indexed_double_namelist_unit alps_io::get_indexed_double_namelist_unit proc~get_indexed_double_namelist_unit->proc~alps_error proc~get_indexed_double_namelist_unit->proc~get_unused_unit proc~input_unit alps_io::input_unit proc~get_indexed_double_namelist_unit->proc~input_unit proc~input_unit_exist alps_io::input_unit_exist proc~get_indexed_double_namelist_unit->proc~input_unit_exist proc~get_indexed_namelist_unit alps_io::get_indexed_namelist_unit proc~get_indexed_namelist_unit->proc~alps_error proc~get_indexed_namelist_unit->proc~get_unused_unit proc~get_indexed_namelist_unit->proc~input_unit proc~get_indexed_namelist_unit->proc~input_unit_exist proc~get_indexed_namelist_unit~2 generate_distribution::get_indexed_namelist_unit proc~get_unused_unit~3 generate_distribution::get_unused_unit proc~get_indexed_namelist_unit~2->proc~get_unused_unit~3 proc~input_unit_exist~3 generate_distribution::input_unit_exist proc~get_indexed_namelist_unit~2->proc~input_unit_exist~3 proc~input_unit~3 generate_distribution::input_unit proc~get_indexed_namelist_unit~2->proc~input_unit~3 proc~get_runname~2 interpolate::get_runname proc~get_runname~3 generate_distribution::get_runname proc~get_unused_unit~2 interpolate::get_unused_unit proc~init_param alps_io::init_param proc~init_param->proc~bm_read proc~init_param->proc~fit_read proc~init_param->proc~get_indexed_double_namelist_unit proc~init_param->proc~get_indexed_namelist_unit proc~init_param->proc~get_unused_unit proc~map_read alps_io::map_read proc~init_param->proc~map_read proc~scan_read alps_io::scan_read proc~init_param->proc~scan_read proc~solution_read alps_io::solution_read proc~init_param->proc~solution_read proc~spec_read alps_io::spec_read proc~init_param->proc~spec_read proc~input_unit_exist~2 interpolate::input_unit_exist proc~input_unit~2 interpolate::input_unit proc~int_ee_rel->proc~determine_sproc_rel proc~int_t_rel->proc~bessj proc~int_t_res alps_fns::int_T_res proc~int_t_res_rel alps_fns_rel::int_T_res_rel proc~int_t_res_rel->proc~cbessj proc~integrate->proc~int_t proc~resu alps_fns::resU proc~integrate->proc~resu proc~integrate_res->proc~funct_g proc~integrate_res->proc~integrate proc~integrate_res_rel->proc~determine_sproc_rel proc~integrate_resu_rel alps_fns_rel::integrate_resU_rel proc~integrate_res_rel->proc~integrate_resu_rel proc~integrate_resu_rel->proc~alps_error proc~integrate_resu_rel->proc~int_t_rel proc~principal_integral_rel alps_fns_rel::principal_integral_rel proc~integrate_resu_rel->proc~principal_integral_rel proc~resu_rel alps_fns_rel::resU_rel proc~integrate_resu_rel->proc~resu_rel proc~isnancheck alps_io::isnancheck proc~landau_integrate->proc~eval_fit proc~landau_integrate->proc~int_t_res proc~landau_integrate_rel->proc~determine_sproc_rel proc~landau_integrate_rel->proc~eval_fit proc~landau_integrate_rel->proc~int_t_res_rel proc~lm_nonlinear_fit->proc~alps_error proc~lm_nonlinear_fit->proc~determine_jt proc~lm_nonlinear_fit->proc~fit_function dgetrf dgetrf proc~lm_nonlinear_fit->dgetrf dgetri dgetri proc~lm_nonlinear_fit->dgetri proc~map_search alps_fns::map_search proc~map_search->proc~disp proc~map_search->proc~find_minima proc~map_search->proc~get_unused_unit proc~refine_guess alps_fns::refine_guess proc~map_search->proc~refine_guess proc~map_search->mpi_bcast proc~normalize interpolation.f90::normalize proc~om_double_scan alps_fns::om_double_scan proc~om_double_scan->proc~alps_error proc~om_double_scan->proc~calc_eigen proc~om_double_scan->proc~determine_bessel_array proc~om_double_scan->proc~determine_nmax proc~om_double_scan->proc~disp proc~om_double_scan->proc~get_unused_unit proc~om_double_scan->proc~isnancheck proc~secant alps_fns::secant proc~om_double_scan->proc~secant proc~split_processes alps_fns::split_processes proc~om_double_scan->proc~split_processes proc~om_double_scan->mpi_barrier proc~om_double_scan->mpi_bcast proc~om_scan alps_fns::om_scan proc~om_scan->proc~alps_error proc~om_scan->proc~calc_eigen proc~om_scan->proc~determine_bessel_array proc~om_scan->proc~determine_nmax proc~om_scan->proc~disp proc~om_scan->proc~get_unused_unit proc~om_scan->proc~isnancheck proc~om_scan->proc~secant proc~om_scan->proc~split_processes proc~om_scan->mpi_barrier proc~om_scan->mpi_bcast proc~output_fit->proc~alps_error proc~output_fit->proc~eval_fit proc~output_fit->proc~isnancheck proc~pass_distribution alps_com::pass_distribution proc~pass_distribution->mpi_bcast proc~pass_instructions alps_com::pass_instructions proc~pass_instructions->mpi_bcast dgesv dgesv proc~polyharmonic_spline->dgesv proc~polyharmonic_spline~2 interpolation.f90::polyharmonic_spline proc~polyharmonic_spline~2->dgesv proc~principal_integral_rel->proc~funct_g_rel proc~read_f0 alps_io::read_f0 proc~read_f0->proc~get_unused_unit proc~read_in_params interpolate::read_in_params proc~read_in_params->proc~get_runname~2 proc~read_in_params->proc~get_unused_unit~2 proc~read_in_params~2 generate_distribution::read_in_params proc~read_in_params~2->proc~get_indexed_namelist_unit~2 proc~read_in_params~2->proc~get_runname~3 proc~read_in_params~2->proc~get_unused_unit~3 proc~spec_read~2 generate_distribution::spec_read proc~read_in_params~2->proc~spec_read~2 proc~refine_guess->proc~disp proc~refine_guess->proc~get_unused_unit proc~refine_guess->proc~secant proc~refine_guess->mpi_barrier proc~secant->proc~disp proc~secant->mpi_barrier dcos dcos proc~wofz->dcos dexp dexp proc~wofz->dexp dsin dsin proc~wofz->dsin idint idint proc~wofz->idint idnint idnint proc~wofz->idnint program~alps alps program~alps->proc~alps_error program~alps->proc~alps_error_init program~alps->proc~derivative_f0 program~alps->proc~determine_bessel_array program~alps->proc~determine_nmax program~alps->proc~determine_param_fit program~alps->proc~display_credits program~alps->proc~init_param program~alps->proc~map_search program~alps->proc~om_double_scan program~alps->proc~om_scan program~alps->proc~output_time program~alps->proc~pass_distribution program~alps->proc~pass_instructions program~alps->proc~read_f0 program~alps->proc~refine_guess program~alps->proc~split_processes program~generate_distribution generate_distribution program~generate_distribution->proc~distribution_analyt program~generate_distribution->proc~read_in_params~2 program~interpolate interpolate program~interpolate->proc~normalize program~interpolate->proc~polyharmonic_spline~2 program~interpolate->proc~read_in_params
Help