secant_osc Subroutine

public subroutine secant_osc(om, in)

Uses

  • proc~~secant_osc~~UsesGraph proc~secant_osc secant_osc module~alps_var alps_var proc~secant_osc->module~alps_var mpi mpi proc~secant_osc->mpi module~alps_var->mpi

Secant method with adaptive damping and Newton search fallback. Ensure we don’t divide by a tiny value Check convergence Detect oscillations over last four iterations Reduce step size If oscillation persists, use finite-difference Newton step Regularized Newton step Apply jump limits to prevent large jumps Cap jump at 10% of om If |D| increases, reduce jump Update previous values Apply update

Arguments

Type IntentOptional Attributes Name
double complex, intent(inout) :: om

Complex wave frequency .
integer, intent(in) :: in

Root number

Calls

proc~~secant_osc~~CallsGraph proc~secant_osc secant_osc mpi_barrier mpi_barrier proc~secant_osc->mpi_barrier proc~disp disp proc~secant_osc->proc~disp proc~disp->mpi_barrier mpi_bcast mpi_bcast proc~disp->mpi_bcast mpi_reduce mpi_reduce proc~disp->mpi_reduce proc~calc_chi calc_chi proc~disp->proc~calc_chi proc~determine_resonances determine_resonances proc~disp->proc~determine_resonances proc~full_integrate full_integrate proc~disp->proc~full_integrate proc~int_ee int_ee proc~disp->proc~int_ee proc~int_ee_rel int_ee_rel proc~disp->proc~int_ee_rel proc~besseli besselI proc~calc_chi->proc~besseli proc~calc_chi_cold calc_chi_cold proc~calc_chi->proc~calc_chi_cold proc~calc_ypsilon calc_ypsilon proc~calc_chi->proc~calc_ypsilon proc~integrate integrate proc~full_integrate->proc~integrate proc~integrate_res integrate_res proc~full_integrate->proc~integrate_res proc~integrate_res_rel integrate_res_rel proc~full_integrate->proc~integrate_res_rel proc~landau_integrate landau_integrate proc~full_integrate->proc~landau_integrate proc~landau_integrate_rel landau_integrate_rel proc~full_integrate->proc~landau_integrate_rel proc~determine_sproc_rel determine_sproc_rel proc~int_ee_rel->proc~determine_sproc_rel proc~bessi BESSI proc~besseli->proc~bessi proc~calc_ypsilon->proc~besseli proc~dispfunct dispfunct proc~calc_ypsilon->proc~dispfunct proc~int_t int_T proc~integrate->proc~int_t proc~resu resU proc~integrate->proc~resu proc~integrate_res->proc~integrate proc~funct_g funct_g proc~integrate_res->proc~funct_g proc~integrate_res_rel->proc~determine_sproc_rel proc~integrate_resu_rel integrate_resU_rel proc~integrate_res_rel->proc~integrate_resu_rel proc~eval_fit eval_fit proc~landau_integrate->proc~eval_fit proc~int_t_res int_T_res proc~landau_integrate->proc~int_t_res proc~landau_integrate_rel->proc~determine_sproc_rel proc~landau_integrate_rel->proc~eval_fit proc~int_t_res_rel int_T_res_rel proc~landau_integrate_rel->proc~int_t_res_rel dfloat dfloat proc~bessi->dfloat float float proc~bessi->float proc~bessi0~2 BESSI0 proc~bessi->proc~bessi0~2 proc~bessi1~2 BESSI1 proc~bessi->proc~bessi1~2 proc~wofz WOFZ proc~dispfunct->proc~wofz proc~distribution_analyt distribution_analyt proc~eval_fit->proc~distribution_analyt proc~fit_function fit_function proc~eval_fit->proc~fit_function proc~fit_function_poly fit_function_poly proc~eval_fit->proc~fit_function_poly proc~funct_g->proc~int_t proc~cbessj CBESSJ proc~int_t_res_rel->proc~cbessj proc~alps_error alps_error proc~integrate_resu_rel->proc~alps_error proc~int_t_rel int_T_rel proc~integrate_resu_rel->proc~int_t_rel proc~principal_integral_rel principal_integral_rel proc~integrate_resu_rel->proc~principal_integral_rel proc~resu_rel resU_rel proc~integrate_resu_rel->proc~resu_rel mpi_abort mpi_abort proc~alps_error->mpi_abort proc~output_time output_time proc~alps_error->proc~output_time proc~fact Fact proc~cbessj->proc~fact proc~bessj BESSJ proc~int_t_rel->proc~bessj proc~funct_g_rel funct_g_rel proc~principal_integral_rel->proc~funct_g_rel 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 proc~bessj->float proc~bessj0 BESSJ0 proc~bessj->proc~bessj0 proc~bessj1 BESSJ1 proc~bessj->proc~bessj1 proc~funct_g_rel->proc~int_t_rel

Called by

proc~~secant_osc~~CalledByGraph proc~secant_osc secant_osc proc~om_double_scan om_double_scan proc~om_double_scan->proc~secant_osc proc~om_scan om_scan proc~om_scan->proc~secant_osc proc~refine_guess refine_guess proc~refine_guess->proc~secant_osc proc~map_search map_search proc~map_search->proc~refine_guess program~alps alps program~alps->proc~om_double_scan program~alps->proc~om_scan program~alps->proc~refine_guess program~alps->proc~map_search