pod.f

pod.f#

These are the subroutines in pod.f.

Subroutines and functions

subroutine  setbases()#

set ub, vb, wb, tb in /morbasis/ with pod basis

Called from:

cres(), c_avei(), set_sigma(), set_nres(), read_sigma(), c_sigma(), set_theta_poisson(), set_theta_heat(), set_theta_ad(), set_theta_ns(), rom_poisson(), setops_laplace(), seta_laplace(), setb_laplace(), rom_solve_laplace(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup()

Call to:

loadbases(), pod(), set0flow(), vnorm(), snorm(), cnorm(), setab(), abm_shuffle(), orthonormb(), ps2b(), uip(), breduce(), pv2b(), sip(), vip(), cip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), regularev(), dump_serial(), gengram(), genevec()

subroutine  ps2k(ck, ux, uub)#

set snapshot coefficients for a given scalar basis

Parameters:

  • ck (1 + nb,ls) [real] :: coefficients

  • ux (lt,ls) [real] :: snapshots

  • uub (lt,1 + nb) [real] :: basis functions

Called from:

setops_laplace(), seta_laplace(), setb_laplace(), rom_solve_laplace()

Call to:

ps2b(), uip(), breduce(), pv2b(), sip(), vip(), cip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), regularev(), dump_serial(), gengram(), genevec()

subroutine  p2k(ck, usnap0, sb, mdim, wk)#

set snapshot coefficients for a given vector basis

Parameters:

  • ck (1 + nb,1) [real] :: coefficients

  • usnap0 (lt,mdim,ls) [real] :: snapshots (assume 0th mode is subtracted)

  • sb (lt,mdim,1 + nb) [real] :: basis functions

  • mdim [integer] :: 1 -> scalar, ndim -> vector

  • wk (ns) [real]

Called from:

rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc()

Call to:

uip(), breduce(), pv2b(), sip(), vip(), cip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), regularev(), dump_serial(), gengram(), genevec()

subroutine  pv2k(ck, usnap, uub, vvb, wwb)#

set snapshot coefficients for a given vector basis

Parameters:

  • ck (1 + nb,ls) [real] :: coefficients

  • usnap (lt,ldim,ls) [real] :: snapshots

  • uub (lt,1 + nb) [real] :: vvb,wwb, := x,y,z components of basis functions

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

Called from:

setops_laplace(), seta_laplace(), setb_laplace(), rom_solve_laplace()

Call to:

pv2b(), sip(), vip(), cip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

subroutine  ps2b(coef, tt, sb)#

get coordinates of a scalar field for a given basis

ck := coordinates of <ux> in <uub> ux := FOM scalar field uub := basis functions

Parameters:

  • coef (1 + nb) [real]

  • tt (lt) [real]

  • sb (lt,1 + nb) [real]

Called from:

