New answers tagged

4 votes
Accepted

How do I evaluate the numerical solution at a point that isn't on the mesh

I'll give a more general answer independent of implementation. TL:DR The point $x$ you seek the value $u$ at lies in a simplex with values at the nodes $u_i$. $x$ has the barycentric coordinates $\...
Sardine's user avatar
  • 306
1 vote

How do I evaluate the numerical solution at a point that isn't on the mesh

After the suggestion from @whpowell96 I used interpolation like this ...
tareqath's user avatar
2 votes
Accepted

Possible bug with scipy.optimize SHGO sobol: TypeError: <lambda>() takes 1 positional argument but 3 were given

As also done here I report a solution for the sake of completeness. The issue discussed is an unsolved (yet) bug of shgo and it's related to how the algorithm feeds ...
Sasche's user avatar
  • 31
1 vote

How to program the convergence of a sequence of systems of integral equations using Scipy

You are computing convolutions. These are best evaluated as convolutions, rather than as sequences of integrals. A simple way of doing this is to discretize and then end up with a matrix-vector ...
Wolfgang Bangerth's user avatar
1 vote
Accepted

How to program the convergence of a sequence of systems of integral equations using Scipy

What you are doing is essentially a fixed-point iteration. That is, you are computing $\lim_{n\to\infty} U_n$, where $U_{n+1} = F(U_n)$. Here $U_n = (u_n,v_n)^T$ and $F$ is a vector containing your ...
whpowell96's user avatar
  • 1,532

Top 50 recent answers are included