The tag has no wiki summary.

learn more… | top users | synonyms

0
votes
0answers
17 views

Bootstrap for a histogram

I create a set of $T$ trajectories with $P$ positions, $\{x_j\colon 0 \leq j < M\}$, with a Monte Carlo simulation. From this data, I calculate quantities like $\langle x \rangle$ and $\langle x^2 ...
12
votes
4answers
254 views

Estimating hardware error probability

Say I run a supercomputer computation on 100k cores for 4 hours on http://www.nersc.gov/users/computational-systems/edison/configuration, exchanging about 4 PB of data over the network and performing ...
4
votes
0answers
31 views

Methods for integrating black box functions on a non-uniform grid

If i have some function expressed as points on a non-uniform grid (I'm specifically interested in logarithmic grids, but general results are also interesting), and I want to integrate it, I believe ...
3
votes
1answer
39 views

a posteriori error estimation for skewed elements

I'm working with error estimates for Poisson's equation of the form $$\mathcal{E}^2_T = h_T^2\|-\Delta u - f\|_{L^2(T)} + \sum_{e\in \partial T} h_e\|n\cdot \nabla u\|_{L^2(e)}$$ where $T$ is an ...
0
votes
0answers
41 views

Approximating forward Error function

and i have i question. i was given equation $$f(x) = 0$$ $$f(x) = cos(\frac{x}{50}) - \frac{1}{\sqrt{2}}$$ and the approximation root $$x_a $$ such that $$\vert f(x_a)\vert < \epsilon ...
6
votes
2answers
177 views

Proving convergence of adaptive finite elements - min res FEM?

There's a body of work out there dealing with the discrete convergence of adaptive finite element methods using error estimators. Most deal with proving the property $\|u-u_{k+1}\|_U \leq (1-\alpha) ...
2
votes
1answer
106 views

Error norm in finite difference calculation

I've used an explicit finite difference scheme to model the 1D time dependent temperature distribution in a friction weld. I want to now verify the consistency and convergence of my algorithm. I have ...
6
votes
1answer
163 views

Computing $\sin(\pi/2)$ numerically

I'm trying to understand the types of numerical errors, to do this I want to calculate $\sin(\pi/2)=1$ numerically. To do this I use the Taylor series of $\sin(x)$ in 0: ...
2
votes
2answers
55 views

Normalizing error when data passes through zero

Given a time-dependent measure of error, $e(t)$, and a time-dependent function to normalize it with, $f(t)$, one would report the normalized error as $$E(t) = \left|\frac{e(t)}{f(t)}\right|.$$ ...
4
votes
2answers
177 views

How can you calculate percent error in tensor approximations?

I have a matrix A which is an approximation to the known matrix B. Both matrices are square, 3x3 matrices and, in this case, are ...
3
votes
1answer
30 views

Known Algorithm to compute errors of given nature

I have three sets of data, measured by three different devices: A,B and C of air balloon whose fall is influenced by wind and ...
7
votes
1answer
127 views

Bounding the relative error of derivative given relative error of the function

Suppose a function $f$ can be computed such that the bound on the relative error is $R$ i.e. $f^-(x) = f(x)(1+r)$ where $f^-$ and $f$ are respectively the computed and exact value $f$ and $|r| \leq ...
5
votes
1answer
67 views

Can we compare the speed of convergence of two different iteration methods of same order looking at their error estimates?

I have a two iterative method for approximating the inverse of given square matrix $A$ whose error terms are given as follows Error estimate of method $1$: $\lVert A^{-1} - X_{k}\rVert \leq ...
4
votes
1answer
105 views

Error bars for pair-correlation function

I have obtained some data from neutron diffraction for some material samples. The "rawest" form of the data is the structure function $S(Q)$. We can choose a variety of different Q-maxes when ...
1
vote
1answer
99 views

How can I quantify the error of FFT-based poisson solvers?

I have an FFT code that solves a particular case of the steady Euler equations where a Poisson equation is solved, what is a good way to quantify the error? Is what I am doing ok? Since I do not have ...
4
votes
2answers
149 views

Problem Condition and Algorithm Stability

Consider 2 mathematical problems: $$ f_1(x) = a - x \\ f_2(x) = e^x -1 $$ The condition number for a function is defined as follows: $$ k(f) = \left| x \cdot \frac{f'}{f} \right| $$ Lets analyze ...
16
votes
2answers
354 views

Why do equi-spaced points behave badly?

Description of experiment: In Lagrange interpolation, the exact equation is sampled at $N$ points (polynomial order $N - 1$) and it is interpolated at 101 points. Here $N$ is varied from 2 to 64. ...
1
vote
1answer
95 views

Error analysis of WENO scheme

I have three questions regarding WENO schemes 1) How to actually compute the smoothness indicators $\beta_j$ for required order of polynomial? Any reference which explains the algorithm will be ...
12
votes
3answers
304 views

Numeric Quadrature with Derivatives

Most numerical methods for quadrature treat the integrand as a black-box function. What if we have more information? In particular, what benefit, if any, can we derive from knowing the first few ...
2
votes
0answers
31 views

Propogated Error in Mesh Interpolation

I am working with a code that solves diffusion/reaction equations on a 2D unstructured mesh. Due to the stiffness of some of the processes, I start with time steps near 1e-13, and end with a final ...
2
votes
1answer
102 views

Diffusion-Transport problem FEM

I was looking at a book of FEM on problems of Diffusion-Transport. $$-div(\mu \nabla u) + b \cdot \nabla u + \gamma u = f \qquad in~\Omega$$ $$u = 0 \qquad in~\partial\Omega\text{ (in the ...
31
votes
4answers
772 views

Scientific standards for numerical errors

In my field of research the specification of experimental errors is commonly accepted and publications which fail to provide them are highly criticized. At the same time I often find that results of ...
5
votes
1answer
106 views

a priori error analysis of cell-centered finite-volume methods

I'm using the cell-centered finite-volume method (for example Morton, Numerical solution of convection-diffusion problems, Chapman&Hall, 1996) to discretize the advection-diffusion equation and ...
2
votes
0answers
76 views

function over conditional probability

I need to create a scoring model out of estimated conditional probability functions for two events, A and B. Let 0.5 be the threshold value. Ideally, the probability is in the interval $[0,0.5)$ for A ...
1
vote
1answer
267 views

How to find the number of principal components that lead to the smallest generalization error?

I am working on a paper part of which is the application of validation rules to find how many principal components give us the least generalization error. The concept goes more or less like this: ...
5
votes
2answers
855 views

When to stop Gauss-Seidel-iterations?

I want to have an estimation, that my solution has an error, let's say less than 1e-8. Usually, I stop the Gauss-Seidel algorithm, when the residual is "small enough" and this is already the problem. ...
5
votes
3answers
1k views

What norm to choose when?

Recently, I saw this question: how to measure the error of a finite difference method I am student of simulation sciences and unfortunately, for me, it's totally unclear, what norm to use in what ...
3
votes
1answer
161 views

Integral average approximation and error bounds

I'm looking into integrals of the form: $$\int_a^b {f(x)g(x)dx}$$ Where $f(x)$ is unknown, but it's integral is: $$\int_a^b {f(x)dx}=F$$ It's been suggested to me that one could approximate this ...
9
votes
2answers
174 views

What about this simple error estimate for linear PDE?

Let $\Omega$ be a convex polygonally bounded Lipschitz domain in $\mathbb R^2$, let $f \in L^2(\Omega)$. Then the solution of the Dirichlet problem $\Delta u = f$ in $\Omega$, $\operatorname{trace} u ...
4
votes
1answer
187 views

finite difference methods and global error

I was going through my notes on different finite difference methods and came across something I don't quite understand. I have code that will calculate an approximate solution we can call this ...