rom_analysis(), fom_analysis(), sample_stats(), ce0n(), bases_analysis(), rom_sample(), snap_analysis(), pod_proj(), pod_df(), setdf(), set_les_imp(), setops_laplace(), seta_laplace(), setb_laplace(), rom_solve_laplace(), setr_laplace(), setbases(), ps2k(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu()

Call to:

sip(), vip(), cip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

subroutine  ps2b1(coef, tt, sb, nb2)#

get coordinates of a scalar field for a given basis w/o 0th mode

ck := coordinates of <ux> in <uub> ux := FOM scalar field uub := basis functions

Parameters:

  • coef (nb) [real]

  • tt (lt) [real]

  • sb (lt,nb) [real]

  • nb2 [integer]

Call to:

sip(), vip(), cip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

subroutine  pv2b(coef, ux, uy, uz, uub, vvb, wwb)#

get coordinates of a vector field for a given basis

Parameters:

  • coef (1 + nb) [real] :: coordinates of <ux,uy,uz> in <uub,vvb,wwb>

  • ux (lt) [real] :: uy,uz := x,y,z components of FOM field

  • uub (lt,1 + nb) [real] :: vvb,wwb := x,y,z components of basis functions

  • uy (lt) [real]

  • uz (lt) [real]

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

Called from:

rom_analysis(), fom_analysis(), sample_stats(), ce0n(), bases_analysis(), rom_sample(), snap_analysis(), pod_proj(), pod_df(), setdf(), set_les_imp(), setbases(), ps2k(), p2k(), pv2k(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu()

Call to:

vip(), cip(), sip(), wl2sip(), h10sip(), hlmsip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

subroutine  pc2b(cfu, cft, ux, uy, uz, tt, uub, vvb, wwb, ttb)#

get coordinates of a combined field for a given basis

Parameters:

  • cfu (1 + nb) [real] :: cft := coords of <ux,uy,uz,tt> in <uub,vvb,wwb,ttb>

  • ux (lt) [real] :: uy,uz,tt := vx,vy,vz,t components of combined FOM field

  • uub (lt,1 + nb) [real] :: vvb,wwb,ttb := vx,vy,vz,t components of basis functions

  • cft (1 + nb) [real]

  • uy (lt) [real]

  • uz (lt) [real]

  • tt (lt) [real]

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

  • ttb (lt,1 + nb) [real]

Called from:

rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu()

Call to:

cip(), sip(), wl2sip(), h10sip(), hlmsip(), vip(), wl2vip(), h10vip(), hlmvip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

function  sip(t1, t2)#

return inner-product of scalar fields

Parameters:

  • t1 (lt) [real] ::

    t2 := scalar fields ips := inner product type

    (L2 = L_2, H10 = H^1_0, HLM = Helmholtz)

  • t2 (lt) [real]

Return:

sip [real]

Called from:

opadd3(), recont(), recont_rms(), reconu_rms(), reconu_rm2(), reconv(), ctke_fom(), ctke(), lints(), ltruncr(), shift(), shift3(), lap2d(), lap3d(), push_op(), pop_op(), push_sol(), pop_sol(), zero_sol(), add_sol(), scale_sol(), copy_sol(), gradp(), setdps(), setconvbases(), setdtbases(), hyperpar(), opbinv1_nom(), binv1_nom(), binv1(), invmat(), cpsi(), csga(), emask(), set_surf(), set_gradn(), dump_sfld(), sol_intp_xline(), sol_intp_xline_qoi(), sol_intp_xline_fqoi(), interp_vec_prop2(), interp_vec2(), interp_sca_prop2(), interp_sca2(), count_gal(), cdump(), set_binv2(), diag(), checker(), checkera(), k_mean(), c_distortion_measure(), projtoprerb(), evaldut(), evalut(), reconv_wo0(), divm1(), recont_wo0(), svint(), checkaeq(), checkaeqp(), checkaeqd(), setdiff(), projvecm(), setpsi(), splitvec(), combvec(), tfunc(), evalf(), evalcflds(), evalaflds(), set0flow(), p2b(), p2b_helper(), fip(), find_cerr(), find_cerr2(), genrij(), orthonormb(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram(), genevec(), vnorm_(), vnorm(), cnorm(), snorm()

Call to:

wl2sip(), h10sip(), hlmsip(), vip(), wl2vip(), h10vip(), hlmvip(), cip(), sip(), h10sip_vd(), h10vip_vd(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

function  vip(t1, t2, t3, t4, t5, t6)#

return inner-product of vector fields

Parameters:

  • t1 (lt) [real] :: t2,t3 := x,y,z components of field 1

  • t4 (lt) [real] ::

    t5,t6 := x,y,z components of field 2 ips := inner product type

    (L2 = L_2, H10 = H^1_0, HLM = Helmholtz)

  • t2 (lt) [real]

  • t3 (lt) [real]

  • t5 (lt) [real]

  • t6 (lt) [real]

Return:

vip [real]

Called from:

opadd3(), recont(), recont_rms(), reconu_rms(), reconu_rm2(), reconv(), ctke_fom(), ctke(), lints(), ltruncr(), shift(), shift3(), lap2d(), lap3d(), push_op(), pop_op(), push_sol(), pop_sol(), zero_sol(), add_sol(), scale_sol(), copy_sol(), gradp(), setdps(), setconvbases(), setdtbases(), hyperpar(), opbinv1_nom(), binv1_nom(), binv1(), invmat(), cpsi(), csga(), emask(), set_surf(), set_gradn(), dump_sfld(), sol_intp_xline(), sol_intp_xline_qoi(), sol_intp_xline_fqoi(), interp_vec_prop2(), interp_vec2(), interp_sca_prop2(), interp_sca2(), count_gal(), cdump(), set_binv2(), diag(), checker(), checkera(), k_mean(), c_distortion_measure(), projtoprerb(), evaldut(), evalut(), reconv_wo0(), divm1(), recont_wo0(), svint(), checkaeq(), checkaeqp(), checkaeqd(), setdiff(), projvecm(), setpsi(), splitvec(), combvec(), tfunc(), evalf(), evalcflds(), evalaflds(), set0flow(), p2b(), p2b_helper(), fip(), find_cerr(), find_cerr2(), genrij(), orthonormb(), dump_global(), idump_serial(), dump_serial(), dump_global_helper(), idump_serial_helper(), dump_serial_helper(), dump_all(), dump_ops(), dump_bas(), dump_cbas(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram(), genevec(), vnorm_(), vnorm()

Call to:

wl2vip(), h10vip(), hlmvip(), cip(), vip(), sip(), h10sip_vd(), h10sip(), h10vip_vd(), hlmsip(), wl2sip(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), dump_serial(), gengram(), genevec()

function  cip(t1, t2, t3, t4, t5, t6, t7, t8)#

return inner-product of vector fields

Parameters:

  • t1 (lt) [real] :: t2,t3,t4 := vx,vy,vz,t components of field 1

  • t5 (lt) [real] :: t6,t7,t8 := vx,vy,vz,t components of field 2

  • t2 (lt) [real]

  • t3 (lt) [real]

  • t4 (lt) [real]

  • t6 (lt) [real]

  • t7 (lt) [real]

  • t8 (lt) [real]

Return:

cip [real]

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram(), genevec(), vnorm_(), vnorm(), cnorm()

Call to:

vip(), sip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), wl2sip(), hlmvip(), wl2vip(), wl2sip_vd(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  h10sip_vd(t1, t2, vd)#

return inner-product of scalar fields using the H^1_0 inner-product with an arbitrary diffusivity field

Parameters:

  • t1 (lt) [real] :: t2 := scalar fields

  • vd (lt) [real] :: variable diffusivity

  • t2 (lt) [real]

Return:

h10sip_vd [real]

Called from:

setops_laplace(), seta_laplace(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae()

Call to:

h10sip(), h10vip(), h10vip_vd(), hlmsip(), wl2sip(), hlmvip(), wl2vip(), wl2sip_vd(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  h10sip(t1, t2)#

return inner-product of scalar fields using the H^1_0 inner-product

Parameters:

  • t1 (lt) [real] :: t2 := scalar fields

  • t2 (lt) [real]

Return:

h10sip [real]

Called from:

cres(), c_avei(), set_sigma(), set_nres(), read_sigma(), c_sigma(), set_theta_poisson(), set_theta_heat(), set_theta_ad(), set_theta_ns(), rom_poisson(), setr_poisson(), set_alphaj(), set_betaj(), set_residual(), csigma_u(), csigma_t(), set_sigma_new(), crd(), cres_new(), cres_new1(), crd_test(), steady_stokes_solve(), crd_divf(), set_residual_unsteady(), set_theta_uns(), cres_uns(), c_rieszrd_uns(), set_alphaj_new(), set_betaj_new(), set_rom_residual(), get_dual_norm(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), set_xi_poisson(), set_xi_heat(), set_xi_ad(), set_xi_ns(), set_rr_ns(), set_sns_divfrr(), set_uns_divfrr(), set_riesz_reps(), set_riesz_one(), set_xi_a(), set_xi_b(), set_xi_c(), binva(), binvb(), binvc(), binv(), c_riesz_norm(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta()

Call to:

h10vip(), h10vip_vd(), hlmsip(), h10sip(), wl2sip(), hlmvip(), wl2vip(), wl2sip_vd(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  h10vip(t1, t2, t3, t4, t5, t6)#

return inner-product of vector fields using the H^1_0 inner-product

Parameters:

  • t1 (lt) [real] :: t2,t3; t4,t5,t6 := x,y,z components of vector fields

  • t2 (lt) [real]

  • t3 (lt) [real]

  • t4 (lt) [real]

  • t5 (lt) [real]

  • t6 (lt) [real]

Return:

h10vip [real]

Called from:

cres(), c_avei(), set_sigma(), set_nres(), read_sigma(), c_sigma(), set_theta_poisson(), set_theta_heat(), set_theta_ad(), set_theta_ns(), rom_poisson(), setr_poisson(), set_alphaj(), set_betaj(), set_residual(), csigma_u(), csigma_t(), set_sigma_new(), crd(), cres_new(), cres_new1(), crd_test(), steady_stokes_solve(), crd_divf(), set_residual_unsteady(), set_theta_uns(), cres_uns(), c_rieszrd_uns(), set_alphaj_new(), set_betaj_new(), set_rom_residual(), get_dual_norm(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), set_xi_poisson(), set_xi_heat(), set_xi_ad(), set_xi_ns(), set_rr_ns(), set_sns_divfrr(), set_uns_divfrr(), set_riesz_reps(), set_riesz_one(), set_xi_a(), set_xi_b(), set_xi_c(), binva(), binvb(), binvc(), binv(), c_riesz_norm(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta()

Call to:

h10vip_vd(), hlmsip(), h10sip(), wl2sip(), hlmvip(), h10vip(), wl2vip(), wl2sip_vd(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  h10vip_vd(t1, t2, t3, t4, t5, t6, vd)#

return inner-product of vector fields using the H^1_0 inner-product with arbitrary diffusivity fields

Parameters:

  • t1 (lt) [real] :: t2,t3; t4,t5,t6 := x,y,z components of vector fields

  • vd (lt,ldim) [real] :: variable diffusivity in multiple dimensions

  • t2 (lt) [real]

  • t3 (lt) [real]

  • t4 (lt) [real]

  • t5 (lt) [real]

  • t6 (lt) [real]

Return:

h10vip_vd [real]

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae()

Call to:

hlmsip(), h10sip(), wl2sip(), hlmvip(), h10vip(), wl2vip(), wl2sip_vd(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  hlmsip(t1, t2)#

return inner-product of scalar fields using the Helmholtz inner-product

Parameters:

  • t1 (lt) [real] :: t2 := scalar fields

  • t2 (lt) [real]

Return:

hlmsip [real]

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd()

Call to:

h10sip(), wl2sip(), hlmvip(), h10vip(), wl2vip(), wl2sip_vd(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  hlmvip(t1, t2, t3, t4, t5, t6)#

return inner-product of vector fields using the Helmholtz inner-product

Parameters:

  • t1 (lt) [real] :: t2,t3; t4,t5,t6 := x,y,z components of vector fields

  • t2 (lt) [real]

  • t3 (lt) [real]

  • t4 (lt) [real]

  • t5 (lt) [real]

  • t6 (lt) [real]

Return:

hlmvip [real]

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip()

Call to:

h10vip(), wl2vip(), wl2sip_vd(), wl2sip(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  wl2sip_vd(t1, t2, rho)#

return inner-product of scalar fields using the L^2 inner-product

Parameters:

  • t1 (lt) [real] :: t2 := scalar fields

  • t2 (lt) [real]

  • rho (lt) [real]

Return:

wl2sip_vd [real]

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb()

Call to:

wl2sip(), wl2vip(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  wl2sip(t1, t2)#

return inner-product of scalar fields using the L^2 inner-product

Parameters:

  • t1 (lt) [real] :: t2 := scalar fields

  • t2 (lt) [real]

Return:

wl2sip [real]

Called from:

cres(), c_avei(), set_sigma(), set_nres(), read_sigma(), c_sigma(), set_theta_poisson(), set_theta_heat(), set_theta_ad(), set_theta_ns(), rom_poisson(), setr_poisson(), set_alphaj(), set_betaj(), set_residual(), csigma_u(), csigma_t(), set_sigma_new(), crd(), cres_new(), cres_new1(), crd_test(), steady_stokes_solve(), crd_divf(), set_residual_unsteady(), set_theta_uns(), cres_uns(), c_rieszrd_uns(), set_alphaj_new(), set_betaj_new(), set_rom_residual(), get_dual_norm(), setops_laplace(), seta_laplace(), setb_laplace(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), set_xi_poisson(), set_xi_heat(), set_xi_ad(), set_xi_ns(), set_rr_ns(), set_sns_divfrr(), set_uns_divfrr(), set_riesz_reps(), set_riesz_one(), set_xi_a(), set_xi_b(), set_xi_c(), binva(), binvb(), binvc(), binv(), c_riesz_norm(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu(), setf(), setupvp(), setfluc()

Call to:

wl2vip(), sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

function  wl2vip(t1, t2, t3, t4, t5, t6)#

return inner-product of vector fields using the L^2 inner-product

Parameters:

  • t1 (lt) [real] :: t2,t3; t4,t5,t6 := x,y,z components of vector fields

  • t2 (lt) [real]

  • t3 (lt) [real]

  • t4 (lt) [real]

  • t5 (lt) [real]

  • t6 (lt) [real]

Return:

wl2vip [real]

Called from:

cres(), c_avei(), set_sigma(), set_nres(), read_sigma(), c_sigma(), set_theta_poisson(), set_theta_heat(), set_theta_ad(), set_theta_ns(), rom_poisson(), setr_poisson(), set_alphaj(), set_betaj(), set_residual(), csigma_u(), csigma_t(), set_sigma_new(), crd(), cres_new(), cres_new1(), crd_test(), steady_stokes_solve(), crd_divf(), set_residual_unsteady(), set_theta_uns(), cres_uns(), c_rieszrd_uns(), set_alphaj_new(), set_betaj_new(), set_rom_residual(), get_dual_norm(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), set_xi_poisson(), set_xi_heat(), set_xi_ad(), set_xi_ns(), set_rr_ns(), set_sns_divfrr(), set_uns_divfrr(), set_riesz_reps(), set_riesz_one(), set_xi_a(), set_xi_b(), set_xi_c(), binva(), binvb(), binvc(), binv(), c_riesz_norm(), rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu(), setf(), setupvp(), setfluc()

Call to:

sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

subroutine  hlmgg(gram, s, ms, mdim)#

set the Gramian based on the Helmholtz inner-product

Parameters:

  • gram (ms,ms) [real] :: Gramian

  • s (lt,mdim,ms) [real] :: snapshots

  • ms [integer] :: Form the Gramian, U=U_K^T A U_K using L2 Norm

  • mdim [integer] :: vector dimension

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram()

Call to:

sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

subroutine  h10gg(gram, s, ms, mdim)#

set the Gramian based on the H^1_0 inner-product

Parameters:

  • gram (ms,ms) [real] :: Gramian

  • s (lt,mdim,ms) [real] :: snapshots

  • ms [integer] :: Form the Gramian, U=U_K^T A U_K using L2 Norm

  • mdim [integer] :: vector dimension

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram()

Call to:

sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

subroutine  wl2gg(gram, s, ms, mdim)#

set the Gramian based on the L^2 inner-product (duplicate, will be deprecated)

Parameters:

  • gram (ms,ms) [real] :: Gramian

  • s (lt,mdim,ms) [real] :: snapshots

  • ms [integer] :: Form the Gramian, U=U_K^T A U_K using L2 Norm

  • mdim [integer] :: vector dimension

Call to:

sip(), vip(), gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), cip(), dump_serial(), gengram(), genevec()

subroutine  gengraml2(gram, s, ms, mdim)#

set the Gramian based on the L^2 inner-product

Parameters:

  • gram (ms,ms) [real] :: Gramian

  • s (lt,mdim,ms) [real] :: snapshots

  • ms [integer] :: Form the Gramian, U=U_K^T A U_K using L2 Norm

  • mdim [integer] :: vector dimension

Called from:

setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram()

Call to:

gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), vip(), cip(), sip(), dump_serial(), gengram(), genevec()

subroutine  gengram(gram, s, ms, mdim, cips, nbat)#

set the Gramian based on the inner-product set by ips

Parameters:

  • gram (1) [real] :: Gramian

  • s (1) [real] :: snapshots

  • ms [integer] :: number of snapshots

  • mdim [integer] :: vector dimension

  • cips [character] :: inner-product space specifier

  • nbat [integer] :: number of inner-products in batch

Called from:

rom_analysis(), fom_analysis(), split_domain(), domain_mask(), domain_mask_p(), pod_part(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram(), genevec(), vnorm_(), vnorm(), cnorm(), snorm(), h10pv2b(), hlmpv2b(), regularev(), cnmax(), cenpm(), pod()

Call to:

gengraml2(), h10gg(), hlmgg(), breduce(), regularev(), vip(), cip(), sip(), dump_serial(), gengram(), genevec()

subroutine  genevec(vec, val, ms, mb, ifld)#

solve eigensystem based on the given Gramian (vec)

Parameters:

  • vec (ms,ms) [real] :: eigenvectors (initially Gramian)

  • val (ms) [real] :: eigenvalues

  • ms [integer] :: number of snapshots

  • mb [integer] :: number of basis

  • ifld [integer] :: field number

Called from:

rom_analysis(), fom_analysis(), split_domain(), domain_mask(), domain_mask_p(), pod_part(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram(), genevec(), vnorm_(), vnorm(), cnorm(), snorm(), h10pv2b(), hlmpv2b(), regularev(), cnmax(), cenpm(), pod()

Call to:

regularev(), vip(), cip(), sip(), dump_serial(), gengram(), genevec()

subroutine  vnorm_(uvwbb)#

Normalizes the vector field uvwbb.

Parameters:

uvwbb (lt,ldim,1 + nb) [real] :: The vector field of shape (lt,ldim,nb)

Call to:

vip(), cip(), sip(), dump_serial(), gengram(), genevec()

subroutine  vnorm(uub, vvb, wwb)#

normalizes vector field

Parameters:

  • uub (lt,1 + nb) [real] :: vvb,wwb := x,y,z components of vector field

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

Called from:

setbases()

Call to:

vip(), cip(), sip(), dump_serial(), gengram(), genevec()

subroutine  cnorm(uub, vvb, wwb, ttb)#

normalizes combined field

Parameters:

  • uub (lt,1 + nb) [real] :: vvb,wwb,ttb := vx,vy,vz,t components of vector field

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

  • ttb (lt,1 + nb) [real]

Called from:

setbases()

Call to:

cip(), sip(), dump_serial(), gengram(), genevec()

subroutine  snorm(ssb)#

normalizes scalar field

Parameters:

ssb (lt,1 + nb) [real] :: scalar field

Called from:

setops_laplace(), seta_laplace(), setb_laplace(), rom_solve_laplace(), setbases()

Call to:

sip(), dump_serial(), gengram(), genevec()

subroutine  h10pv2b(coef, ux, uy, uz, uub, vvb, wwb)#

get coordinates of a vector field for a H^1_0 orthogonal basis

Parameters:

  • coef (1 + nb) [real] :: coordinates of <ux,uy,uz> in <uub,vvb,wwb>

  • ux (lt) [real] :: uy,uz := x,y,z components of FOM field

  • uub (lt,1 + nb) [real] :: vvb,wwb := x,y,z components of basis functions

  • uy (lt) [real]

  • uz (lt) [real]

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

Called from:

rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu()

Call to:

dump_serial(), gengram(), genevec()

subroutine  hlmpv2b(coef, ux, uy, uz, uub, vvb, wwb)#

get coordinates of a vector field for a Helmholtz orthogonal basis

Parameters:

  • coef (1 + nb) [real] :: coordinates of <ux,uy,uz> in <uub,vvb,wwb>

  • ux (lt) [real] :: uy,uz := x,y,z components of FOM field

  • uub (lt,1 + nb) [real] :: vvb,wwb := x,y,z components of basis functions

  • uy (lt) [real]

  • uz (lt) [real]

  • vvb (lt,1 + nb) [real]

  • wwb (lt,1 + nb) [real]

Called from:

rom_update(), offline_mode(), online_mode(), recon_mode(), rom_setup(), asnap(), setops(), setqoi(), update_k(), update_hyper(), setmisc(), mor_init_params(), mor_set_params_rea(), mor_set_params_uni_pre(), mor_set_params_uni_post(), mor_show_params(), mor_init_fields(), setc(), seta(), sets(), setb(), setae(), setu()

Call to:

dump_serial(), gengram(), genevec()

subroutine  regularev(a, lam, n, wk, lwork)#

Solve the eigenvalue problem A x = lam x

A := symmetric matrix

“SIZE” is included here only to deduce WDSIZE, the working precision, in bytes, so as to know whether dsygv or ssygv should be called.

Parameters:

  • a (n,n) [real]

  • lam (n) [real]

  • n [integer]

  • wk (n,n) [real]

  • lwork [integer]

Called from:

opadd3(), recont(), recont_rms(), reconu_rms(), reconu_rm2(), reconv(), ctke_fom(), ctke(), lints(), ltruncr(), shift(), shift3(), lap2d(), lap3d(), push_op(), pop_op(), push_sol(), pop_sol(), zero_sol(), add_sol(), scale_sol(), copy_sol(), gradp(), setdps(), setconvbases(), setdtbases(), hyperpar(), opbinv1_nom(), binv1_nom(), binv1(), invmat(), cpsi(), csga(), emask(), set_surf(), set_gradn(), dump_sfld(), sol_intp_xline(), sol_intp_xline_qoi(), sol_intp_xline_fqoi(), interp_vec_prop2(), interp_vec2(), interp_sca_prop2(), interp_sca2(), count_gal(), cdump(), set_binv2(), diag(), setbases(), ps2k(), p2k(), pv2k(), ps2b(), ps2b1(), pv2b(), pc2b(), sip(), vip(), cip(), h10sip_vd(), h10sip(), h10vip(), h10vip_vd(), hlmsip(), hlmvip(), wl2sip_vd(), wl2sip(), wl2vip(), hlmgg(), h10gg(), wl2gg(), gengraml2(), gengram(), genevec()

Call to:

dump_serial(), gengram(), genevec()

subroutine  cnmax(val, fname, ifld)#

compute maximum number of POD modes to be used

Parameters:

  • val (ls) [real] :: eigenvalues based on the generated Gramians

  • fname [character] :: filename that will be created

  • ifld [integer] :: field number

Call to:

dump_serial(), gengram(), genevec()

subroutine  cenpm(val, fname, ifld)#

Compute the percentage that each POD mode represents in the total averaged energy. The energy is defined based on the inner-product set by ips

Parameters:

  • val (ls) [real] :: eigenvalues based on the generated Gramians

  • fname [character] :: filename that will be created

  • ifld [integer] :: field number

Call to:

dump_serial(), gengram(), genevec()

subroutine  pod(basis, eval, gram, snaps, mdim, cips, nb, ns, ifpod, cop, nbat)#
Parameters:

  • basis (lt,mdim,nb) [real]

  • eval (ns) [real]

  • gram (ns,ns) [real]

  • snaps (lt,mdim,ns) [real]

  • mdim [integer]

  • cips [character]

  • nb [integer]

  • ns [integer]

  • ifpod [logical]

  • cop [character]

  • nbat [integer]

Called from:

dump_global(), idump_serial(), dump_serial(), dump_global_helper(), idump_serial_helper(), dump_serial_helper(), dump_all(), dump_ops(), dump_bas(), dump_cbas(), setbases()

Call to:

gengram(), dump_serial(), genevec()