Questions tagged [scipy]
SciPy is an open source library of algorithms and mathematical tools for the Python programming language.
87 questions
2 votes
2 answers
107 views
Heat Diffusion Equation Solver
I'm implementing a numerical solver for the 1D diffusion (heat) equation using an iterative method. The initial temperature distribution is a Gaussian profile, which should theoretically smooth out ...
6 votes
1 answer
460 views
Implementing a joint differential equation and eigenvalue solver
I've written a Python script to solve a system of coupled differential equations and then use the results to solve an eigenvalue problem. The full problem is: \$\frac{d^2}{dr^2} K(r) = \frac{K(r)}{r^2}...
9 votes
2 answers
515 views
Correctly sum pixel values into bins of angle relative to center
A quick way to try to do this would be to use scipy.ndimage.map_coordinates, mapping radius to one axis of a rectangular array, and angle to the other, then summing along the radius axis. But, I think ...
- 605
5 votes
2 answers
319 views
Srivastava multivariate Fox H function in MATLAB
I was trying to rewrite the Python code in MATLAB. The result is consistent. But, the MATLAB code is so slow. Any help would be appreciated. Ref python code link The MATLAB code written by me is as ...
- 469
2 votes
1 answer
164 views
Multithreaded sparse linear solver in Python
I would like to know if there are any multithreaded sparse linear solvers for LU decomposed sparse matrices. I want to increase the solving speed of a system like: ...
- 29
3 votes
1 answer
129 views
"Flattening" a 2D STM image by subtracting an Nth order polynomial without typing out all the terms for each N?
I have some scanning tunneling microscopy (STM) topography images that I need to flatten. I have found that for some data I need a 4th order polynomial, but in other cases a lower order works better. ...
- 605
4 votes
2 answers
221 views
Matrix Factorisation class packaging methods for factorisation of explicit & implicit data matrices using Gradient Descent, SGD and ALS
Attached below, and also as this GitHub gist, is code for a Python class I wrote as part of a personal learning/portfolio project on collaborative-filtering recommender systems via matrix ...
- 43
1 vote
2 answers
165 views
Python 3.10+ deconstructing an over engineered solution to better understand how metaclasses work with properties, static methods, and classmethods
TL;DR This question examines an over-engineered example of python metaclasses and dataclasses to create a LiteralEnum (for validating a stringly-typed keyword ...
- 291
2 votes
1 answer
128 views
Local solver built into a global one
I am trying to optimize writing a script that has 2 minimizations, one dependent on the other. My code is a bit bloated, and I find the global parameters that I am solving for depend quite strongly on ...
- 491
8 votes
2 answers
557 views
Fourier Series of a given function
This is a very simple code that expresses a function in terms of Trigonometric Fourier Series and generates a animation based on the harmonics. I would like to know some ways to improve the ...
- 83
1 vote
1 answer
94 views
Fixing math library functions in Black-Scholes options pricing model
I've amended a code for the Black-Scholes formula for European pricing options found here at the bottom of the page and fixed the math functions accordingly. Unlike the code on the website, mine has ...
- 113
9 votes
1 answer
4k views
2 votes
2 answers
2k views
Solving the TDoA multilateration problem in 3-dimensions
Background. I've written an algorithm to solve the Time Difference of Arrival (TDoA) multilateration problem in 3-dimensions. That is, given the known coordinates of ...
- 305
2 votes
1 answer
277 views
Integration loop over multiple doping and temperature levels
I want to perform some calculations on a large dataset. The code can be found below, where I want to calculate the values for 'results_nr' over a large loop (1000 x 910) values. Can you help me out ...
- 21
1 vote
1 answer
140 views
Optimization using a differential evolution technique, ideally with a million iterations
I'm working on an optimization code using differential evolution to solve. However, it is taking a long time to get the solution. See that I have a variable called number_of_iterations that must equal ...
- 11