quadrature/s_quadrature_functions.f90
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | !> @file s_quadrature_functions.f90 | ||
| 2 | !> @brief Implementations of Gauss-Legendre quadrature functions | ||
| 3 | !> | ||
| 4 | !> @note This file is auto-generated by quadrature_generate.jl. Do not edit manually. | ||
| 5 | submodule(m_quadrature_functions) s_quadrature_functions | ||
| 6 | implicit none | ||
| 7 | |||
| 8 | contains | ||
| 9 | |||
| 10 | 42564 | module procedure quadrature_1 | |
| 11 | real(wp) :: h, xmid | ||
| 12 | |||
| 13 | 42564 | h = (xright - xleft) / 2 | |
| 14 | 42564 | xmid = (xleft + xright) / 2 | |
| 15 | |||
| 16 | 42564 | ans = gl_weights_1(1) * f(xmid + h * gl_nodes_1(1)) | |
| 17 | 42564 | ans = ans * h | |
| 18 | 42564 | end procedure | |
| 19 | |||
| 20 | 70771 | module procedure quadrature_2 | |
| 21 | real(wp) :: h, xmid | ||
| 22 | |||
| 23 | 70771 | h = (xright - xleft) / 2 | |
| 24 | 70771 | xmid = (xleft + xright) / 2 | |
| 25 | |||
| 26 | ans = gl_weights_2(1) * f(xmid + h * gl_nodes_2(1)) & | ||
| 27 | 70771 | & + gl_weights_2(2) * f(xmid + h * gl_nodes_2(2)) | |
| 28 | 70771 | ans = ans * h | |
| 29 | 70771 | end procedure | |
| 30 | |||
| 31 | 139646 | module procedure quadrature_3 | |
| 32 | real(wp) :: h, xmid | ||
| 33 | |||
| 34 | 139646 | h = (xright - xleft) / 2 | |
| 35 | 139646 | xmid = (xleft + xright) / 2 | |
| 36 | |||
| 37 | ans = gl_weights_3(1) * f(xmid + h * gl_nodes_3(1)) & | ||
| 38 | & + gl_weights_3(2) * f(xmid + h * gl_nodes_3(2)) & | ||
| 39 | 139646 | & + gl_weights_3(3) * f(xmid + h * gl_nodes_3(3)) | |
| 40 | 139646 | ans = ans * h | |
| 41 | 139646 | end procedure | |
| 42 | |||
| 43 | 638075 | module procedure quadrature_4 | |
| 44 | real(wp) :: h, xmid | ||
| 45 | |||
| 46 | 638075 | h = (xright - xleft) / 2 | |
| 47 | 638075 | xmid = (xleft + xright) / 2 | |
| 48 | |||
| 49 | ans = gl_weights_4(1) * f(xmid + h * gl_nodes_4(1)) & | ||
| 50 | & + gl_weights_4(2) * f(xmid + h * gl_nodes_4(2)) & | ||
| 51 | & + gl_weights_4(3) * f(xmid + h * gl_nodes_4(3)) & | ||
| 52 | 638075 | & + gl_weights_4(4) * f(xmid + h * gl_nodes_4(4)) | |
| 53 | 638075 | ans = ans * h | |
| 54 | 638075 | end procedure | |
| 55 | |||
| 56 | 2505391 | module procedure quadrature_5 | |
| 57 | real(wp) :: h, xmid | ||
| 58 | |||
| 59 | 2505391 | h = (xright - xleft) / 2 | |
| 60 | 2505391 | xmid = (xleft + xright) / 2 | |
| 61 | |||
| 62 | ans = gl_weights_5(1) * f(xmid + h * gl_nodes_5(1)) & | ||
| 63 | & + gl_weights_5(2) * f(xmid + h * gl_nodes_5(2)) & | ||
| 64 | & + gl_weights_5(3) * f(xmid + h * gl_nodes_5(3)) & | ||
| 65 | & + gl_weights_5(4) * f(xmid + h * gl_nodes_5(4)) & | ||
| 66 | 2505391 | & + gl_weights_5(5) * f(xmid + h * gl_nodes_5(5)) | |
| 67 | 2505391 | ans = ans * h | |
| 68 | 2505391 | end procedure | |
| 69 | |||
| 70 | 65161 | module procedure quadrature_6 | |
| 71 | real(wp) :: h, xmid | ||
| 72 | |||
| 73 | 65161 | h = (xright - xleft) / 2 | |
| 74 | 65161 | xmid = (xleft + xright) / 2 | |
| 75 | |||
| 76 | ans = gl_weights_6(1) * f(xmid + h * gl_nodes_6(1)) & | ||
| 77 | & + gl_weights_6(2) * f(xmid + h * gl_nodes_6(2)) & | ||
| 78 | & + gl_weights_6(3) * f(xmid + h * gl_nodes_6(3)) & | ||
| 79 | & + gl_weights_6(4) * f(xmid + h * gl_nodes_6(4)) & | ||
| 80 | & + gl_weights_6(5) * f(xmid + h * gl_nodes_6(5)) & | ||
| 81 | 65161 | & + gl_weights_6(6) * f(xmid + h * gl_nodes_6(6)) | |
| 82 | 65161 | ans = ans * h | |
| 83 | 65161 | end procedure | |
| 84 | |||
| 85 | 76389 | module procedure quadrature_7 | |
| 86 | real(wp) :: h, xmid | ||
| 87 | |||
| 88 | 76389 | h = (xright - xleft) / 2 | |
| 89 | 76389 | xmid = (xleft + xright) / 2 | |
| 90 | |||
| 91 | ans = gl_weights_7(1) * f(xmid + h * gl_nodes_7(1)) & | ||
| 92 | & + gl_weights_7(2) * f(xmid + h * gl_nodes_7(2)) & | ||
| 93 | & + gl_weights_7(3) * f(xmid + h * gl_nodes_7(3)) & | ||
| 94 | & + gl_weights_7(4) * f(xmid + h * gl_nodes_7(4)) & | ||
| 95 | & + gl_weights_7(5) * f(xmid + h * gl_nodes_7(5)) & | ||
| 96 | & + gl_weights_7(6) * f(xmid + h * gl_nodes_7(6)) & | ||
| 97 | 76389 | & + gl_weights_7(7) * f(xmid + h * gl_nodes_7(7)) | |
| 98 | 76389 | ans = ans * h | |
| 99 | 76389 | end procedure | |
| 100 | |||
| 101 | 98079 | module procedure quadrature_8 | |
| 102 | real(wp) :: h, xmid | ||
| 103 | |||
| 104 | 98079 | h = (xright - xleft) / 2 | |
| 105 | 98079 | xmid = (xleft + xright) / 2 | |
| 106 | |||
| 107 | ans = gl_weights_8(1) * f(xmid + h * gl_nodes_8(1)) & | ||
| 108 | & + gl_weights_8(2) * f(xmid + h * gl_nodes_8(2)) & | ||
| 109 | & + gl_weights_8(3) * f(xmid + h * gl_nodes_8(3)) & | ||
| 110 | & + gl_weights_8(4) * f(xmid + h * gl_nodes_8(4)) & | ||
| 111 | & + gl_weights_8(5) * f(xmid + h * gl_nodes_8(5)) & | ||
| 112 | & + gl_weights_8(6) * f(xmid + h * gl_nodes_8(6)) & | ||
| 113 | & + gl_weights_8(7) * f(xmid + h * gl_nodes_8(7)) & | ||
| 114 | 98079 | & + gl_weights_8(8) * f(xmid + h * gl_nodes_8(8)) | |
| 115 | 98079 | ans = ans * h | |
| 116 | 98079 | end procedure | |
| 117 | |||
| 118 | 121759 | module procedure quadrature_9 | |
| 119 | real(wp) :: h, xmid | ||
| 120 | |||
| 121 | 121759 | h = (xright - xleft) / 2 | |
| 122 | 121759 | xmid = (xleft + xright) / 2 | |
| 123 | |||
| 124 | ans = gl_weights_9(1) * f(xmid + h * gl_nodes_9(1)) & | ||
| 125 | & + gl_weights_9(2) * f(xmid + h * gl_nodes_9(2)) & | ||
| 126 | & + gl_weights_9(3) * f(xmid + h * gl_nodes_9(3)) & | ||
| 127 | & + gl_weights_9(4) * f(xmid + h * gl_nodes_9(4)) & | ||
| 128 | & + gl_weights_9(5) * f(xmid + h * gl_nodes_9(5)) & | ||
| 129 | & + gl_weights_9(6) * f(xmid + h * gl_nodes_9(6)) & | ||
| 130 | & + gl_weights_9(7) * f(xmid + h * gl_nodes_9(7)) & | ||
| 131 | & + gl_weights_9(8) * f(xmid + h * gl_nodes_9(8)) & | ||
| 132 | 121759 | & + gl_weights_9(9) * f(xmid + h * gl_nodes_9(9)) | |
| 133 | 121759 | ans = ans * h | |
| 134 | 121759 | end procedure | |
| 135 | |||
| 136 | 407 | module procedure quadrature_10 | |
| 137 | real(wp) :: h, xmid | ||
| 138 | |||
| 139 | 407 | h = (xright - xleft) / 2 | |
| 140 | 407 | xmid = (xleft + xright) / 2 | |
| 141 | |||
| 142 | ans = gl_weights_10(1) * f(xmid + h * gl_nodes_10(1)) & | ||
| 143 | & + gl_weights_10(2) * f(xmid + h * gl_nodes_10(2)) & | ||
| 144 | & + gl_weights_10(3) * f(xmid + h * gl_nodes_10(3)) & | ||
| 145 | & + gl_weights_10(4) * f(xmid + h * gl_nodes_10(4)) & | ||
| 146 | & + gl_weights_10(5) * f(xmid + h * gl_nodes_10(5)) & | ||
| 147 | & + gl_weights_10(6) * f(xmid + h * gl_nodes_10(6)) & | ||
| 148 | & + gl_weights_10(7) * f(xmid + h * gl_nodes_10(7)) & | ||
| 149 | & + gl_weights_10(8) * f(xmid + h * gl_nodes_10(8)) & | ||
| 150 | & + gl_weights_10(9) * f(xmid + h * gl_nodes_10(9)) & | ||
| 151 | 407 | & + gl_weights_10(10) * f(xmid + h * gl_nodes_10(10)) | |
| 152 | 407 | ans = ans * h | |
| 153 | 407 | end procedure | |
| 154 | |||
| 155 | 451 | module procedure quadrature_11 | |
| 156 | real(wp) :: h, xmid | ||
| 157 | |||
| 158 | 451 | h = (xright - xleft) / 2 | |
| 159 | 451 | xmid = (xleft + xright) / 2 | |
| 160 | |||
| 161 | ans = gl_weights_11(1) * f(xmid + h * gl_nodes_11(1)) & | ||
| 162 | & + gl_weights_11(2) * f(xmid + h * gl_nodes_11(2)) & | ||
| 163 | & + gl_weights_11(3) * f(xmid + h * gl_nodes_11(3)) & | ||
| 164 | & + gl_weights_11(4) * f(xmid + h * gl_nodes_11(4)) & | ||
| 165 | & + gl_weights_11(5) * f(xmid + h * gl_nodes_11(5)) & | ||
| 166 | & + gl_weights_11(6) * f(xmid + h * gl_nodes_11(6)) & | ||
| 167 | & + gl_weights_11(7) * f(xmid + h * gl_nodes_11(7)) & | ||
| 168 | & + gl_weights_11(8) * f(xmid + h * gl_nodes_11(8)) & | ||
| 169 | & + gl_weights_11(9) * f(xmid + h * gl_nodes_11(9)) & | ||
| 170 | & + gl_weights_11(10) * f(xmid + h * gl_nodes_11(10)) & | ||
| 171 | 451 | & + gl_weights_11(11) * f(xmid + h * gl_nodes_11(11)) | |
| 172 | 451 | ans = ans * h | |
| 173 | 451 | end procedure | |
| 174 | |||
| 175 | 1 | module procedure quadrature_12 | |
| 176 | real(wp) :: h, xmid | ||
| 177 | |||
| 178 | 1 | h = (xright - xleft) / 2 | |
| 179 | 1 | xmid = (xleft + xright) / 2 | |
| 180 | |||
| 181 | ans = gl_weights_12(1) * f(xmid + h * gl_nodes_12(1)) & | ||
| 182 | & + gl_weights_12(2) * f(xmid + h * gl_nodes_12(2)) & | ||
| 183 | & + gl_weights_12(3) * f(xmid + h * gl_nodes_12(3)) & | ||
| 184 | & + gl_weights_12(4) * f(xmid + h * gl_nodes_12(4)) & | ||
| 185 | & + gl_weights_12(5) * f(xmid + h * gl_nodes_12(5)) & | ||
| 186 | & + gl_weights_12(6) * f(xmid + h * gl_nodes_12(6)) & | ||
| 187 | & + gl_weights_12(7) * f(xmid + h * gl_nodes_12(7)) & | ||
| 188 | & + gl_weights_12(8) * f(xmid + h * gl_nodes_12(8)) & | ||
| 189 | & + gl_weights_12(9) * f(xmid + h * gl_nodes_12(9)) & | ||
| 190 | & + gl_weights_12(10) * f(xmid + h * gl_nodes_12(10)) & | ||
| 191 | & + gl_weights_12(11) * f(xmid + h * gl_nodes_12(11)) & | ||
| 192 | 1 | & + gl_weights_12(12) * f(xmid + h * gl_nodes_12(12)) | |
| 193 | 1 | ans = ans * h | |
| 194 | 1 | end procedure | |
| 195 | |||
| 196 | 73 | module procedure quadrature_13 | |
| 197 | real(wp) :: h, xmid | ||
| 198 | |||
| 199 | 73 | h = (xright - xleft) / 2 | |
| 200 | 73 | xmid = (xleft + xright) / 2 | |
| 201 | |||
| 202 | ans = gl_weights_13(1) * f(xmid + h * gl_nodes_13(1)) & | ||
| 203 | & + gl_weights_13(2) * f(xmid + h * gl_nodes_13(2)) & | ||
| 204 | & + gl_weights_13(3) * f(xmid + h * gl_nodes_13(3)) & | ||
| 205 | & + gl_weights_13(4) * f(xmid + h * gl_nodes_13(4)) & | ||
| 206 | & + gl_weights_13(5) * f(xmid + h * gl_nodes_13(5)) & | ||
| 207 | & + gl_weights_13(6) * f(xmid + h * gl_nodes_13(6)) & | ||
| 208 | & + gl_weights_13(7) * f(xmid + h * gl_nodes_13(7)) & | ||
| 209 | & + gl_weights_13(8) * f(xmid + h * gl_nodes_13(8)) & | ||
| 210 | & + gl_weights_13(9) * f(xmid + h * gl_nodes_13(9)) & | ||
| 211 | & + gl_weights_13(10) * f(xmid + h * gl_nodes_13(10)) & | ||
| 212 | & + gl_weights_13(11) * f(xmid + h * gl_nodes_13(11)) & | ||
| 213 | & + gl_weights_13(12) * f(xmid + h * gl_nodes_13(12)) & | ||
| 214 | 73 | & + gl_weights_13(13) * f(xmid + h * gl_nodes_13(13)) | |
| 215 | 73 | ans = ans * h | |
| 216 | 73 | end procedure | |
| 217 | |||
| 218 | 1 | module procedure quadrature_14 | |
| 219 | real(wp) :: h, xmid | ||
| 220 | |||
| 221 | 1 | h = (xright - xleft) / 2 | |
| 222 | 1 | xmid = (xleft + xright) / 2 | |
| 223 | |||
| 224 | ans = gl_weights_14(1) * f(xmid + h * gl_nodes_14(1)) & | ||
| 225 | & + gl_weights_14(2) * f(xmid + h * gl_nodes_14(2)) & | ||
| 226 | & + gl_weights_14(3) * f(xmid + h * gl_nodes_14(3)) & | ||
| 227 | & + gl_weights_14(4) * f(xmid + h * gl_nodes_14(4)) & | ||
| 228 | & + gl_weights_14(5) * f(xmid + h * gl_nodes_14(5)) & | ||
| 229 | & + gl_weights_14(6) * f(xmid + h * gl_nodes_14(6)) & | ||
| 230 | & + gl_weights_14(7) * f(xmid + h * gl_nodes_14(7)) & | ||
| 231 | & + gl_weights_14(8) * f(xmid + h * gl_nodes_14(8)) & | ||
| 232 | & + gl_weights_14(9) * f(xmid + h * gl_nodes_14(9)) & | ||
| 233 | & + gl_weights_14(10) * f(xmid + h * gl_nodes_14(10)) & | ||
| 234 | & + gl_weights_14(11) * f(xmid + h * gl_nodes_14(11)) & | ||
| 235 | & + gl_weights_14(12) * f(xmid + h * gl_nodes_14(12)) & | ||
| 236 | & + gl_weights_14(13) * f(xmid + h * gl_nodes_14(13)) & | ||
| 237 | 1 | & + gl_weights_14(14) * f(xmid + h * gl_nodes_14(14)) | |
| 238 | 1 | ans = ans * h | |
| 239 | 1 | end procedure | |
| 240 | |||
| 241 | 1 | module procedure quadrature_15 | |
| 242 | real(wp) :: h, xmid | ||
| 243 | |||
| 244 | 1 | h = (xright - xleft) / 2 | |
| 245 | 1 | xmid = (xleft + xright) / 2 | |
| 246 | |||
| 247 | ans = gl_weights_15(1) * f(xmid + h * gl_nodes_15(1)) & | ||
| 248 | & + gl_weights_15(2) * f(xmid + h * gl_nodes_15(2)) & | ||
| 249 | & + gl_weights_15(3) * f(xmid + h * gl_nodes_15(3)) & | ||
| 250 | & + gl_weights_15(4) * f(xmid + h * gl_nodes_15(4)) & | ||
| 251 | & + gl_weights_15(5) * f(xmid + h * gl_nodes_15(5)) & | ||
| 252 | & + gl_weights_15(6) * f(xmid + h * gl_nodes_15(6)) & | ||
| 253 | & + gl_weights_15(7) * f(xmid + h * gl_nodes_15(7)) & | ||
| 254 | & + gl_weights_15(8) * f(xmid + h * gl_nodes_15(8)) & | ||
| 255 | & + gl_weights_15(9) * f(xmid + h * gl_nodes_15(9)) & | ||
| 256 | & + gl_weights_15(10) * f(xmid + h * gl_nodes_15(10)) & | ||
| 257 | & + gl_weights_15(11) * f(xmid + h * gl_nodes_15(11)) & | ||
| 258 | & + gl_weights_15(12) * f(xmid + h * gl_nodes_15(12)) & | ||
| 259 | & + gl_weights_15(13) * f(xmid + h * gl_nodes_15(13)) & | ||
| 260 | & + gl_weights_15(14) * f(xmid + h * gl_nodes_15(14)) & | ||
| 261 | 1 | & + gl_weights_15(15) * f(xmid + h * gl_nodes_15(15)) | |
| 262 | 1 | ans = ans * h | |
| 263 | 1 | end procedure | |
| 264 | |||
| 265 | 1 | module procedure quadrature_16 | |
| 266 | real(wp) :: h, xmid | ||
| 267 | |||
| 268 | 1 | h = (xright - xleft) / 2 | |
| 269 | 1 | xmid = (xleft + xright) / 2 | |
| 270 | |||
| 271 | ans = gl_weights_16(1) * f(xmid + h * gl_nodes_16(1)) & | ||
| 272 | & + gl_weights_16(2) * f(xmid + h * gl_nodes_16(2)) & | ||
| 273 | & + gl_weights_16(3) * f(xmid + h * gl_nodes_16(3)) & | ||
| 274 | & + gl_weights_16(4) * f(xmid + h * gl_nodes_16(4)) & | ||
| 275 | & + gl_weights_16(5) * f(xmid + h * gl_nodes_16(5)) & | ||
| 276 | & + gl_weights_16(6) * f(xmid + h * gl_nodes_16(6)) & | ||
| 277 | & + gl_weights_16(7) * f(xmid + h * gl_nodes_16(7)) & | ||
| 278 | & + gl_weights_16(8) * f(xmid + h * gl_nodes_16(8)) & | ||
| 279 | & + gl_weights_16(9) * f(xmid + h * gl_nodes_16(9)) & | ||
| 280 | & + gl_weights_16(10) * f(xmid + h * gl_nodes_16(10)) & | ||
| 281 | & + gl_weights_16(11) * f(xmid + h * gl_nodes_16(11)) & | ||
| 282 | & + gl_weights_16(12) * f(xmid + h * gl_nodes_16(12)) & | ||
| 283 | & + gl_weights_16(13) * f(xmid + h * gl_nodes_16(13)) & | ||
| 284 | & + gl_weights_16(14) * f(xmid + h * gl_nodes_16(14)) & | ||
| 285 | & + gl_weights_16(15) * f(xmid + h * gl_nodes_16(15)) & | ||
| 286 | 1 | & + gl_weights_16(16) * f(xmid + h * gl_nodes_16(16)) | |
| 287 | 1 | ans = ans * h | |
| 288 | 1 | end procedure | |
| 289 | |||
| 290 | end submodule | ||
| 291 |