Questions tagged [error-estimation]
For questions about determining the error caused by specific computational procedures, approximations, or numerical representations.
168
questions
0
votes
0
answers
39
views
What is the minimum error achievable using gaussian process emulation?
I am interested in using Gaussian processes as emulators for other computational models, and I would like to characterize the expected numerical precision of the emulator.
Specifically, how small can $...
4
votes
2
answers
425
views
Understanding the diffusion error of numerical schemes
In chapter 7 of Numerical Computation of Internal and External Flows (Second Edition) (https://www.sciencedirect.com/science/article/abs/pii/B9780750665940500497) the author describes the diffusion ...
0
votes
1
answer
60
views
Reverse engineering phase shift and numerical damping
I've been trying to validate the physics behind a particle system framework, but I'm having some difficulties.
A particle system is a set of lumped masses connected by spring-damper elements. Linear ...
4
votes
1
answer
132
views
The error propagation in calculating the inverse using a matrix decomposition
I have been trying to calculate the matrix inverse of some large matrix with entries ranging by orders of magnitude. I tried to use the matrix decomposition to simplify the computation, where a matrix
...
1
vote
1
answer
158
views
Storing Raw Simulation Data or Truncated Data?
I have a simulation that can generate quite a bit of data when it runs, for example $650\cdot 400 \cdot 400$ floating point numbers. Without compression, that's a few gigabytes worth if I want to save ...
0
votes
1
answer
72
views
Finite difference approximation error
I was reading Scientific Computin, An Introductory Survey, by Michael Heath. In the Example 1.11, he madr a Finite Difference Aproximation, with the usual approxination : $f’(x)\neq \frac{f(x+h)-f(x)}{...
1
vote
0
answers
116
views
Which dense matrices are hard to invert?
Suppose I'm solving $Ax=b$ for dense $m\times d$ matrix $A$. For which $A$ is this hard to do?
More concretely, is there any work on estimating the error after $k$ steps of iterative solver, $k\le d$, ...
1
vote
1
answer
77
views
Geometrically nonlinear finite element problem and mesh distortion
In my research on topology optimization of fluid-structure interaction problems (2D), I am using a geometrically nonlinear model to represent the structure. Body/surface forces are extracted from the ...
0
votes
0
answers
41
views
How to compare the bias between the two contour plots?
I wish to compare the value of the normalized concentration (c+), between two contour plots (one is experimental and another is a simulation). If I don't have a c+ value from the same point on both ...
0
votes
0
answers
55
views
Finding optimal values from multiple parameter estimation runs
I've performed a parameter estimation repeat (i.e. 1000 parallel runs with the same
initial values of parameters). I am trying to estimate ~20 parameters using measurements from experiments.
After ...
0
votes
1
answer
64
views
Convergence of FEM on curved boundaries, and inhomogenous boundary data
In a smooth domain in $\mathbb{R}^{2}$ or $\mathbb{R}^{3}$ let's consider $-\Delta u = f$ with $u=g$ on a part of the boundary and $\partial_\nu u = w$ on another part of the boundary, which is far ...
0
votes
1
answer
87
views
finding discretization error in Burger equation
I was reading the paper given in the link http://www.unige.ch/~hairer/preprints/parareal.pdf and I have a problem in understanding in page 10 for the Burger equation on implementing Parareal method ...
7
votes
1
answer
305
views
How does non-dimensionalization improve the behavior of ODE solvers?
I have a set of coupled ODEs that I'm solving numerically. The independent variable is time and runs from values of $10^{15}$ to $10^{17}$ in units of seconds. The state variables in their usual ...
1
vote
1
answer
77
views
A priori estimates in finite elements for inhomogeneous heat equation
Consider the problem
$$\partial_t u-\Delta u = f\\
u(\Sigma_1)=f_D\\
\partial_\nu u (\Sigma_2)=f_N\\u(0)=u_0$$
where the sides of the space-time cylinder $\Sigma_i$ are disjoint (one of them could be ...
3
votes
0
answers
60
views
How to maximize effectiveness of 2 check bits overseeing 6 data bits?
Having a single parity bit gives a Hamming distance of 2, so only one single bit can be corrupted to reliably detect an error.
The next step I know of would be a Hamming code, which uses at least 3 ...
2
votes
1
answer
130
views
"A posteriori" estimates for finite difference methods
Suppose I have a PDE in a rectangular domain that I am solving numerically via a finite difference method. How do I answer the question, "How fine do I need to make the grid to so that my ...
2
votes
1
answer
182
views
Accuracy loss in single-precision Euclidean norm computation
I do hydrodynamics simulations with Fortran and recently I met with this issue:
I have a single-precision array b of length ...
-1
votes
1
answer
46
views
What is the ERRCON parameter in rkqs?
Ive take a course in computational physics and was asked to implement some numerical methods to solve ODES. I was reading up on the algorithms described in the textbook: NUMERICAL RECIPES IN FORTRAN.
...
0
votes
2
answers
89
views
calculating the Laplacian of the field variable in estimating the local residual error in the finite element method
to perform adaptive refinement in the finite element method according to the explicit residual method, the quantity
$$\eta_K^2=h_K^2\left\lVert r\right\rVert_{L_2(K)}^2+h_K\left\lVert R\right\rVert_{...
0
votes
1
answer
55
views
Numerical Error source when dealing with integer series
I am currently trying to compute the value of the first Fibonacci number recursively. the idea is as follow:
Compute $f_{n}$ and $f_{n-1}$ for $n = 2,...,100$,
Compute $f_k$ for $k = n−2, n−3, \dots, ...
8
votes
0
answers
102
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 ...
4
votes
1
answer
97
views
Roundoff errors in FEM computations - generalized eigenvalues
This is a continuation of my previous question. I am trying to effectively compute a bound for the roundoff errors in some FEM computation (2d polygons, triangular meshes). Below I will write some of ...
2
votes
2
answers
546
views
Different sources of error in Finite Element computations
Consider the problem $-\Delta u = f$ in $\Omega$, with $u=0$ on $\partial \Omega$. Suppose that $\Omega$ is a polygon and that we approximate the solutions to the previous problem using Lagrange ...
3
votes
1
answer
191
views
How to measure the error of Finite Element approximation in satisfying the PDE?
In Galerkin methods, we seldom can measure the accuracy of an approximation by tracking the value of the residual. For example, take the wave equation:
$ u_{tt} = u_{xx}$, $(x,t) \in (0,L) \times (0,T)...
1
vote
1
answer
287
views
Question about step in the proof of standard discrete trace inequality
I'm studying from Guermond lecture notes available at https://www.math.tamu.edu/~guermond/M661_FALL_2019/chap12.pdf (see Lemma 12.8( Discrete trace inequality).)
Consider the simple case $p=r$, i.e. ...
4
votes
0
answers
85
views
Global reconstruction defined elementwise in a-posteriori error estimator
This question is a follow-up of this previous one. In "Error Control for Discontinuous Galerkin Methods for First Order Hyperbolic Problems" by Georgoulis et al., an error estimator is ...
0
votes
1
answer
107
views
How is the integral of a projection over an element $T$ computed in practice? (deal.II related)
I'm studying an error estimator for the equation $\nabla\cdot(\beta u) + cu = f$ and it contains the following term $$||f - cU_h - \Pi(f-c U_h) ||_T$$
where :
$\Pi$ is the local orthogonal $L^2$ ...
2
votes
1
answer
4k
views
Error in Simpson's 3/8 rule is higher than that of Simpson's 1/3 rule
For a given function $f(x)$, I have tried to find its numerical integral using Simpson's 1/3 and Simpson's 3/8 rules.
I then compare the solution from the numerical quadratures to the analytical ...
3
votes
0
answers
72
views
Typo in a-priori error estimate in a Discontinuous Galerkin paper
I'm looking at this famous paper which is available in the link below:
Franco Brezzi, LD Marini, Endre Süli, Discontinuous Galerkin methods for first-order hyperbolic problems, Mathematical Models ...
1
vote
1
answer
116
views
Finding the source of numerical instability in a electrostatic problem solved by conformal mapping
I'm using conformal mapping to solve a 2D electrostatic problem (calculating the potential $u(x,y)$ in the plane). Let $C_1$ and $C_2$ be two circles at an electric potential $U_1$ and $U_2$, ...
3
votes
1
answer
151
views
Proof of R. Verfürth paper on adaptive mesh and bubble functions
I'm studying adaptive meshes, and my professor wrote the following property for a bubble function ( see this scicomp post for the definition I'm using)$b_T$ defined on a triangle $T$.
$$||b_T \phi ||_{...
1
vote
1
answer
104
views
Classical global estimate for $H^1$ error
I'm having lots of troubles in understanding the proof the estimation of the classical $H^1$ error using finite elements of degree $r$.
$$||u-u_h||_{H^1(\Omega)} \leq \frac{M}{\alpha} C h^r |u|_{H^{r+...
3
votes
1
answer
275
views
How to find the optimum finite difference method for derivatives?
Related to: What are the negatives of using higher order finite diference schemes?
Problem: I have some discrete data of a trajectory $x_t$ with errors $\delta x_t$ of a physical system sampled at ...
4
votes
1
answer
145
views
Computing the residual in a Dual Weighted Residual (DWR) method
I am in the process if computing the Dual-Weighted Residual (DWR) for a linear PDE with a linear functional but I am struggling with the residual part of the calculation.
For example suppose we want ...
1
vote
2
answers
214
views
Stability condition FCTS method
The FTCS method comes from the discretization of a diffusion PDE like this:
$$
a^{2} \frac{u_{i+1}^{k}-2 u_{i}^{k}+u_{i-1}^{k}}{\Delta x^{2}}=\frac{u_{i}^{k+1}-u_{i}^{k}}{\Delta t}
$$
If I have the ...
0
votes
1
answer
52
views
Round-off error step choice
In the Numerical Recipes in section 5.7.- Numerical derivatives the choice of the step size $h$ in the numerical derivative should lead to a difference between $x$ and $x+h$ representable by an exact ...
0
votes
1
answer
204
views
Taylor expansion round-off error
In the Numerical Recipes in section 5.7.- Numerical derivatives it's introduced de roundoff error of:
$$
f^{\prime}(x) \approx \frac{f(x+h)-f(x)}{h}
$$
as (with $h$ an "exact" number):
$$
\...
2
votes
1
answer
65
views
Initial condition precision
Is there a way to have an estimate of the error propagated on an ode numerical solution by the error of the initial conditions? I suppose this depend on the numerical method used and on the problem ...
1
vote
1
answer
279
views
FEM solution for Poisson is not exact at nodes
Let's consider the usual Poisson problem for FEM $$-u''(x)=1 \quad x \in [0,1]$$ with homogeneous Dirichlet boundary conditions. The solution is $u(x)=-\frac{1}{2}x(x-1)$
I know that the FEM solution (...
0
votes
1
answer
201
views
Perturbation problem using Runge-Kutta 4
I'm trying to evaluate the perturbations magnitude between 2 body orbiting a central one in three dimensions. In order to do this I need to have an estimate of the error, which I did using Richardson ...
2
votes
1
answer
228
views
Comparing numerical solutions with very different time grids
I've read an article (Long-term integrations and stability of planetary orbits in our Solar system) in which the authors solved the problem of the absence of an analytical solution for the solar ...
0
votes
1
answer
229
views
Error too large in leapfrog method for solving the wave equation of a vibrating string
I have been trying to figure out what I did wrong for the last two days. I do not know if I actually did something wrong or if the error is supposed to be this large in usual leapfrog problems. I ...
0
votes
1
answer
200
views
What are the advantages and disadvantages of using norm error control in the MATLAB ODE suit?
In MATLAB's ODE suit, there seem to be two basic methods of controlling the Local Truncation Error (LTE) of the ODE which the user can choose from, namely:
The absolute error control (default), ...
4
votes
2
answers
89
views
Dividing functions over a wide range-
I try to solve a system of coupled equations, where a very nasty division operation occurs. In fact, I need to compute a derivative of two exponential decaying functions. Let's illustrate this with ...
1
vote
1
answer
523
views
How to compute the $L^{2}$ error of the gradient in the Finite Element Method
Let $\Omega\subset \mathbb{R}^{2}$ and $\tau_{h} = \{\Omega_{k}\}_{k=1}^{N}$ be a triangulation of $\Omega$. The $L^2$ error for a FEM approximation $u_{h}$ is given by:
$ || u-u_{h} ||_{L^2} = \sqrt{ ...
3
votes
1
answer
106
views
Cauchy Lorentzian simulation on FFT with oscillation
Recently I do simulation on Lorentzian Function with FFT
Lorentzian Function is 2a/(x**2+a**2)
...
0
votes
0
answers
133
views
Error in Monte Carlo integration
I am looking for a concise description to help me understand the error for Monte Carlo Integration using Uniform Sampling and Importance Sampling
For Importance Sampling I have that the error is just ...
3
votes
1
answer
857
views
Calculating Error for Poisson Equation using Successive Over-Relaxation technique, Python
I am trying to solve the Poisson Equation
$\frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} = 32(x(x-1) + y(y-1))$
for a 61x61 grid using Python3 with boundary conditions being $T=...
1
vote
1
answer
69
views
Recovery of smoothed continuous stresses using the Z^2 error estimator
I have trouble implementing the $Z^2$ error estimator (The superconvergent patch recovery and a posteriori error estimates. Part 1: the recovery technique by Zienkiewicz and Zhu).
For this, I am ...
2
votes
0
answers
411
views
Implementation of Z^2 error estimator in Abaqus for adaptive mesh refinement
Currently, I am working on a remeshing routine for my simulations (Abaqus 6.14-1) using python scripts. The simulation deals with the Brinell indentation test and as the remeshing software I use Gmsh ...