All Questions

Filter by
Sorted by
Tagged with
6
votes
1answer
26 views

Generalization of eigendecomposition problem

Let $A\in \mathbb{R}^{n\times n}$ and $v \in \mathbb{R}^n$. We recognize $Av=\lambda v$ for some scalar $\lambda$ as an eigendecomposition problem. Suppose $\mu \in \mathbb{R}^n$, and let $\odot$ ...
0
votes
0answers
13 views

Mineral dissolution and solute transport around a solid

I am trying to simulate solute transport of acid (HCl) and consequent mineral dissolution around a grain (calcite). The governing equation for transport is the advection-diffusion equation, given as: ...
1
vote
3answers
102 views

Hypergeometric function $_2F_1(z)$ with $|z| > 1$ in GSL

I need to evaluate the hypergeometric function $_2F_1$ with $|z| > 1$ as in Wolfram Language with GSL but the GSL documentation says the $_2F_1$ needs $|z| < 1$. Is there any way I can use GSL ...
3
votes
0answers
59 views

How to position oneself, if one has little practical engineering background or function, but may have a broader insight into mathematical modeling?

How to position oneself, if one has little practical engineering background or function, but may have a broader insight into mathematical modeling? Or, lets say: I could write models to prove ...
0
votes
0answers
16 views

Find tuples of points from multiple sets

