Top new questions this week:
|
|
I just for the first time saw the function softmax() in this SO answer to How do I use a minimization function in scipy with constraints and was intrigued. Another way of weighting variables where …
|
|
I am trying to use the Kramers-Kronig algorithm to transform the real and imaginary contributions to the anomalous scattering factor from a diffraction anomalous fine structure (DAFS) experiment. I …
|
|
Given the discrete Lyapunov equation $$AXA^T – X + Q = 0$$ how can I solve for $X$ as a function of the eigenvectors of some matrix $H$? More precisely, in the case of the continuous Lyapunov equation …
|
|
I am trying to solve a coupled system of ODE’s using the solve_ivp function from scipy. The general form of the equation is given via $$\dot{y}(t) = M(t)y(t).$$ The time dependence of matrix is …
|
|
I am looking up celestial navigation, and according to https://youtu.be/-ARXW8InStY?t=3320 a specific sun angle reading (sun angle above the horizon) will be the same on a small-circle with the centre …
|
|
I have experience with classical linear algebra packages in C++ like Eigen, Blaze etc. I have never wrote my own PetSC back-end solver but I used/modified several of them. In a new project I would …
|
|
Given the Algebratic Riccati Equation (ARE) $$A^T X + XA + XRX + Q = 0$$ where $A,R,Q \in \mathbb R^{n \times n}$, we are interested in the matrix $X$ that solves this equation. If we define the $2n \…
|
Greatest hits from previous weeks:
|
|
Going to teach students of undergraduate level a course titled Introduction to Computer Programming. I am confused a bit. In Computational Physics scientists use C/C++ or Python or Fortran,CUDA etc…….
|
|
Is there a Fortran subroutine which performs linear interpolation in one-dimenional data? I need something similar to MATLAB function interp1.
|
|
What is the preferred and efficient approach for interpolating multidimensional data? Things I’m worried about: performance and memory for construction, single/batch evaluation handling dimensions …
|
|
I know that the Nitsche’s method is a very attractive methods since it allows to take into account Dirichlet type boundary conditions or contact with friction boundary conditions in a weak way without …
|
|
Is there a way, using some established Python package (e.g. SciPy) to define my own probability density function (without any prior data, just $f(x) = a x + b$), so I can then make calculations with …
|
|
Overview My understanding is that one should use a time step $\Delta t < \frac{h}{v}$ (where h – smallest mesh element, v – velocity) to get an accurate result. But how important is this really …
|
|
I want to plot a numerical integral function of some function $f$ using scipy and matplotlib. How can I do this? I tried the …
|
