Orthogonal Basis Function Approximation of Particle Distributions In Numerical Simulations of Beams
|
|
- Earl Melvin Page
- 5 years ago
- Views:
Transcription
1 Orthogonal Basis Function Approximation of Particle Distributions In Numerical Simulations of Beams Balša Beam Physics and Astrophysics Group Northern Illinois University Jefferson Lab Beam Seminar April 3, 008
2 Motivation Studying dynamics of multi particle systems (charged particle beams, plasma, galaxies...) heavily relies on N body simulations It is important for N body codes to: be as efficient as possible, without compromising accuracy minimize numerical noise due to Nsimulation<< Nphysical account for multiscale dynamics for some applications: have a compact representation of history We present two orthonormal bases which, as a part of an N body code, address these requirements wavelet basis scaled Gauss Hermite basis
3 Outline of the Talk Algorithms for N body simulations Wavelet basis brief overview of wavelets wavelet based Poisson equation solver advantages applications Scaled Gauss Hermite basis mathematical formalism Poisson equation solver applications Discussion of further work 3
4 Algorithms for N body Simulations Direct summation: CPU cost scales as N Tree: direct summation nearby and statistical treatment farther away Particle In Cell (PIC): particles binned in cells (grid) system at t=t0 bin macroparticles to obtain particle distribution on a finite Nx Ny Nz grid N macroparticles Newton's equations advance particles by t system at t=t0+ t interpolate to find F for each particle F= solve the Poisson equation = on a finite Nx Ny Nz grid 4
5 Algorithms for N body Simulations Direct summation: CPU cost scales as N Tree: direct summation nearby and statistical treatment farther away Particle In Cell (PIC): particles binned in cells (grid) system at t=t0 bin macroparticles to obtain particle distribution on a finite Nx Ny Nz grid N macroparticles Newton's equations advance particles by t system at t=t0+ t interpolate to find F for each particle F= solve the Poisson equation = on a finite Nx Ny Nz grid use wavelets 5
6 Algorithms for N body Simulations Alternative N body algorithm: analytical function approximation analytical functions form a finite orthogonal basis N macroparticles, but no grid system at t=t0 expand particle distribution in a finite orthogonal basis (x,y,z)= ψ (x,y,z) N macroparticles ijk Newton's equations advance particles by t differentiate to get F for each particle F = system at t=t0+ t ijk ijk solve the Poisson equation = in the same basis (x,y,z)= ψ (x,y,z) ijk ijk ijk 6
7 Algorithms for N body Simulations Alternative N body algorithm: analytical function approximation analytical functions form a finite orthogonal basis N macroparticles, but no grid use scaled Gauss Hermite basis system at t=t0 expand particle distribution in a finite orthogonal basis (x,y,z)= ψ (x,y,z) N macroparticles ijk Newton's equations advance particles by t differentiate to get F for each particle F = system at t=t0+ t ijk ijk solve the Poisson equation = in the same basis (x,y,z)= ψ (x,y,z) ijk ijk ijk 7
8 Wavelets Wavelets: orthogonal basis composed of scaled and translated versions of the same localized mother wavelet ψ(x) and the scaling function ф(x): ψik (x)=k/ψ (kx i) f(x) =s00ф00 (x) + i k dikψik (x) Discrete Wavelet Transfrom (DFT) iteratively separates scales ~O(MN) operation, M size of the wavelet filter, N size of the signal Advantages: simultaneous localization in both space and frequency compact representation of data, enabling compression (FBI fingerprints) signal denoising: natural setting in which noise can be partially removed denoised simulation simulation with more macroparticles 8
9 Numerical Noise in PIC Simulations Any N body simulation will have numerical noise Sources of numerical noise in PIC simulations: graininess of the distribution function: Nsimulation<< Nphysical discreteness of the computational domain: and specified on a finite grid Each macroparticle is deposited onto a finite grid by either: Nearest Grid Point (NGP) dep. scheme Cloud In Cell (CIC) dep. scheme NGP CIC x marcoparticle location gridpoint location 9
10 Numerical Noise in PIC Simulations For NGP, at each gridpoint, particle dist. is Poissonian: n nj is the expected number in jth cell; P= n! n j e 1 n j For CIC, at each gridpoint, particle dist. is contracted Poissonian: n P= n! an j e 1 n integer an j D/ a= / 3 ~0.54 3D, 0.67 D, 0.8 1D Measure of error (noise) in depositing macroparticles onto a grid: = N grid 1 N grid Var qi i=1 NGP = Q total N N grid CIC = a Q total N N grid where qi= (Qtotal/N)ni, Qtotal total charge; Ngrid number of gridpoints (For more details see, Pogorelov & Bohn 007, PR STAB, 10, 03401) This error/noise estimate is crucial for optimal wavelet denoising IDEA: Solve the Poisson equation in such a way so as to minimize numerical noise USE WAVELETS 10
11 Numerical Noise in PIC Simulations In wavelet space: signal noise few large wavelet coefficients cij many small wavelet coefficients cij Denoising by wavelet thresholding: if cij < T, set to cij =0 (choose threshold T carefully!) A great deal of study has been devoted to estimating optimal T T = log N grid (σ was estimated earlier), Pogorelov & Bohn 007, PR STAB, 10,
12 Wavelet Denoising and Compression Whenever the discrete signal is analytically known, one can compute the Signal to Noise Ratio (SNR) which measures its quality SNR ~ Nppc ANALYTICAL Nppc: avg. # of particles per cell Nppc= 3 SNR =.0 Nppc = N/Ncells Nppc= 05 SNR =
13 Wavelet Denoising and Compression Whenever the discrete signal is analytically known, one can compute the Signal to Noise Ratio (SNR) which measures its quality SNR ~ Nppc Nppc: avg. # of particles per cell Nppc = N/Ncells D superimposed Gaussians on grid ANALYTICAL Nppc=3 SNR=.0 Nppc=05 SNR=
14 Wavelet Denoising and Compression Whenever the discrete signal is analytically known, one can compute the Signal to Noise Ratio (SNR) which measures its quality SNR ~ Nppc Nppc: avg. # of particles per cell Nppc = N/Ncells D superimposed Gaussians on grid ANALYTICAL Nppc=3 SNR=.0 Nppc=05 SNR=
15 Wavelet Denoising and Compression Whenever the discrete signal is analytically known, one can compute the Signal to Noise Ratio (SNR) which measures its quality SNR ~ Nppc Nppc: avg. # of particles per cell Nppc = N/Ncells D superimposed Gaussians on grid ANALYTICAL Nppc=3 SNR=.0 Nppc=05 SNR=
16 Wavelet Denoising and Compression Whenever the discrete signal is analytically known, one can compute the Signal to Noise Ratio (SNR) which measures its quality SNR ~ Nppc Nppc: avg. # of particles per cell Nppc = N/Ncells D superimposed Gaussians on grid ANALYTICAL Nppc=3 SNR=.0 Nppc=05 SNR=16.89 denoising by wavelet thresholding: if cij < T, set to 0 16
17 Wavelet Denoising and Compression Whenever the discrete signal is analytically known, one can compute the Signal to Noise Ratio (SNR) which measures its quality SNR ~ Nppc Nppc: avg. # of particles per cell Nppc = N/Ncells D superimposed Gaussians on grid COMPACT: only 0.1% of coeffs WAVELET THRESHOLDING ANALYTICAL Nppc=3 SNR=.0 Nppc=05 SNR=16.89 denoising by wavelet thresholding: if cij < T, set to 0 Advantages: DENOISED Nppc= 3 SNR = increase in SNR by c c more macroparticles (here c=8.3, c=69) compact storage in wavelet space (in this example 79/65536: 0.1%) 17
18 Wavelet Based Poisson Equation Solver Poisson equation in physical space u=f BCs: using Green's functions ubnd = g Preconditioned Conjugate Gradient in wavelet space A X B ε B DWT 1 solution U on the Nx Ny Nz grid transform discretized Poisson eq. to wavelet space discretize Poisson equation on a Nx Ny Nz grid LU=F DWT Lw Uw = Fw k (L) ~ O(Nx) wavelet threshold source Fw, operator Lw AX=B (PLwP) P 1Uw = PFw A compressed operator PLwP Precondition Laplacian Lw B denoised distribution PFw with diagonal preconditioner X denoised solution P 1Uw P k (Lw) ~ O(Nx) k (PLwP) ~ O(Nx) physical space wavelet space 18
19 Preconditioned Conjugate Gradient (PCG) i convergence rate depends on condition number k u u i k 1 u k ~ O(Nx) preconditioning (diagonal in wavelet space): good initial approximation: solution at previous time step k 1 k ~ O(Nx ) 19
20 Preconditioned Conjugate Gradient (PCG) i convergence rate depends on condition number k u u i k 1 u k ~ O(Nx) preconditioning (diagonal in wavelet space): good initial approximation: solution at previous time step k 1 k ~ O(Nx ) typical realistic beam simulation no preconditioning Ui=0 initial guess average over step run 75. 0
21 Preconditioned Conjugate Gradient (PCG) i convergence rate depends on condition number k u u i k 1 u k ~ O(Nx) preconditioning (diagonal in wavelet space): good initial approximation: solution at previous time step k 1 k ~ O(Nx ) typical realistic beam simulation average over step run no preconditioning Ui=0 initial guess 75. preconditioned Ui=0 initial guess
22 Preconditioned Conjugate Gradient (PCG) i convergence rate depends on condition number k u u i k 1 u k ~ O(Nx) preconditioning (diagonal in wavelet space): good initial approximation: solution at previous time step k 1 k ~ O(Nx ) typical realistic beam simulation no preconditioning Ui=0 initial guess 75. preconditioned Ui=0 initial guess 60.7 no preconditioning Ui=Ui 1 initial guess average over step run 4.8
23 Preconditioned Conjugate Gradient (PCG) i convergence rate depends on condition number k u u i k 1 u k ~ O(Nx) preconditioning (diagonal in wavelet space): good initial approximation: solution at previous time step k 1 k ~ O(Nx ) typical realistic beam simulation average over step run no preconditioning Ui=0 initial guess 75. preconditioned Ui=0 initial guess 60.7 no preconditioning Ui=Ui 1 initial guess 4.8 preconditioned Ui=Ui 1 initial guess.4 considerable computational speedup (, Pogorelov & Bohn 007) 3
24 Using PCG in Numerical Simulations Our goal: develop wavelet based Poisson solver which can easily be integrated into existing PIC codes First: test the PCG as a stand alone solver on examples from: beam dynamics galactic dynamics Second: insert the PCG Poisson solver into an existing PIC code (IMPACT T) IMPACT T and run realistic charged particle beam simulations (, Pogorelov & Bohn 007, PR STAB, 10, 03401) compare (conventional FFT based) IMPACT T Vs. IMPACT T with PCG: rms properties level of detail computational speed 4
25 Conventional IMPACT T vs. IMPACT T with PCG Code Comparison: rms Properties Fermilab/NICADD photoinjector grid 1 nc charge rms beam radius rms normalized transverse emittance rms bunch length rms normalized longitudinal emittance Good agreement to a few percent 5
26 Conventional IMPACT T vs. IMPACT T with PCG Code Comparison: Level of Detail & Speed IMPACT T transverse charge distribution for the Fermilab/NICADD photoinjector simulation very non axisymmetric beam grid, N=00000 very good agreement in detail initial distribution z = 1.0 m Speed comparison: IMPACT T w/ PCG ~ 10% faster than the conventional serial IMPACT T IMPACT T w/ PCG z =.0 m z = 4.0 m 6
27 Data Compression with PCG PCG provides excellent compression of data and operators in wavelet space Fermilab/NICADD photoinjector: Real Simulations particle distribution Laplacian operator grid, N=15 000, Nppc=4.58 ~ 3.5% coefficients retained on average grid, N= , Nppc=4.58 ~ 1.75% coefficients retained on average compact storage of beam's distribution history needed for CSR simulations compact storage of beam's potential needed for modeling halo formation 7
28 Ongoing Project: Improving the PCG Solver Currently, we (graduate student Ben Sprague and I) are working on a number of improvements to the wavelet based Poisson equation solver: change from fixed to adaptive grid simplify BC computation (currently a computational bottleneck) further exploit sparsity of operators and data sets use a non standard operator form to better separate scales use more sophisticated wavelet families (biorthogonal, lifted) explore other preconditioners parallelize and optimize Possible future applications of PCG solver: CSR simulations: computation of retarded potentials requires integration over history of the system compactly represented in wavelet space develop a new PIC code to simulate self gravitating systems 8
29 Scaled Gauss Hermite Basis 1 x ψ n x = n H n x e n! π Gauss Hermite orthonormal basis: (solution to quantum harmonic oscillator) Orthonormal: ψ l x ψ m x e x ψ m x dx=δ m dx=δ l m Basis functions: oscillatory exponentially decaying δlm, δm Kronecker delta Infinite expansion (D): 0, 1,, 4, 5 finite: f x, y = al m ψ l x ψ m y l =0 m=0 L M l=0 m=0 l =0 m=0 Scaled and translated version: L f M x y x, y = a l m ψ l x ψ m y α1 α l=0 m=0 1 1 α 1=, α = σx σy σ x x x f x, y dx dy y y { } { } σ y = x = x f x, y dx dy y y {} {} 9
30 Scaled Gauss Hermite Basis N Define collocation points: {γ j } j =0 roots of HL+1(x) { β k }km=0 roots of HM+1(y) L At collocation points: f γ j, β k = a l m ψ l γ j ψ m β k l=0 m =0 Take advantage of the relation for Hermite polynomials: n to obtain al m = H k x H k y k k! k =0 L C jk = = H n 1 x H n y H n x H n 1 y n 1 n! x y 1 C f γ j, β k ψ l γ j ψ m β k j=0 k=0 jk 0 l L, 0 m M M [ψ n γ j ] [ψ m β k ] l =0 γj γ j = x α1 βk β k = y α M L M 0 j L, 0 k M m=0 This formalism is general and can easily be extended to higher dimensions 30
31 Poisson Equation in Scaled Gauss Hermite Basis Poisson equation: x y x y x y ΔΦ x, y = [ x y ]Φ x, y = κf x, y α1 α α1 α α1 α x y Φ x, y = α1 α b l m ψ l x ψ m y l =0 m =0 where blm are given by the difference relation: α 1 l l 1 b n m α m m 1 b l m = κa l m with boundary coefficients: κ a l m, α m m 1 b n m= κ al m, α 1 l l 1 { l=0,1 m, m=0,1 l, } No need to invert the difference equation: compute boundary first, and then work inside computationally simple and efficient 31
32 Simulating Multiparticle Systems with Scaled Gauss Hermite Expansion N body realization of the discrete particle distribution: 1 f x, y = N N δ x x i δ y y i i=1 Expanding f(x,y) in scaled Gauss Hermite basis reduces to the following steps: L 1. tabulate the unchanging quantities: C jk = 1 N 1 y = N x = N xi i=1 N yi α1 = α = i=1 3. evaluate f γ j, β k at the nodes 4. compute coefficients [ψ l γ j ] [ψ m β k ] l =0 pl m j k =. compute x, y, α 1, α : al m = L [ [ N N N x i x i=1 N y i y i =1 1/ ] ] 1/ M m=0 ψ l γ j ψ m β k C jk γj x α1 β β k = k y α γ j = M p l m j k f γ j, β k j=0 k=0 3
33 Simulating Multiparticle Systems with Scaled Gauss Hermite Expansion N body realization of the discrete particle distribution: 1 f x, y = N N δ x x i δ y y i i=1 Expanding f(x,y) in scaled Gauss Hermite basis reduces to the following steps: L 1. tabulate the unchanging quantities: C jk = 1 N 1 y = N x = N xi i=1 N yi α1 = α = i=1 3. evaluate f γ j, β k at the nodes 4. compute coefficients [ψ l γ j ] [ψ m β k ] l =0 pl m j k =. compute x, y, α 1, α : al m = L [ [ N N N x i x i=1 N y i y i =1 1/ ] ] 1/ M M m=0 ψ l γ j ψ m β k C jk γj x α1 β β k = k y α γ j = function estimation from a discrete sample p l m j k f γ j, β k j=0 k=0 33
34 Simulating Multiparticle Systems with Scaled Gauss Hermite Expansion Evaluate f γ j, β k from a discrete sample (nonparametric density estimation) hx hy f γ l, β m = f γ l x, β m y d y d x shifted histogram estimator h x h y hx hy with window [ hx,hx] [ hy,hy] d y d x h x h y Optimal size of the window: 9 [ 1/ 5 hx f ' ' x d x h x ] N 1 /5 Integrated means square error (IMSE) IMSE = h opt = 1 /5 5 4 /5 1/ / 5 hx [ h x f ' ' x d x ] N 4/ 5 ~ N 4 / 5 Other, more sophisticated estimators are available: kernel, adaptive estimators, and others... 34
35 Applying Scaled Gauss Hermite Expansion Future application of scaled Gauss Hermite approximation: D CSR code of Bassi, Ellison, Heinemann and Warnock: particle distribution is sampled by N macroparticles distribution is approximated at each timestep with a cosine expansion beam self forces are computed from the analytic expansion Problems: unphysical wiggles in the tails of the distribution computational speed: each coefficient requires N cosine evaluations Problems resolved (?) with scaled Gauss Hermite: no wiggles basis functions are exponentially decaying computing coefficients scales more favorably and does not involve evaluation of any expensive function (addition & multiplication) 35
36 Applying Scaled Gauss Hermite Expansion Typical simulation: N=106 Wiggles in cosine expan. cosine Reduced by orders of mag. in scaled GH Nx=Ny=4 scaled GH Nx=Ny=10 fraction of negative volume Nx=Ny=16 18x18 grid 36
37 Applying Scaled Gauss Hermite Expansion Scaled Gauss Hermite expansion is computationally appreciably faster: cosine expansion: tcos ~ O(LMNpart) scaled Gauss Hermite: tsgh ~ O((L+M)Npart) + O(LM) ratio tcos/tsgh ~ O(L) + O(Npart/L) (assume L=M) 37
38 Applying Scaled Gauss Hermite Expansion Scaled Gauss Hermite expansion is computationally appreciably faster: cosine expansion: tcos ~ O(LMNpart) scaled Gauss Hermite: tsgh ~ O((L+M)Npart) + O(LM) ratio tcos/tsgh ~ O(L) + O(Npart/L) (assume L=M) ~0 times faster 38
39 Scaled Gauss Hermite Expansion: Loose Ends There are several issues with the scaled Gauss Hermite expansion that we are still exploring/resolving: convergence different estimators for evaluating f γ j, β k optimal number of basis functions (when is more less?) what do we lose by using an analytic expansion? avoid danger of smoothing over physical small scale structures adequate resolution: can this approach resolve physical small scale structure? When these issues are properly addressed, we will have another tool with which to attack CSR and integration over beam's history 39
40 Summary Designed an iterative wavelet based Poisson solver (PCG) wavelet compression and denoising achieves computational speedup preconditioning and sparsity of operators and data in wavelet space reduce CPU load integrated PCG into a PIC code (IMPACT T) for beam dynamics simulations current efforts: adaptive grid, parallelization, optimization future uses: probe usefulness of wavelet methodology in CSR simulations simulate self gravitating systems Developed a scaled Gauss Hermite approximation (still a prototype): efficient representation of particle distribution Poisson equation solved directly at a marginal cost current efforts: resolving issues of convergence, truncation of expansion future uses (?): in Bassi et al.'s D CSR code in Rui Li's D CSR code 40
41 Auxiliary Slides 41
42 Computing BCs: Current Implementation Green's function corresponding to a grounded rectangular pipe V=0 y computational grid EIGENFUNCTIONS sin l x pipe walls sin m y V=0 V=0 EIGENVALUES l =l / a V=0 x pipe walls a b 4 lm z = 0 0 x, y, z sin l x sin m y dx dy ab lm z z m=m / b lm = l m lm lm z = lm z Nx Ny lm x, y, z = l =1 m =1 z sin l x sin m y 4
43 Computing BCs: Adaptive Grid V=0 computational grid V=0 V=0 V=0 43
44 Computing BCs: Adaptive Grid V=0 computational grid V=0 V=0 V=0 44
45 Computing BCs: Adaptive Grid V=0 computational grid V=0 V=0 V=0 45
46 Preconditioner in Wavelet Space
47 Laplacian in Standard Form in Wavelet Space 47
48 Laplacian in Non Standard Form in Wavelet Space 48
49 Wavelet Decomposition The continuous wavelet transform of a function f (t) is s, = f t s, t dt 1 s, t = t s s t mother wavelet with scale and translation dimensions s and respectively 49
50 How Do Wavelets Work? Wavelet analysis (wavelet transform): S signal A approximation D detail Approximation apply low pass filter to Signal and down sample Detail apply high pass filter to Signal and down sample Wavelet synthesis (inverse wavelet transform): up sampling & filtering Complexity: 4MN, M the size of the wavelet, N number of cells Recall: FFT complexity 4N logn 50
51 Wavelets Wavelet transform separates scales superimposed gaussians FFT wavelet transform 51
52 Numerical Noise in PIC Simulations In wavelet space: signal noise few large wavelet coefficients cij many small wavelet coefficients cij Anscombe transformation Poissonian noise Denoising by wavelet thresholding: Gaussian noise if cij < T, set to cij =0 (choose threshold T carefully!) A great deal of study has been devoted to estimating optimal T Nppc=5 T = log N grid DF(T) (σ was estimated earlier) Ni=64 NGP Nppc=5 DF(T) CIC Ni=64 w/ Anscombe transformation w/out Anscombe transformation JLab Beam Terzi Seminar, April 008 & Bohn 007, PR STAB, 10, ć, Pogorelov
53 3D Plummer sphere on a 3x3x3 grid BCs: open in all directions (analytically specified) 53
54 3D 'fuzzy cigar' on a 3x3x3 grid BCs: grounded rectangular pipe (V=0 on the sides), open in z direction (Green's function solution) 54
55 D EXAMPLE: MICROBUNCHED DENSITY DISTRIBUTION (with U=0 starting guess) BCs: open in all directions (analytically specified) 55
56 IMPACT T Vs. IMPACT T with PCG: rms properties Fermilab/NICADD photoinjector (3 3 3 grid) 1 nc charge bunch compressor at z=5.65m from cathode rms beam radius rms normalized transverse emittance rms bunch length rms normalized longitudinal emittance 56
57 IMPACT T Vs. IMPACT T with PCG: level of detail IMPACT T transverse charge distribution for the Fermilab/NICADD photoinjector simulation 5 beamlets (masked photocathode) grid, N=00000 very good agreement in detail IMPACT T w/ PCG initial distribution z =.0 m z = 4.0 m z = 6.0 m 57
58 IMPACT T Vs. IMPACT T with PCG: level of detail Impact IMPACT T IMPACT Tresidue w/ PCG Impact with CPCG transverse charge distribution for the AES/JLab low charge photoinjector simulation very non axisymmetric beam grid, N=00000 very good agreement in detail initial distribution z = 0.5 m z = 1.9 m z =.5 m 58
59 IMPACT T Vs. IMPACT T with PCG: level of detail Impact IMPACT T IMPACT Tresidue w/ PCG Impact with CPCG transverse charge distribution for the AES/JLab low charge photoinjector simulation very non axisymmetric beam grid, N=00000 very good agreement in detail initial distribution z = 0.5 m z = 1.9 m z =.5 m 59
60 Scaled Gauss Hermite Basis 1 x ψ n x = n H n x e n! π Gauss Hermite orthonormal basis: (solution to quantum harmonic oscillator) Hermite polynomials Hn(x) relations: Orthonormal: ψ l x ψ m x e x δlm, δm Basis functions: oscillatory exponentially decaying 0, 1,, 4, 5 Infinite expansion (D): f x, y = al m ψ l x ψ m y Scaled and translated version: L M x y f x, y = a l m ψ l x ψ m y α1 α l=0 m=0 α 1= 1 1, α = σx σy σ x finite: l =0 m=0 Kronecker delta H n 1 x = x H n x n H n 1 x H n ' x = x H n 1 x ψ m x dx=δ m dx=δ l m L M l=0 m=0 l =0 m=0 x x f x, y dx dy y y { } { } σ y = x = y {} x f x, y dx dy y {} 60
Particle In Cell Beam Dynamics Simulations With a Wavelet based Poisson Solver
Particle In Cell Beam Dynamics Simulations With a Wavelet based Poisson Solver Balša Terzić Beam Physics and Astrophysics Group, Department of Physics Northern Illinois University in collaboration with
More informationApplication of Wavelets to N body Particle In Cell Simulations
Application of Wavelets to N body Particle In Cell Simulations Balša Terzić (NIU) Northern Illinois University, Department of Physics April 7, 2006 Motivation Gaining insight into the dynamics of multi
More informationProgress Report: New Particle-In-Cell Code For Numerical Simulations of Coherent Synchrotron Radiation
Progress Report: New Particle-In-Cell Code For Numerical Simulations of Coherent Synchrotron Radiation Balša Terzić and Rui Li CASA, Jefferson Lab CASA Group Meeting May 26, 2010 Notice: Authored by Jefferson
More informationParticle-in-cell beam dynamics simulations with a wavelet-based Poisson solver
PHYSICAL REVIEW SPECIAL TOPICS - ACCELERATORS AND BEAMS, 3 (7) Particle-in-cell beam dynamics simulations with a wavelet-based Poisson solver Balša Terzić, Ilya V. Pogorelov, and Courtlandt L. Bohn Beam
More informationParticle-in-Cell Beam Dynamics Simulations With A Wavelet-Based Poisson Solver
Particle-in-Cell Beam Dynamics Simulations With A Wavelet-Based Poisson Solver Balša Terzić, Ilya V. Pogorelov, Courtlandt L. Bohn Beam Physics and Astrophysics Group, Northern Illinois University, DeKalb,
More informationWavelet-Based Poisson Solver for Use in Particle-in-Cell Simulations
Wavelet-Based Poisson Solver for Use in Particle-in-Cell Simulations Balša Terzić Department of Physics, Northern Illinois University, DeKalb, Illinois 65 Ilya V. Pogorelov Accelerator and Fusion Research
More informationFast Space Charge Calculations with a Multigrid Poisson Solver & Applications
Fast Space Charge Calculations with a Multigrid Poisson Solver & Applications Gisela Pöplau Ursula van Rienen Rostock University DESY, Hamburg, April 26, 2005 Bas van der Geer Marieke de Loos Pulsar Physics
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationABSTRACT. Title: Wavelet-space solution of the Poisson equation: An algorithm for use in particle-incell simulations
ABSTRACT Name: Benjamin Sprague Department: Physics Title: Wavelet-space solution of the Poisson equation: An algorithm for use in particle-incell simulations Major: Physics Degree: Master of Science Approved
More informationEFFICIENT 3D POISSON SOLVERS FOR SPACE-CHARGE SIMULATION
Abstract EFFICIENT 3D POISSON SOLVERS FOR SPACE-CHARGE SIMULATION Three-dimensional Poisson solver plays an important role in the self-consistent space-charge simulation. In this paper, we present several
More informationSimulations of the Microbunching Instability in FEL Beam Delivery Systems
Simulations of the Microbunching Instability in FEL Beam Delivery Systems Ilya Pogorelov Tech-X Corporation Workshop on High Average Power & High Brightness Beams UCLA, January 2009 Outline The setting:
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More information3D Space Charge Routines: The Software Package MOEVE and FFT Compared
3D Space Charge Routines: The Software Package MOEVE and FFT Compared Gisela Pöplau DESY, Hamburg, December 4, 2007 Overview Algorithms for 3D space charge calculations Properties of FFT and iterative
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing (Wavelet Transforms) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids
More informationCSR calculation by paraxial approximation
CSR calculation by paraxial approximation Tomonori Agoh (KEK) Seminar at Stanford Linear Accelerator Center, March 3, 2006 Short Bunch Introduction Colliders for high luminosity ERL for short duration
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
More informationModule 4 MULTI- RESOLUTION ANALYSIS. Version 2 ECE IIT, Kharagpur
Module 4 MULTI- RESOLUTION ANALYSIS Lesson Theory of Wavelets Instructional Objectives At the end of this lesson, the students should be able to:. Explain the space-frequency localization problem in sinusoidal
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Multi-dimensional signal processing Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione dei Segnali Multi-dimensionali e
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationAdaptive Nonparametric Density Estimators
Adaptive Nonparametric Density Estimators by Alan J. Izenman Introduction Theoretical results and practical application of histograms as density estimators usually assume a fixed-partition approach, where
More informationSparse Matrix Techniques for MCAO
Sparse Matrix Techniques for MCAO Luc Gilles lgilles@mtu.edu Michigan Technological University, ECE Department Brent Ellerbroek bellerbroek@gemini.edu Gemini Observatory Curt Vogel vogel@math.montana.edu
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationWavelets and Signal Processing
Wavelets and Signal Processing John E. Gilbert Mathematics in Science Lecture April 30, 2002. Publicity Mathematics In Science* A LECTURE SERIES FOR UNDERGRADUATES Wavelets Professor John Gilbert Mathematics
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationRepresentation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices
Representation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices Andreas Rhodin, Harald Anlauf German Weather Service (DWD) Workshop on Flow-dependent aspects of data assimilation,
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More information1 The Continuous Wavelet Transform The continuous wavelet transform (CWT) Discretisation of the CWT... 2
Contents 1 The Continuous Wavelet Transform 1 1.1 The continuous wavelet transform (CWT)............. 1 1. Discretisation of the CWT...................... Stationary wavelet transform or redundant wavelet
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More informationSimulations of single bunch collective effects using HEADTAIL
Simulations of single bunch collective effects using HEADTAIL G. Rumolo, in collaboration with E. Benedetto, O. Boine-Frankenheim, G. Franchetti, E. Métral, F. Zimmermann ICAP, Chamonix, 02.10.2006 Giovanni
More informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationMultiresolution Analysis
Multiresolution Analysis DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Frames Short-time Fourier transform
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationMachine vision. Summary # 4. The mask for Laplacian is given
1 Machine vision Summary # 4 The mask for Laplacian is given L = 0 1 0 1 4 1 (6) 0 1 0 Another Laplacian mask that gives more importance to the center element is L = 1 1 1 1 8 1 (7) 1 1 1 Note that the
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationSolving Boundary Value Problems (with Gaussians)
What is a boundary value problem? Solving Boundary Value Problems (with Gaussians) Definition A differential equation with constraints on the boundary Michael McCourt Division Argonne National Laboratory
More informationA Vlasov-Maxwell Solver to Study Microbunching Instability in the First Bunch Compressor System
Microbunching in the FERMI@ELETTRA First Bunch Compressor System / Gabriele Bassi Page 1 A Vlasov-Maxwell Solver to Study Microbunching Instability in the FERMI@ELETTRA First Bunch Compressor System Gabriele
More informationNUMERICAL CALCULATION OF RANDOM MATRIX DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS. Sheehan Olver NA Group, Oxford
NUMERICAL CALCULATION OF RANDOM MATRIX DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS Sheehan Olver NA Group, Oxford We are interested in numerically computing eigenvalue statistics of the GUE ensembles, i.e.,
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationDiscrete Orthogonal Harmonic Transforms
Discrete Orthogonal Harmonic Transforms Speaker: Chun-Lin Liu, Advisor: Soo-Chang Pei Ph. D Image Processing Laboratory, EEII 530, Graduate Institute of Communication Engineering, National Taiwan University.
More informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationMATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials
MATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials. Hermite Equation In the study of the eigenvalue problem of the Hamiltonian for the quantum harmonic oscillator we have
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationOne-dimensional Schrödinger equation
Chapter 1 One-dimensional Schrödinger equation In this chapter we will start from the harmonic oscillator to introduce a general numerical methodology to solve the one-dimensional, time-independent Schrödinger
More informationMachine vision, spring 2018 Summary 4
Machine vision Summary # 4 The mask for Laplacian is given L = 4 (6) Another Laplacian mask that gives more importance to the center element is given by L = 8 (7) Note that the sum of the elements in the
More informationNumerical Linear Algebra Primer. Ryan Tibshirani Convex Optimization /36-725
Numerical Linear Algebra Primer Ryan Tibshirani Convex Optimization 10-725/36-725 Last time: proximal gradient descent Consider the problem min g(x) + h(x) with g, h convex, g differentiable, and h simple
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationExternal injection of electron bunches into plasma wakefields
External injection of electron bunches into plasma wakefields Studies on emittance growth and bunch compression p x x Timon Mehrling, Julia Grebenyuk and Jens Osterhoff FLA, Plasma Acceleration Group (http://plasma.desy.de)
More informationModelling the Microbunching Instability for Bunch Compressor Systems
Modelling the Microbunching Instability for Bunch Compressor Systems Gabriele Bassi Page Modelling the Microbunching Instability for Bunch Compressor Systems Gabriele Bassi Department of Physics, University
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationDESY/TEMF Meeting - Status 2011
Erion Gjonaj Technische Universität Darmstadt, Computational Electromagetics Laboratory (TEMF) Schlossgartenstr. 8 64289 Darmstadt, Germany DESY/TEMF Meeting - Status 211 DESY, Hamburg, 16.12.211 Contents
More informationTime resolved transverse and longitudinal phase space measurements at the high brightness photo injector PITZ
Time resolved transverse and longitudinal phase space measurements at the high brightness photo injector PITZ 1. Motivation 2. Transverse deflecting structure 3. Longitudinal phase space tomography 4.
More informationMath 115 ( ) Yum-Tong Siu 1. Derivation of the Poisson Kernel by Fourier Series and Convolution
Math 5 (006-007 Yum-Tong Siu. Derivation of the Poisson Kernel by Fourier Series and Convolution We are going to give a second derivation of the Poisson kernel by using Fourier series and convolution.
More informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationNonparametric Methods
Nonparametric Methods Michael R. Roberts Department of Finance The Wharton School University of Pennsylvania July 28, 2009 Michael R. Roberts Nonparametric Methods 1/42 Overview Great for data analysis
More informationCHAPTER 6 Quantum Mechanics II
CHAPTER 6 Quantum Mechanics II 6.1 6.2 6.3 6.4 6.5 6.6 6.7 The Schrödinger Wave Equation Expectation Values Infinite Square-Well Potential Finite Square-Well Potential Three-Dimensional Infinite-Potential
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationIterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems Luca Bergamaschi e-mail: berga@dmsa.unipd.it - http://www.dmsa.unipd.it/ berga Department of Mathematical Methods and Models for Scientific Applications University
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationarxiv: v1 [physics.acc-ph] 21 Nov 2011
arxiv:1111.4971v1 [physics.acc-ph] 21 Nov 2011 On FFT based convolutions and correlations, with application to solving Poisson s equation in an open rectangular pipe Abstract Robert D. Ryne Lawrence Berkeley
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationThe Quantum Harmonic Oscillator
The Classical Analysis Recall the mass-spring system where we first introduced unforced harmonic motion. The DE that describes the system is: where: Note that throughout this discussion the variables =
More informationTwo-Dimensional Signal Processing and Image De-noising
Two-Dimensional Signal Processing and Image De-noising Alec Koppel, Mark Eisen, Alejandro Ribeiro March 12, 2018 Until now, we considered (one-dimensional) discrete signals of the form x : [0, N 1] C of
More information1 Infinite-Dimensional Vector Spaces
Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators
More information(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y
() Consider A = { q Q : q 2 2} as a subset of the metric space (Q, d), where d(x, y) = x y. Then A is A) closed but not open in Q B) open but not closed in Q C) neither open nor closed in Q D) both open
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationmodeling of space charge effects and CSR in bunch compression systems
modeling of space charge effects and CSR in bunch compression systems SC and CSR effects are crucial for the simulation of BC systems CSR and related effects are challenging for EM field calculation non-csr
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationInterpolation & Polynomial Approximation. Hermite Interpolation I
Interpolation & Polynomial Approximation Hermite Interpolation I Numerical Analysis (th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationModeling of Space Charge Effects and Coherent Synchrotron Radiation in Bunch Compression Systems. Martin Dohlus DESY, Hamburg
Modeling of Space Charge Effects and Coherent Synchrotron Radiation in Bunch Compression Systems Martin Dohlus DESY, Hamburg SC and CSR effects are crucial for design & simulation of BC systems CSR and
More informationPDF Estimation via Characteristic Function and an Orthonormal Basis Set
PDF Estimation via Characteristic Function and an Orthonormal Basis Set ROY M. HOWARD School of Electrical Engineering and Computing Curtin University of Technology GPO Bo U987, Perth 6845. AUSTRALIA r.howard@echange.curtin.edu.au
More informationWavelets Marialuce Graziadei
Wavelets Marialuce Graziadei 1. A brief summary 2. Vanishing moments 3. 2D-wavelets 4. Compression 5. De-noising 1 1. A brief summary φ(t): scaling function. For φ the 2-scale relation hold φ(t) = p k
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationChapter 7 Wavelets and Multiresolution Processing. Subband coding Quadrature mirror filtering Pyramid image processing
Chapter 7 Wavelets and Multiresolution Processing Wavelet transform vs Fourier transform Basis functions are small waves called wavelet with different frequency and limited duration Multiresolution theory:
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More information