- finished a serial implementation of extrapolation based on Euler’s method with a convergence test. However the convergence test still shows that the order is 1 for any value of p. I don’t know why yet.
- ideas to improve this basic serial implementation: 1) use memoization technique for the evaluations of
f
. 2) make the computation of the first loop in ex_one_step()
parallel by assigning each iteration of the loop to a thread. The computation of each iteration is independent of the other, so synchronization between threads is easy.