2,410 reputation
720
bio website
location Berkeley, CA
age
visits member for 2 years, 3 months
seen yesterday

Mar
18
answered NIntegrate(with singularities) gives different results, which one should I trust?
Feb
12
awarded  Informed
Jan
28
answered NIntegrate::slwcon Problem
Jan
28
comment NIntegrate::slwcon Problem
Unfortunately, this is actually an example of LocalAdaptive failing: It incorrectly believes it has converged to a precision of 11. But the correct answer is zero, so any finite answer (such as -1.77636*10^-14) is incorrect to any precision.
Jan
17
awarded  Yearling
Jan
7
awarded  Good Answer
Oct
25
comment How to collect result continuously (interruptible calculation) when running parallel calculations?
@xzczd The Reap & ParallelSow are used to preserve partially computed results if the computation is aborted before finishing.
Aug
25
comment why maxrecursion didn't work if I specify singularity in NIntegrate
Isn't the exact value of the integral 1? With MaxRecursion -> 5, and splitting the integration region first at x == 2, NIntegrate converges to 1. Thus 5 recursions are sufficient. And any MaxRecursion -> n for n > 5 should be expected to be identical.
Aug
25
comment why maxrecursion didn't work if I specify singularity in NIntegrate
What are you expecting to happen? What you are seeing is that 5 recursions (within each of the sub-intervals {-Infinity, 2} and {2, Infinity}) are sufficient to solve the integral to within PrecisionGoal -> Automatic.
Aug
6
awarded  Nice Answer
Aug
5
revised Elegant way of obtaining the envelope of oscillating function
A little explanation about how it works
Jul
26
answered NIntegrate extremely piecewised functions
Jul
9
revised Speed of convergence for NIntegrate
Clarify slwcon criterion
Jul
9
answered Speed of convergence for NIntegrate
Jul
9
comment Speed of convergence for NIntegrate
WorkingPrecision -> 10 amounts to a decrease of working precision from ~16 (the value of $MachinePrecision). It also decreases the value of PrecisionGoal -> Automatic, to ~4 I think (from a default of ~6).
Jun
29
answered Elegant way of obtaining the envelope of oscillating function
Jun
14
comment Exporting multiple arrays in a single HDF5 file in Mathematica 7
Can you give a specific example of the values of the arrays data, x1, and y1? I tried the last example from the HDF5 documentation in Mathematica 9 and it seemed to work correctly exporting two matrices: Export["test.h5", {{{1, 2}, {2, 3}}, {{10, 11}, {12, 13}}}, {"Datasets", {"m1", "m2"}}]
Jun
13
answered Determining which rule NIntegrate selects automatically
Jun
13
answered NIntegrate converging too slowly when increasing size of array
Jun
11
awarded  Citizen Patrol