All Questions

1,949 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
4
votes
0answers
452 views

Finite volume method

I have question connected with finite volume method. Consider equation $$\frac{\partial u}{\partial t}=\operatorname{div}A\nabla u +f . \quad x\in \overline{B}_{1} (0)\subset \mathbb{R}^3 -\text{unit ...
4
votes
0answers
494 views

full rank update to cholesky decomposition for multivariate normal distribution

This question is a specialization of full rank update to cholesky decomposition, to which I hope to get a more positive answer. When calculating the minus log of the multivariate normal distribution, ...
4
votes
0answers
86 views

How do I perform chebyshev interpolation from a to b with custom angle range?

Typically Chebyshev interpolation from $-1$ to $1$ with angle from $0$ to $\pi$: $\xi_j=\cos \left ({\pi j \over N}\right )$ $x_j=(1+\xi_j) * {L \over 2}$ $w$: $w_0=\pi/(2N)$ $w_{1,...,N-1}=\pi/(N)$...
4
votes
0answers
82 views

Computing linear combinations of sines and cosines (phasors)

I have a finite series that looks like this: $f(t) = \sum^n_{i=0} A_i cos(\Theta_i + \omega_i t) + B_i sin(\Theta_i + \omega_i t)$ That is, a finite series of pairs of phasors. What's the state of ...
4
votes
0answers
235 views

Sound Waves Simulation in 3D Environment

I want to do a simulation of sound waves including wave propagation, absorption, and reflection in 3D space. I did some research and I found this question in stackoverflow but it talks about ...
4
votes
0answers
654 views

Finding roots of systems of equations with a Jacobian that is singular everywhere

Given $a\in\mathbb{R}^{mn\times n}$, find a $C\in\mathbb{R}^{n}$, $x\in\mathbb{R}^{m\times n}$ such that $$ 0 = f_{k}(\boldsymbol{C}, \boldsymbol{x}):=\sum_{i=1}^{m} C_{i} \left(\prod_{j=1}^{n} a_{kj}^...
4
votes
0answers
232 views

Negative viscosity stabilized by fourth order terms

I am trying to solve a "Navier-Stokes"-type problem where the viscosity is negative. Of course this renders the equation unstable and thus I add a fourth order term, so the entire equation becomes: $$...
4
votes
0answers
113 views

What are the governing equations solved in coupled atmosphere-ocean models?

In (hydrostatic) atmospheric general circulation models, for example the so-calle Primitive Equations, consisting of the horizontal momentum equation, the hydrostatic balance, the continuity equation, ...
4
votes
0answers
184 views

Time-stepping for coupled nonlinear PDEs

What are good references for time-stepping of the coupled incompressible Navier-Stokes-heat equation (Boussinesq flow), $$ \begin{cases} \rho\left(\dot{\mathbf{u}} + \mathbf{u}\cdot\nabla \mathbf{u}\...
4
votes
0answers
333 views

Understanding the meaning of Computational Order of Convergence

I am a postgraduate student with interest in numerical methods for solving nonlinear systems of equations. I have read some papers that discussed about 'computational order of convergence' for some ...
4
votes
0answers
1k views

Using MINPACK for curve fitting: implementation?

I need to implement a non-linear fitting algorithm in Fortran and chose to use MINPACK's flavor of the Levenberg-Marquardt algorithm as a basis for the least-squares stuff. However, I seem to ...
4
votes
0answers
77 views

Order of convergence of Scrodinger eq. with CN scheme

I'm trying to solve numerically the 1-dim time dependent Schrodinger equation using the Crank Nicolson scheme and the Thomas algorithm to solve the tridiagonal matrix. The physical system consists of ...
4
votes
0answers
683 views

Generating pseudo-random orthonormal bases for random projection

I am performing series of random projections i.e. projecting the input matrix onto randomly generated orthonormal bases (of much lower dimensionality). The projection is just a matrix multiplication ...
4
votes
0answers
83 views

Non-convergance when calculating temperature/heat flows through a section of rock

I am attempting to calculate temperature of section of rock in the earth as a function of vertical position in the rock and time. Along with it I am calculating the heat flow through the rock as a ...
4
votes
0answers
505 views

Fenics: Result of Steady state dynamic linear elastic doesn't match with actual values

I solved the steady state dynamic linear elastic model in a solid. My equation is a function of frequency and the strong form is: $$\operatorname{div}(\operatorname{stress}(\vec x, w)) + w^2 \rho u( \...
4
votes
0answers
223 views

Adaptive mesh data structure for Fast Marching Method to overcome RAM limit

On an uniform mesh of positions in space $\ (x_i,y_j,z_k)$: $$\ x_i = x_0 + i\Delta x,\quad i=0,\ldots,n_x$$ $$\ y_j = y_0 + j\Delta y,\quad j=0,\ldots,n_y$$ $$\ z_k = z_0 + k\Delta z,\quad k=0,\...
4
votes
0answers
170 views

Solving diffusion PDE using finite differences

I need some hints on how to solve this diffusion equation ($\alpha, k_1,k_2$ and $k_3$ are constants): $$ {\partial P \over \partial y} + k_1 {\partial P \over \partial t} + \alpha P = {1 \over k_2} ...
4
votes
0answers
174 views

Find maximum distance between elements given constraints on some

I have a list of numbered elements 1 to N that fit into positions on a number line starting with 1. I also have constraints for these elements: The element 1 is in position 1, and element N must be ...
4
votes
0answers
106 views

Existing software/scripts for spiral graphs?

I am looking for existing software or scripts to generate spiral graphs from cyclical (time) data, as presented by Webber and Muller. The graphs shown in the paper look like a great means of ...
4
votes
0answers
657 views

Perron-Frobenius theorem on general real symmetric matrices

From the Perron-Frobenius theorem, it might be concluded that the spectral radius is the largest eigenvalue for positive matrices, ie, matrices with strictly positive entries. In other words, the ...
4
votes
0answers
99 views

Time-stable spectral decomposition algorithm

Consider an $n \times n$ real, time-dependent matrix $A(t)$ such that $A(t) = A(t)^{T} > 0$ and $A(t)$ is continuous on $[a,b]$. Then it is posible to specify a matrix $S(t) \in SO(n)$ such that $S(...
4
votes
0answers
413 views

Why is my lower convex hull extraction algorithm not working?

Recently, I wrote an algorithm to obtain a delaunay triangulation of a random point set in $I=[-10,10]$x$[-10,10] \subset R^2$ by projecting these points onto the 3 dimensional paraboloid $z=x^2+y^2$, ...
4
votes
0answers
134 views

Probabilistic algorithms for matrix approximation

Considering regular matrix approximation inequality || $A - QQ^TA $|| < e where we try to approximate matrix $A$ by a lower rank orthonormal matrix $Q$. I've read an article on probabilistic ...
4
votes
0answers
254 views

Convergence rate of Monte-Carlo variance estimate

What is the convergence rate for Monte-Carlo variance estimate for a random variable $X \in {L^q}(\Omega ,R),2 < q < 4$?
3
votes
0answers
23 views

Numerical Issues with DDE from SEIRU$\delta$ model

I'm new in this community. I moved this question from Math community. I'm reading the following article Article Here and my target is to replicate the results for a project. SEIRU Model: I obtained ...
3
votes
1answer
84 views

Solve two-player game - minimize the l-infinity norm of a matrix-vector product

I have a matrix $M$ with non-negative real entries, and I would like to minimize the objective function $$\Phi(v) = \|Mv\|_\infty,$$ where $v$ is constrained to be a probability vector, i.e., $v_1+\...
3
votes
0answers
51 views

Choosing good modelling method for solving Boltzmann equation

I'm writing a solver for Boltzmann Equations (BE) including a force term in rarefied plasma, for my PhD. The aim is to see if an instability occurs inside an electric streamer (theoretically it should,...
3
votes
1answer
107 views

FEM with elastic inhomogeneous properties leads to mesh-induced anisotropy

I'm solving an elastic homogenization problem and I'm having problems with mesh artifacts. I would like to first give a brief summary of what I do: I have a system with inhomogeneous (but isotropic) ...
3
votes
0answers
49 views

Some formulations of domain coupling lead to saddle point problems. Is this merely an artifact?

Background I wanted to learn how to couple FEM and BEM (for the Poisson equation), because I wanted to better understand how open boundary conditions look like. Therefore I worked through the relevant ...
3
votes
0answers
86 views

Compute Nullspace of Sparse Matrix

I am computing the nullspace of a sparse rectangular $m$ x $n$ matrix $A$, where $m$ << $n$. I do this by computing the QR decomposition of $A^T$ and extract the $n-m$ right-most columns of the ...
3
votes
0answers
39 views

Numerical calculation of the Berry connection

I'm doing some numerical calculations involving Hermitian matrices, and derivatives of the eigenvectors. Essentially, I have an n x n, Hermitian matrix H(x), which is dependent on some continuous ...
3
votes
0answers
85 views

Is the matrix exponential and the Jordan canonical form actually useful for solving differential equations?

All of my yearlong graduate-level Linear Algebra course notes from my professor—an algebraist/representation theorist—shows his love for the exponential map $e^A$ and the Jordan canonical form—and one ...
3
votes
0answers
56 views

How to construct a Fortin Operator for Crouzeix-Raviart Element?

I want to prove the LBB condition for the Stokes Equations discretised by the Crouzeix-Raviart element. The continuous Stokes Equation in the weak formulation is Find $u \in H_0^1(\Omega, \mathbb{R}^...
3
votes
0answers
44 views

Least-squares fit of explicit parabolic sheet to data points

For a given set of data points $$\{(x_i, y_i, z_i)\}$$ there exists some $$f_{ABC}(x,y)=Ax^2+Bxy+Cy^2$$ that minimizes $$\sum_i(f_{ABC}(x_i,y_i)-z_i)^2$$ $A$, $B$, and $C$ can be found quickly ...
3
votes
0answers
43 views

Time independent Runge Kutta integration of SDE

I am trying to compare the result of numerical integration of time independent Runge_Kutta, github page for stochastic differential equations with the analytical solution. True answer match the ...
3
votes
0answers
61 views

Conjugate Gradient for nonlinear equation system

Is it possible to apply adaptions of the conjugate gradient algorithm i.e. Fletcher-Reeves, Polak-Ribere or others to systems of nonlinear equations? How should the equation system be adjusted so one ...
3
votes
0answers
92 views

What is the algorithm to convert an adjacency matrix to a block diagonal structure?

I'm a little at loss as to what my supervisor meant by the "hierarchical block decomposition" of a matrix, but the goal is to put a sparse symmetric adjacency matrix into a block diagonal structure to ...
3
votes
0answers
45 views

Spectral solver on em-pic

I'm recently studying for the spectral solver to implement EM-PIC code. I read an article and have some questions. Many PIC codes uses spectral solver to overcome numerical artifacts on FDTD. In the ...
3
votes
0answers
125 views

Automatically generate constraints for trajectory optimization

This is a follow up to my previous post here I'm interested in performing trajectory optimization from the problem mentioned in abov link. I want to supply the following as dynamical constraints to ...
3
votes
0answers
62 views

how to construct chebyshev differentiation matrix of even/odd polynomial

I know there is efficient way to construct chebyshev differentiation matrix $D_N$. The idea is to fit the function $f(x)$ with a polynomial of $p_N(x)$ of order $N$, and use the derivative $p_N'(x)$ ...
3
votes
0answers
32 views

Is the given equation for Fractional Fourier Transform wrong?

In order to compute $$ \hat{f}_k = \sum_{m = 0}^{M - 1} e^{-2\pi i k m \theta} f_m, \ \ k = 0, ..., M - 1$$ for any $\theta$, my book states that this can be done using fractional Fourier transform ...
3
votes
0answers
65 views

Solving saddle point problem having non-invertible top-left block with a PETSc nested matrix

My system is a symmetric FE problem with lagrange multipliers: $Z=\begin{pmatrix}A & C^T \\ C & 0\end{pmatrix}$ The matrix $A$ is positive semi-definite, non-invertible. The whole matrix is ...
3
votes
0answers
49 views

Hit-n-Run Monte Carlo on convex polytope

So, I'm currently trying to implement a MCMC to uniformly sampling hyper-points from the polytope defined as $\mathbb{K}=\{x\in\mathbb{R}^{n}\;\;\text{s.t.}\;\; A\,x=b \}$ in the specific case where, ...
3
votes
0answers
37 views

Integer partition algorithms

I am familiar with and have written MathCad algorithms for the partition functions 𝑝(𝑛,𝑘),which gives the number of ways of partitioning 𝑛 into 𝑘 parts, 𝑞(𝑛,𝑘), which gives the number of ways ...
3
votes
0answers
57 views

Levinson Recursion for Non Square Toeplitz Matrices

Given a rectangular Toeplitz Matrix $ H $, how could one solve: $$ y = H x $$ For instance, $ H $ can be Linear Convolution Matrix of the filter $ h $: $$ H = \begin{bmatrix} {h}_{1} & 0 & ...
3
votes
0answers
59 views

What is the reason for this finite-difference high errors on non-uniform grid?

tl;dr Using a Taylor-matched method to find coefficients for the discretized equation $ \mathbf{A} \vec{f}'' = \mathbf{B} \vec{f} $, a Fortran code has been implemented to find the second derivative ...
3
votes
0answers
68 views

Computation of Troullier-Martins pseudowavefunctions

The computation of Troullier-Martins pseudowavefunctions has been described in [1]. The pseudowavefunction $R^{\textrm{PP}}_l$ is defined by $$ R^{\textrm{PP}}_l(r) = \left\{ \begin{array}{ll} R^{\...
3
votes
0answers
73 views

Large-scale optimization of nonlinear equations

I'm looking to find a computationally efficient solution to a large system of nonlinear equations. I'm trying to maximize the following function: $$ f(\vec{x}) = \sum_i^N C_i (x_i-A_i)x_i^{\epsilon_{...
3
votes
0answers
51 views

Methods to approximate obective function gradients from point cloud

Problem statement: Assume that I have an objective function $f(x)$ which takes as input a $D$-dimensional vector $x\in\mathbb{R}^D$, and that $f(x)$ is sufficiently smooth. Assume further that I ...
3
votes
0answers
111 views

What is the fastest algorithm for computing log determinant of a PSD matrix? (All possible PSD matrices)

I am diving into some literature to understand which is the best algorithm for computing the log-determinant of a PSD matrix. More generally, I am interested in a list of resources to read, which ...

15 30 50 per page
1
3 4
5
6 7
39