Given n sets of points in general position in dimension 2 (n typically small, 2-6), can one find tuples of points, one from each of the sets, which are close in some sense (the closest, mutual ...
0
votes
0answers
19 views

Can line search solve linear objective with nonlinear constraints?

Consider an optimization problem of the form: $$\max_{f(x)\le K,\\0\le x\le M} c^\top x,$$ where $f(x)$ is nonlinear. Can a line search of the following form be used to solve this problem? $$ \max_{\...
2
votes
2answers
361 views

Modelling question: example of a physical phenomenon with this jump condition at an interface?

in our finite element class we were talking about interface problems our teacher came up with the following, where $K_i$ are two given functions and $u_i$ is the restriction of the solution $u$ to $\...
2
votes
1answer
204 views

Integral image resizing

I'm trying to find an approach for integral image resizing. I found out that I can do it with a bilinear interpolation method, but with this approach I can only resizing by the factor which is a power ...
3
votes
1answer
441 views

Methods to implement floor dirt detection algorithm

I'm trying to detect dirty floor areas in a series of images, using MATLAB and its Image Processing Toolbox, like the one that follows: In the image above, there are two distinct areas, the whiter ...
7
votes
1answer
273 views

Robust smoothers for geometric multigrid

I'm searching for robust smoothers for geometric multigrids. By robust I mean: Effective for high order approximations (say spectral element, spectral Discontinuous Galerkin), Parallel (suitable for ...
0
votes
0answers
24 views

Restriction in (geometric) multigrid for vectors of non-even length

Naive restriction operators in geometric multigrid that I have seen are typically implemented as a convolution and a subsequent averaging of every two entries in a vector $v^h$. For example: $$\tilde{...
14
votes
4answers
2k views

What are some good strategies to test a floating point arithmetic implementation for double numbers?

For IEEE, the single representation is 1-bit sign, 8-bit exponent and 23-bit mantissa. This means that at each exponent value, you can test all 2^23-1 (roughly 9mil cases) possible combination of ...
11
votes
1answer
423 views

Increasing V-cycles for constant Coarsest Grid Size and increasing Fine Grid size

Problem statement I implemented geometric multigrid for $-\nabla^{2}=f$ where $f=\frac{3\pi^{2}}{4}sin \frac{\pi x}{2} sin \frac{\pi y}{2} sin \frac{\pi z}{2}$ on $\Omega \in [0,1]$ on a unit cube. ...
0
votes
0answers
30 views

Differences between Two-Grid Correction and V-Cycle Scheme

I found in A Multigrid Tutorial: Second Edition by William L. Briggs‏, Van Emden Henson‏, Steve F. (here) the following Schemes: Two-Grid Correction Scheme, p. 37: And V-Cycle Scheme, p. 40 I want ...
3
votes
1answer
106 views

Calculations on discontinous grids

Suppose for a grid-based calculation a grid is used such that the grid Jacobian is discontinuous. For example, in 1D, for a domain $x \in$ [0,1], one half of the domain is covered uniformly by twice ...
2
votes
0answers
47 views

Comparing block versus non-block Krylov methods for handling multiple right-hand-sides

Suppose I wish to solve a linear system $AX=B$ iteratively where $A$ is an $m\times m$ matrix and $X,B$ are $m \times s $ matrices (not single vectors). Instead of solving $s$ independent systems I'm ...
2
votes
1answer
47 views

How to stack N boxes of varying heights into M stacks, most evenly

The "standard" box-stacking algorithm(s) AFAIK assume a single stack and try to put the "largest" boxes on the bottom. The case I want to solve is simply to distribute the N boxes ...
6
votes
2answers
3k views

How to find the nearest/a near positive definite from a given matrix?

I'm given a matrix. How do I find the nearest (or a near) positive definite from it? The matrix can have complex eigenvalues, not be symmetric, etc. However, all its entries are real valued. The ...
4
votes
1answer
78 views

Worst Case complexity of a search engine algorithm

Computer make it possible to find information in large databases. However, the results are often too large to be returned in their entirety to the user who requests them. Computer therefore sort the ...
2
votes
0answers
84 views

Solve Rational Equation for Root Music in MATLAB

I'm trying to estimate DOA in the Hybrid architecture using root music so I need to solve the attached equation to find the roots for the Root_Music equation in Matlab. Does anyone have an idea for ...
4
votes
0answers
103 views

Any literature that extensively discusses the ability to strongly inter-/extrapolate computationally on little empirical data?

Any literature that extensively discusses the ability to strongly inter-/extrapolate computationally on little empirical data? This topic has fascinated me, but I find that it seems a bit novel. ...
0
votes
1answer
113 views

deal.ii - ParaView "warp by scalar" of my output is not continuous

During our finite element course, we've solved the linear elasticity problem in 2D on a square (GridGenerator::hyper_cube) with $Q_1$ bilinear finite elements in ...
11
votes
1answer
953 views

Numerical computation of Lyapunov exponent

I'm trying to compute the Lyapunov exponent for a smooth continuous time dynamical system(say, $\dot{\bar{x}} = f(\bar x)$). I using the QR decomposition method. Here are the steps that I follow. ...
-1
votes
0answers
47 views

Numerical scheme HJB equation

Without dwelling on details on how to obtain the HJB equation for this problem, I would like to know if the scheme I wrote for solving it numerically is viable or did I miss something. I need to solve ...
6
votes
2answers
103 views

Approximating the boundary between two sets of points (in 2D): Fitting a region

Given two sets of points $p_{\text{in},i}$ and $p_{\text{out},j}$ inside and outside of what I intuitively call a "region", I would like to estimate and describe the boundary of this region. ...
0
votes
1answer
123 views

Specifying mesh spacing for DFT in numpy

I was testing the .fft package of numpy 1.16.1 in Python 3.7.2. In particular I was trying to verify that the transform resembles the analytical one for: $$f(x) = \mathrm{exp}\left[-\left(\frac{x-5}{2}...
6
votes
0answers
70 views

How do we approximate the numerical error a numerical scheme (e.g Runge Kutta, Euler etc) makes without having access to an analytical solution?

So I recently encountered this question in my head while taking my Scientific Computing class, where the lecturer talked about computing numerical error of a scheme. My guess would be that we take a ...
3
votes
1answer
141 views

Is there any way/any python function to calculate the condition number of the roots of a polynomial directly?

I know that NumPy has linalg.cond(A) to find the condition number of a matrix A. But, if I want to find the condition numbers of the roots of a large polynomial ...
0
votes
1answer
57 views

Finding total derivative of a multivariate function in Maple

In Maple, I have a function $f(x(t),y(t),t)$ that I want to differentiate with respect to $t$. I know the command for partial derivative $\frac{\partial f}{\partial x}$,$\frac{\partial f}{\partial y}$,...
1
vote
1answer
72 views

Can someone explain the equivalence between Oja's rule and PCA in a simple way?

I have to give a presentation on unsupervised learning in 2 days, and I have to explain/show the equivalence between Hebb's learning rule (or Oja's rule to be more specific) and PCA. The thing is that ...
1
vote
1answer
394 views

Suitable finite difference method for a convection-diffusion system?

I am trying to solve a system of PDEs $H_{t} = \frac{0.3}{0.7} - \frac{0.005 B f(h(H))}{\theta} - \frac{0.3 f(h(H))}{0.7} + \frac{500}{0.7} (HH_x)_x + (HH_y)_y$ $N_t = \frac{N_{in} - 0.002 [N] B f(h(...
1
vote
1answer
70 views

Elementary matrix of Raviart-Thomas elements

We can use the $RT0$ to solve the Darcy equation, i.e. $$k^{-1}\mathbf{u}+\nabla p = 0, \text{ in } \Omega,$$ $$-\nabla \cdot \mathbf{u} = 0, \text{ in } \Omega,$$ $$p = p_D \text{ on } \partial\Omega,...
-1
votes
0answers
35 views

Identifying Noetherian symmetries from general Lie symmetries of a differential equation

I know the Lie symmetry group of the harmonic oscillator differential equation from the literature. It is an 8 parameter Lie group. 5 of these generators generate Noether symmetry: $$G_1=\sin(2t)\frac{...
0
votes
0answers
54 views

Contact analysis does not converge due to the projection falls outside valid domain

I implemented Node-To-Surface contact algorithm (Wriggers, Peter, Computational contact mechanics., Berlin: Springer (ISBN 3-540-32608-1/hbk). xii, 518 p. (2006). ZBL1104.74002.). The code is done by ...
0
votes
0answers
25 views

How to use plot command in fenics [closed]

I'm new to Fenics and want to understand the plot command. My code: ...
11
votes
2answers
1k views

How do I find the minimum-area ellipse that encloses a set of points?

I have a set of points that resembles more of an ellipse than a circle. I implemented the optimization formulation below and the solution gives a circle. I tried with various initial values, still to ...
0
votes
0answers
48 views

Solve simultaneous differential equations with embedded functions and a parameter estimation

The aim is to solve the below equations and plot $m$ with time, i.e. $\frac{dm}{dt}$ $k$ is unknown and needs to be estimated. For the parameter estimation, the below values in the table for m versus ...
2
votes
2answers
133 views

Test functions of Raviart-Thomas elements?

The test functions of general finite elements are like interpolation functions (if my understanding is correct). But how about test functions of Raviart-Thomas elements? Let's raise the $RT0$ element ...
2
votes
1answer
34 views

Bounds condition for IFT to obtain a $1/f$ time-series

I am coding a function to obtain a randomized time-series from a given $\frac{1}{f}$ law. The randomization is obtained by introducing a random phase in the function. I experience a problem in the ...
0
votes
1answer
71 views

1D wave equation using Finite difference method MATLAB

I have the wave equation $$u_{tt} = 4 u_{xx}$$ with the boundary conditions $$u(0,t) = u(L,t) = 0\,,\quad x \leq 0 \leq 2\pi \,,\quad t\geq 0$$ and initial conditions $$\begin{align} &u(x,0)=\...
2
votes
0answers
38 views

Appropiate Artificial Boundary Conditions for the radial part of the Klein Gordon equation?

I am trying to simulate the following equation using FDTD $ \left(- \partial^2_t + \partial^2_x + V(x) \right) \psi(x,t) =0 $ subjected to the initial conditions $\psi(x,0) = f(x),~ \partial_t \psi(x,...
0
votes
0answers
165 views

Interpolation question

I have a set of data $(x_i,y_i,y'_i)$, $i=1,\dots,N$ and I want to fit an interpolating curve $f(x)$ which matches both the data $y_i$ and the first derivatives $y'_i$ at the nodes $x_i$, \begin{align}...
5
votes
3answers
228 views

Numerical integration giving incorrect sign

For my research, I need to integrate the following function: $$ W(z)=\int_0^{\infty}dx\ w(x,z)\\ =\int_0^{\infty}dx\frac{e^x}{(e^x+1)^2}\log{\left(\frac{e^{z^2/4x+x}+1}{e^{z^2/4x+x}-e^x}\right)}\\ =\...
1
vote
0answers
50 views

ICCG negative residual products $r^TM^{-1}r$

I have a linear system $Ax=b$ resulting from a finite element discretization of the Poisson equation. I am applying an IC0 (incomplete Cholesky ($LDL^T$) with the same sparsity as the original matrix) ...
1
vote
1answer
55 views

solve_ivp doesn't work with toms748

I have the following code ...
1
vote
1answer
123 views

A better word to indicate slowness/high latency?

We are comparing two techniques in computer science. We want to say X has "significantly high latency" when executed on system Y. Is there a better one-word term we can use for the above to ...
0
votes
0answers
26 views

Conserved current for a 1d oscillator using Maple

I found out the conserved current for harmonic oscillator with angular frequency 1, particle falling in gravity close to ground ($g=1$) using maple. I'm unable to understand the result: $$J[t](t, x(t)...
5
votes
1answer
198 views

Eigen C++ Library Over 4x Slower than Fortran

I am experimenting with ways to call compiled programs from Python. My primary interest is iterative optimization methods, so I'm testing an implementation of Newton's method to solve a logistic ...
1
vote
1answer
78 views

SciPy odeint giving different results with matrix multiplication

I've asked this at stackoverflow but maybe this community will have a better idea of the answer. I'm currently trying to develop a function that performs matrix multiplication while expanding a ...
1
vote
1answer
77 views

Numerical Solution of the Schrödinger equation for hydrogen

I'm trying to solve the Schrödinger equation for the hydrogen atom in the following form numerically: $$\left[-\frac{\hbar^2}{2m}\frac{d^2}{dr^2}+V(r)+\frac{\hbar^2l(l+1)}{2mr^2}\right]R(r)=ER(r).$$ ...

15 30 50 per page
1
2 3 4 5
200