Top new questions this week:
|
|
I need to compute an eigendecomposition of an $n\times n$ matrix $$ D + c vv^\top = Q\Lambda Q^\top \tag{1} $$ in MATLAB, where $D$ is a real diagonal matrix, $c$ is a scalar, and $v$ is a real vector….
|
|
I have a function $f : \mathbb{R}^2 \mapsto \mathbb{R}^+$ and I wish to numerically evaluate the integral below over a finite domain $\Omega \subset \mathbb{R}^2$ $$ I = \int_\Omega e^{i k \cdot f(\…
|
|
I am currently using the protocol described in stackoverflow.com/questions/55599950/computation-of-pseidoinverse-with-svd-in-c-using-blas-and-lapacke to compute the pseudo inverse of a matrix. …
|
|
Notation: Denote $T_{h} = \left\{K\right\}$ to be a face-conforming triangulation of a domain $\Omega$ such that $K_{i} \cap K_{j} = \emptyset$ for $i \neq j.$ Additionally, denote $\mathcal{V}_{h} = …
|
|
Consider the advection/transport equation in 1D with constant velocity $a(x) \equiv 1$ $$u_t(t,x) + u_x(t,x) = 0$$ on a, say, periodic domain. On uniform grids $$ \{x_i\}_{i = 1, \dots, N}, \quad x_{i …
|
Greatest hits from previous weeks:
|
|
While reading some literature on PDE solvers I came across the term pseudo time-stepping today. It seems to be a common term, however I failed to find a good definition or an introductionary article …
|
|
Is there a Fortran subroutine which performs linear interpolation in one-dimenional data? I need something similar to MATLAB function interp1.
|
|
I was comparing a few of my codes to “stock” MATLAB codes. I am surprised at the results. I ran a sample code (Sparse Matrix) …
|
|
Which software provides a good workflow from simple plotting of a few datapoints up to the creation of publication level graphics with detailed styles, mathematical typesetting and “professional …
|
|
I am using image filtering for an image processing algorithm I’m developing. I’m using a predefined Matlab function to do the convolution, but I’d like to know what the computational complexity is for …
|
|
I am getting the warning in the post subject when attempting to optimize a function in Python with the scipy.optimize.fmin_bfgs function. The complete output: Warning: Desired error not necessarily …
|
|
The Fast Fourier Transform algorithm computes a Fourier decomposition under the assumption that its input points are equally spaced in the time domain, $t_k = kT$. What if they’re not? Is there …
|
