19"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
21+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
23-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
25call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
26err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
27if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
28getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
290.592319847946765693983261139952398973,
0.592319847946765693983261139952399069,
0.420679694161706175066609104098715408E-29,
0.962964972193617926527988971292463659E-34 >= 1 (unbiased)? TRUE
30[numFuncEval, numInterval]
34"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
36+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
38-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
40call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
41err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
42if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
43getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
440.592319847946765693983261139952398973,
0.592319847946765693983261139952399069,
0.420679694161706175066609104098715408E-29,
0.962964972193617926527988971292463659E-34 >= 1 (unbiased)? TRUE
45[numFuncEval, numInterval]
53"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
55+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
57-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
59call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
60err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
61if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
62getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
630.592319847946765693983261139952398973,
0.592319847946765693983261139952395699,
0.749057713854147205376412430387292458E-25,
0.327408090545830095019516250239437644E-32 >= 1 (unbiased)? TRUE
64[numFuncEval, numInterval]
68"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
70+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
72-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
74call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
75err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
76if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
77getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
780.592319847946765693983261139952398973,
0.592319847946765693983261139952395699,
0.749057713854147205376412430387292458E-25,
0.327408090545830095019516250239437644E-32 >= 1 (unbiased)? TRUE
79[numFuncEval, numInterval]
87"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
89+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
91-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
93call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
94err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
95if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
96getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
970.592319847946765693983261139952398973,
0.592319847946765693983261139952399069,
0.282354400096319256248938688151661529E-30,
0.962964972193617926527988971292463659E-34 >= 1 (unbiased)? TRUE
98[numFuncEval, numInterval]
102"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
104+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
106-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
108call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
109err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
110if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
111getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
1120.592319847946765693983261139952398973,
0.592319847946765693983261139952399069,
0.282354400096319256248938688151661529E-30,
0.962964972193617926527988971292463659E-34 >= 1 (unbiased)? TRUE
113[numFuncEval, numInterval]
121"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
123+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
125-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
127call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
128err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
129if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
130getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
1310.592319847946765693983261139952398973,
0.592319847946765693983261139952369891,
0.130763857408511882147171505007517289E-23,
0.290815421602472613811452669330324025E-31 >= 1 (unbiased)? TRUE
132[numFuncEval, numInterval]
136"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
138+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
140-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
142call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
143err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
144if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
145getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
1460.592319847946765693983261139952398973,
0.592319847946765693983261139952369891,
0.130763857408511882147171505007517289E-23,
0.290815421602472613811452669330324025E-31 >= 1 (unbiased)? TRUE
147[numFuncEval, numInterval]
155"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
157+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
159-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
161call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
162err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK31, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
163if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
164getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
1650.592319847946765693983261139952398973,
0.592319847946765693983261139952398973,
0.908014202174603824005803987572878882E-26,
0.00000000000000000000000000000000000 >= 1 (unbiased)? TRUE
166[numFuncEval, numInterval]
170"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
172+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
174-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
176call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
177err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK31,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
178if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
179getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
1800.592319847946765693983261139952398973,
0.592319847946765693983261139952398973,
0.908014202174603824005803987572878882E-26,
0.00000000000000000000000000000000000 >= 1 (unbiased)? TRUE
181[numFuncEval, numInterval]
189"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
191+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
193-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
195call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
196err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK21, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
197if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
198getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
1990.592319847946765693983261139952398973,
0.592319847946765693983261139950914562,
0.242973532478881507761272478019151270E-22,
0.148441050463646203374289499924733273E-29 >= 1 (unbiased)? TRUE
200[numFuncEval, numInterval]
204"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
206+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
208-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
210call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
211err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK21,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
212if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
213getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
2140.592319847946765693983261139952398973,
0.592319847946765693983261139950914466,
0.242973532478881507761272477921271703E-22,
0.148450680113368139553554779814446198E-29 >= 1 (unbiased)? TRUE
215[numFuncEval, numInterval]
223"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
225+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
227-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
229call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
230err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK15, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
231if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
232getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
2330.592319847946765693983261139952398973,
0.592319847946765693983261139952415054,
0.311644046752404524231338211060753717E-22,
0.160815150356334193730174158205841431E-31 >= 1 (unbiased)? TRUE
234[numFuncEval, numInterval]
238"int1_type: an algebraic integrand of the form f(x) = x**2 / (x**2 + 1) / (x**2 + 4) for x in (lb, ub)"
240+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
242-4.00000000000000000000000000000000000,
+4.00000000000000000000000000000000000
244call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
245err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK15,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
246if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
247getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
2480.592319847946765693983261139952398973,
0.592319847946765693983261139952415054,
0.311644046752404524231537458497527643E-22,
0.160815150356334193730174158205841431E-31 >= 1 (unbiased)? TRUE
249[numFuncEval, numInterval]
262"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
264+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
266+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
268call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
269err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
270if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
272 **** INTEGRATION FAILED: ERR
= 3 ****
274getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
2750.942809041582063365867792482806465323,
0.942809041582063353517270922122486063,
0.275996032154812596404096004412307632E-15,
0.123505215606839792599902659423179129E-16 >= 1 (unbiased)? TRUE
276[numFuncEval, numInterval]
280"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
282+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
284+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
286call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
287err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
288if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
289getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
2900.942809041582063365867792482806465323,
0.942809041582063365867792482806360745,
0.780290516968488605865629463438283303E-30,
0.104577995980226906820939602282361553E-30 >= 1 (unbiased)? TRUE
291[numFuncEval, numInterval]
299"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
301+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
303+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
305call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
306err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
307if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
309 **** INTEGRATION FAILED: ERR
= 3 ****
311getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
3120.942809041582063365867792482806465323,
0.942809041582063354181389418572630975,
0.395526291759776583849399257993266119E-15,
0.116864030642338343481331615363661636E-16 >= 1 (unbiased)? TRUE
313[numFuncEval, numInterval]
317"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
319+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
321+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
323call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
324err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
325if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
326getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
3270.942809041582063365867792482806465323,
0.942809041582063365867792482806368256,
0.683127351274152557078955376234873720E-30,
0.970668691971166869940212883062803369E-31 >= 1 (unbiased)? TRUE
328[numFuncEval, numInterval]
336"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
338+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
340+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
342call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
343err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
344if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
346 **** INTEGRATION FAILED: ERR
= 3 ****
348getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
3490.942809041582063365867792482806465323,
0.942809041582063353612391313951514454,
0.393612625583696433467216358595637327E-15,
0.122554011688549508692539077120199527E-16 >= 1 (unbiased)? TRUE
350[numFuncEval, numInterval]
354"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
356+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
358+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
360call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
361err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
362if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
363getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
3640.942809041582063365867792482806465323,
0.942809041582063365867792482806399167,
0.654527291600002104661074103787487549E-30,
0.661556935897015515524728423277922534E-31 >= 1 (unbiased)? TRUE
365[numFuncEval, numInterval]
373"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
375+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
377+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
379call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
380err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
381if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
383 **** INTEGRATION FAILED: ERR
= 3 ****
385getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
3860.942809041582063365867792482806465323,
0.942809041582063354504151452306557739,
0.393427741808206315477620382019741615E-15,
0.113636410304999075837096309795641304E-16 >= 1 (unbiased)? TRUE
387[numFuncEval, numInterval]
391"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
393+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
395+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
397call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
398err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
399if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
400getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
4010.942809041582063365867792482806465323,
0.942809041582063365867792482806474856,
0.311326575510196675646498834418853501E-30,
0.953335322471681747262709081579539023E-32 >= 1 (unbiased)? TRUE
402[numFuncEval, numInterval]
410"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
412+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
414+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
416call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
417err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK31, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
418if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
420 **** INTEGRATION FAILED: ERR
= 3 ****
422getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
4230.942809041582063365867792482806465323,
0.942809041582063360076317716910846285,
0.279164711945302751623913154970201256E-15,
0.579147476589561903767056045130549251E-17 >= 1 (unbiased)? TRUE
424[numFuncEval, numInterval]
428"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
430+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
432+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
434call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
435err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK31,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
436if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
437getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
4380.942809041582063365867792482806465323,
0.942809041582063365867792482806439901,
0.205400428568898703728420047576682499E-30,
0.254222752659115132603389088421210406E-31 >= 1 (unbiased)? TRUE
439[numFuncEval, numInterval]
447"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
449+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
451+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
453call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
454err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK21, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
455if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
457 **** INTEGRATION FAILED: ERR
= 3 ****
459getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
4600.942809041582063365867792482806465323,
0.942809041582063358489159110539769991,
0.274189625236721766830515122145815584E-15,
0.737863337226669533225571680774701636E-17 >= 1 (unbiased)? TRUE
461[numFuncEval, numInterval]
465"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
467+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
469+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
471call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
472err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK21,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
473if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
474getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
4750.942809041582063365867792482806465323,
0.942809041582063365867792482806451456,
0.747866619686254516023298199584791888E-27,
0.138666955995880981420030411866114767E-31 >= 1 (unbiased)? TRUE
476[numFuncEval, numInterval]
484"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
486+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
488+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
490call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
491err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK15, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
492if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
494 **** INTEGRATION FAILED: ERR
= 3 ****
496getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
4970.942809041582063365867792482806465323,
0.942809041582063352985651916754259394,
0.107881013975510488055028051825168280E-14,
0.128821405660522059292977722684510948E-16 >= 1 (unbiased)? TRUE
498[numFuncEval, numInterval]
502"int2_type: an algebraic integrand of the form f(x) = 1 / sqrt(a - b * x) for x in (0, a / b) with a > 0 and b > 0 with a singularity at the upper bound of integration"
504+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
506+0.00000000000000000000000000000000000,
+0.666666666666666666666666666666666635
508call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
509err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK15,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
510if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
511getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
5120.942809041582063365867792482806465323,
0.942809041582063365867792482806467056,
0.894595922874628788043750076873935104E-26,
0.173333694994851226775038014832643459E-32 >= 1 (unbiased)? TRUE
513[numFuncEval, numInterval]
526"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
528+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
530+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
532call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
533err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
534if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
535getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
536-3.12249862669012531099145873075559453,
-3.12249862669012531099145703288302519,
0.159554531774406733373916406327365902E-21,
0.169787256934451856030425956397998728E-23 >= 1 (unbiased)? TRUE
537[numFuncEval, numInterval]
541"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
543+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
545+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
547call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
548err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
549if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
550getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
551-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075586917,
0.104192809991349459650328406693844568E-29,
0.274637610069619832645782454612610636E-30 >= 1 (unbiased)? TRUE
552[numFuncEval, numInterval]
560"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
562+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
564+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
566call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
567err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
568if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
569getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
570-3.12249862669012531099145873075559453,
-3.12249862669012531099145675992893943,
0.158821823031096796049113138417448082E-21,
0.197082665510145750102784224886897408E-23 >= 1 (unbiased)? TRUE
571[numFuncEval, numInterval]
575"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
577+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
579+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
581call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
582err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
583if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
584getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
585-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075560878,
0.295052467480124532688175820804010865E-30,
0.142518815884655453126142367751284622E-31 >= 1 (unbiased)? TRUE
586[numFuncEval, numInterval]
594"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
596+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
598+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
600call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
601err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
602if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
603getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
604-3.12249862669012531099145873075559453,
-3.12249862669012531099145637935769671,
0.157863935952385241531199129294364333E-21,
0.235139789781688314990001714829882296E-23 >= 1 (unbiased)? TRUE
605[numFuncEval, numInterval]
609"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
611+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
613+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
615call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
616err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
617if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
618getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
619-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075562303,
0.278104283969516857181283214909263505E-30,
0.285037631769310906252284735502569243E-31 >= 1 (unbiased)? TRUE
620[numFuncEval, numInterval]
628"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
630+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
632+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
634call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
635err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
636if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
637getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
638-3.12249862669012531099145873075559453,
-3.12249862669012531099145581781496689,
0.156352167288494027740417861647194691E-21,
0.291294062764122468285662650331711561E-23 >= 1 (unbiased)? TRUE
639[numFuncEval, numInterval]
643"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
645+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
647+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
649call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
650err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
651if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
652getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
653-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075552443,
0.445660189131206376397153295914152182E-30,
0.701038499756953850512375971100913544E-31 >= 1 (unbiased)? TRUE
654[numFuncEval, numInterval]
662"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
664+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
666+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
668call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
669err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK31, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
670if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
671getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
672-3.12249862669012531099145873075559453,
-3.12249862669012531099145489154805592,
0.153862855021951093878837768152750949E-21,
0.383920753861232441434641228203070917E-23 >= 1 (unbiased)? TRUE
673[numFuncEval, numInterval]
677"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
679+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
681+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
683call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
684err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK31,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
685if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
686getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
687-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075556603,
0.182963344716787406040317904545568095E-30,
0.285037631769310906252284735502569243E-31 >= 1 (unbiased)? TRUE
688[numFuncEval, numInterval]
696"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
698+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
700+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
702call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
703err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK21, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
704if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
705getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
706-3.12249862669012531099145873075559453,
-3.12249862669012531099145729306953740,
0.166762278786010609240106982579981836E-21,
0.143768605713088174426049560344490814E-23 >= 1 (unbiased)? TRUE
707[numFuncEval, numInterval]
711"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
713+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
715+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
717call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
718err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK21,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
719if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
720getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
721-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075553329,
0.638742935999207969456768107027576605E-25,
0.612445722315141001271800985742006887E-31 >= 1 (unbiased)? TRUE
722[numFuncEval, numInterval]
730"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
732+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
734+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
736call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
737err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK15, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
738if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
739getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
740-3.12249862669012531099145873075559453,
-3.12249862669012531099145821665001954,
0.168984966934533206681021840309021027E-21,
0.514105574994657236239111087685608621E-24 >= 1 (unbiased)? TRUE
741[numFuncEval, numInterval]
745"int3_type: an algebraic integrand of the form f(x) = log(x) / sqrt(x) for x in (0, ub) with a singularity at the lower limit of integration"
747+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
749+0.00000000000000000000000000000000000,
+3.00000000000000000000000000000000000
751call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
752err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK15,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
753if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
754getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
755-3.12249862669012531099145873075559453,
-3.12249862669012531099145873075569198,
0.736077254850655834241784402267592237E-24,
0.974520551859941341646324838947973223E-31 >= 1 (unbiased)? TRUE
756[numFuncEval, numInterval]
769"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
771+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
773+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
775call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
776err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
777if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
778getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
779-0.189275187882093321180367135892330336,
-0.189275187882093321180367157318588076,
0.709016787220663753280987987987449705E-23,
0.214262577393144583924325489311691073E-25 >= 1 (unbiased)? TRUE
780[numFuncEval, numInterval]
784"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
786+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
788+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
790call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
791err = getQuadErr(getFunc, lb, ub, abstol, reltol, nodeK, weightK, weightG,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
792if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
793getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
794-0.189275187882093321180367135892330336,
-0.189275187882093321180367124314328969,
0.875190015574302935869778694464445775E-23,
0.115780013668100956229775935139168619E-25 >= 1 (unbiased)? TRUE
795[numFuncEval, numInterval]
803"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
805+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
807+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
809call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
810err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
811if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
812getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
813-0.189275187882093321180367135892330336,
-0.189275187882093321180367149919369956,
0.551012606285227548119424458552174544E-23,
0.140270396193317866081626615262131133E-25 >= 1 (unbiased)? TRUE
814[numFuncEval, numInterval]
818"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
820+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
822+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
824call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
825err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK61,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
826if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
827getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
828-0.189275187882093321180367135892330336,
-0.189275187882093321180367089569735215,
0.151044926215718500490728968093616592E-23,
0.463225951208172227937270578607855780E-25 >= 1 (unbiased)? TRUE
829[numFuncEval, numInterval]
837"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
839+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
841+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
843call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
844err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
845if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
846getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
847-0.189275187882093321180367135892330336,
-0.189275187882093321180367156609609617,
0.990374035889893040148759005915834788E-23,
0.207172792803981521991901235694034717E-25 >= 1 (unbiased)? TRUE
848[numFuncEval, numInterval]
852"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
854+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
856+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
858call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
859err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK51,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
860if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
861getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
862-0.189275187882093321180367135892330336,
-0.189275187882093321180367136737390768,
0.206538524316503455750529471463185887E-23,
0.845060431740830626700484040037212854E-27 >= 1 (unbiased)? TRUE
863[numFuncEval, numInterval]
871"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
873+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
875+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
877call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
878err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
879if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
880getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
881-0.189275187882093321180367135892330336,
-0.189275187882093321180367151548235028,
0.957392833849128474709628758212286240E-23,
0.156559046920400077057876054566418227E-25 >= 1 (unbiased)? TRUE
882[numFuncEval, numInterval]
886"int4_type: an algebraic integrand of the form f(x) = log(x) / (1. + log(x)**2)**2 for x in (0, 1) with a singularity at the lower limit"
888+0.00000000000000000000000000000000000,
+0.559579552072376591733071579965869345E-22
890+0.00000000000000000000000000000000000,
+1.00000000000000000000000000000000000
892call setNodeWeightGK(order
= 35_IK, nodeK
= nodeK, weightK
= weightK, weightG
= weightG)
893err = getQuadErr(getFunc, lb, ub, abstol, reltol,
GK41,
weps, integral, abserr, sinfo, sindex, numFuncEval, numInterval)
894if (err
/= 0_IK)
call disp%show(getStrUpper(SK_
' **** integration failed: err = '//getStr(err)
//SK_
' ****'), tmsize
= 1_IK, bmsize
= 1_IK)
895getStr([truth, integral, abserr, abs(integral
- truth)])
//SK_
' >= 1 (unbiased)? '//getStr(abs(integral
- truth)
<= abserr)
896-0.1892751878820933211803671358