New answers tagged scipy
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 $\...
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
...
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 ...
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 ...
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 ...
Top 50 recent answers are included
Related Tags
scipy × 247python × 152
ode × 44
numpy × 39
optimization × 28
numerics × 22
integration × 22
differential-equations × 19
sparse-matrix × 18
linear-algebra × 17
quadrature × 12
matrix × 11
interpolation × 11
nonlinear-equations × 10
constrained-optimization × 10
eigenvalues × 8
curve-fitting × 8
computational-physics × 7
pde × 6
numerical-modelling × 6
lapack × 6
least-squares × 6
matrix-factorization × 6
rootfinding × 6
linear-solver × 5