70 character(LEN=13),
dimension(0:QUAD_TOT_NB) ::
quad_name 97 logical,
intent(in) :: b
99 integer :: qt, geo, dim, nn
101 if (b)
write(*,*)
"quad_mod : quad_init" 223 quad_coord(qt)%y(1,1) = 0.21132486540518702_rp
227 quad_wgt(qt)%y(1) = 0.44444444444444444444444444444444444_rp
228 quad_wgt(qt)%y(2) = 0.27777777777777777777777777777777778_rp
229 quad_wgt(qt)%y(3) = 0.27777777777777777777777777777777778_rp
231 quad_coord(qt)%y(1,2) = 0.11270166537925852_rp
235 quad_wgt(qt)%y(1) = 0.17392742256872692868653197461099970_rp
236 quad_wgt(qt)%y(2) = 0.17392742256872692868653197461099970_rp
237 quad_wgt(qt)%y(3) = 0.32607257743127307131346802538900030_rp
238 quad_wgt(qt)%y(4) = 0.32607257743127307131346802538900030_rp
239 a = 0.06943184420297371238802675555359525_rp
242 a = 0.33000947820757186759866712044837766_rp
258 & trim(
choral_dir)//
'SRC/quad_data/strang5_w.txt')
261 & trim(
choral_dir)//
'SRC/quad_data/strang5_x.txt')
266 & trim(
choral_dir)//
'SRC/quad_data/strang9_w.txt')
269 & trim(
choral_dir)//
'SRC/quad_data/strang9_x.txt')
273 & trim(
choral_dir)//
'SRC/quad_data/strang10_w.txt')
276 & trim(
choral_dir)//
'SRC/quad_data/strang10_x.txt')
280 & trim(
choral_dir)//
'SRC/quad_data/TOMS584_19_w.txt')
283 & trim(
choral_dir)//
'SRC/quad_data/TOMS584_19_x.txt')
288 & trim(
choral_dir)//
'SRC/quad_data/toms612_28_w.txt')
291 & trim(
choral_dir)//
'SRC/quad_data/toms612_28_x.txt')
295 & trim(
choral_dir)//
'SRC/quad_data/toms706_37_w.txt')
298 & trim(
choral_dir)//
'SRC/quad_data/toms706_37_x.txt')
307 & trim(
choral_dir)//
'SRC/quad_data/keast1_w.txt')
310 & trim(
choral_dir)//
'SRC/quad_data/keast1_x.txt')
314 & trim(
choral_dir)//
'SRC/quad_data/keast6_w.txt')
317 & trim(
choral_dir)//
'SRC/quad_data/keast6_x.txt')
321 & trim(
choral_dir)//
'SRC/quad_data/keast8_w.txt')
324 & trim(
choral_dir)//
'SRC/quad_data/keast8_x.txt')
328 & trim(
choral_dir)//
'SRC/quad_data/keast9_w.txt')
331 & trim(
choral_dir)//
'SRC/quad_data/keast9_x.txt')
integer, parameter quad_gauss_trg_28
integer, parameter cell_tet
Tetrahedron.
integer, parameter quad_gauss_tet_15
type(r_1d), dimension(quad_tot_nb), public quad_wgt
quad weights
integer, parameter quad_gauss_edg_4
QUADRATURE RULES ON REFERENCE CELLS
subroutine def_quad_data()
integer, parameter quad_gauss_tet_4
integer, parameter quad_gauss_trg_3
integer, parameter quad_none
subroutine, public quad_init(b)
initialise QUAD_XXX arrays:
integer, parameter cell_edg
Edge (line segment)
integer, parameter quad_gauss_tet_1
integer, dimension(cell_tot_nb), public cell_dim
Dimension for each cell type.
integer, parameter quad_gauss_tet_45
integer, dimension(quad_tot_nb), public quad_dim
Reference cell dimension for each quad method.
conversion integers or rational to real
integer, dimension(quad_tot_nb), public quad_order
Order for each quad method.
character(len=100), parameter, public choral_dir
CHORAL_DIR = path to choral.
REAL NUMBERS PRECISION IN CHORAL: selects simple/double/quad
IO: module for input/output
type(r_2d), dimension(quad_tot_nb), public quad_coord
quad node coordinates
integer, parameter quad_gauss_edg_3
integer, parameter quad_tot_nb
Total number of quadrature rules.
integer, parameter cell_trg
Triangle.
integer, dimension(quad_tot_nb), public quad_nbnodes
Number of nodes for each quad method.
integer, parameter quad_gauss_trg_1
subroutine def_quad_xxx()
R_1D: this type will allow to define arrays of 1D real arrays.
integer, parameter quad_gauss_trg_19
integer, parameter quad_gauss_trg_37
DEFINITION OF DIRECTORY PATHS
integer, parameter quad_gauss_edg_2
integer, parameter quad_gauss_trg_13
integer, parameter quad_gauss_edg_1
integer, parameter quad_gauss_tet_31
character(len=13), dimension(0:quad_tot_nb), public quad_name
QUAD_ARRAY.
R_2D: this type will allow to define arrays of 2D real arrays.
integer, parameter quad_gauss_trg_6
integer, parameter quad_gauss_trg_12
DEFINITION OF GEOMETRICAL CELLS (for meshes)
integer, dimension(quad_tot_nb), public quad_geo
Reference cell geometry for each quad method.