Parameter Identification in Partial Differential Equations
|
|
- Gertrude Brooks
- 6 years ago
- Views:
Transcription
1 Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation of a function f, Martin Burger Johannes Kepler Universität Linz Find its derivative Winter School Inverse Problems, Geilo 2 Differentiation of data Pointwise noise at measurement points identically normally distributed with mean 0 and variance δ. Differentiation of data Example: Law of large numbers yields Noise satisfies But error in derivative is large Winter School Inverse Problems, Geilo 3 Winter School Inverse Problems, Geilo 4
2 Differentiation of data Conclusion 1 Error: Arbitrarily large since k can be arbitrarily large Winter School Inverse Problems, Geilo 5 Winter School Inverse Problems, Geilo 6 Differentiation of data Additional information: Solution is smooth (e.g. twice differentiable) Regularization: e.g. smoothing by elliptic PDE Differentiaton of data Detailed estimate: Lecture notes, p. 7 Variational principle for elliptic PDEs: minimize Winter School Inverse Problems, Geilo 7 Winter School Inverse Problems, Geilo 8
3 Computerized Tomography (A nice link between Austria & Norway) Problem: reconstruct a spatial density f in a domain D from measurements of X-rays traveling through the domain. Computerized Tomography X-ray beam has intensity I Model: decay of itensity in a small distance t proportional to I, f, and distance X-rays travel along rays, parametrized by distance s from origin and its unit normal vector ω Limit t to zero Winter School Inverse Problems, Geilo 9 Winter School Inverse Problems, Geilo 10 Computerized Tomography Measurement: intensity at emitter and detector Computerized Tomography Integrate ODE Parameter identification problem for system of ordinary differential equations (first-order) Leads to inversion of the Radon transform Overdetermined for given f since initial and final value are known Winter School Inverse Problems, Geilo 11 Winter School Inverse Problems, Geilo 12
4 Computerized Tomography Radon Transform Computerized Tomography Radial symmetry, Exact inversion formula by Johann Radon SVD computed by McCormick et. al. in 1960s, Nobel Prize for Medicine in the 1970 Winter School Inverse Problems, Geilo 13 Winter School Inverse Problems, Geilo 14 Computerized Tomography Radial symmetry, Groundwater Filtration Identification of diffusivity of sediments from an observation of the piezometric head (describing the flow) Abel Integral equation! SVs decay with half speed as for differentiation Knowledge of diffusivity allows to draw conclusions about the structure of the groundwater Winter School Inverse Problems, Geilo 15 Winter School Inverse Problems, Geilo 16
5 Groundwater Filtration Mathematical model Groundwater Filtration Instead of exact data, we only know noisy measurement Diffusivity a, piezometric head u (measured), density of water sources and sinks f Typically the noise n(x) comes from identically normally distributed random variables From a law of large numbers this gives an estimate + Appropriate boundary conditions (e.g. u=0) Winter School Inverse Problems, Geilo 17 Winter School Inverse Problems, Geilo 18 Groundwater Filtration Consider 1D version of inverse problem, for simplicity with We can integrate the equation to obtain Groundwater Filtration Detailed estimate: Lecture notes, p. 9 Hence, if the derivative of u does not vanish, a is determined uniquely (Identifiability) Winter School Inverse Problems, Geilo 19 Winter School Inverse Problems, Geilo 20
6 Groundwater Filtration Some results (from Hanke, 1995) Groundwater Filtration Reconstructions for two different noise levels (1% and 0.1%) Exact phantom Winter School Inverse Problems, Geilo 21 Winter School Inverse Problems, Geilo 22 Groundwater Filtration Reconstruction along the diagonal Groundwater Filtration Reconstructions of piecewise constant parameter, no noise (nonuniqueness curve), fine grid Winter School Inverse Problems, Geilo 23 Winter School Inverse Problems, Geilo 24
7 Groundwater Filtration Reconstructions of piecewise constant parameter, no noise (nonuniqueness curve), hierarchical grids Electrical Impedance Tomography Application in Medical Imaging Set of electrodes placed around human chest, Response to various electrical impulses measured Picture from RPI-EIT Project Winter School Inverse Problems, Geilo 25 Winter School Inverse Problems, Geilo 26 Electrical Impedance Tomography Electrical Impedance Tomography Picture from RPI-EIT Project Pictures from RPI-EIT Project Winter School Inverse Problems, Geilo 27 Winter School Inverse Problems, Geilo 28
8 Electrical Impedance Tomography Mathematical model: Maxwell equations, reduced to potential equation Electrical Impedance Tomography Measurement: current density on the boundary for different voltage patterns u is electrical potential, f applied voltage pattern a denotes the (unknown) conductivity Idealized mathematical model: Dirichlet-Neumann map Winter School Inverse Problems, Geilo 29 Winter School Inverse Problems, Geilo 30 Electrical Impedance Tomography INVERSE CONDUCTIVITY PROBLEM In practice finite number of measurements Electrical Impedance Tomography To compute output, we have to solve solutions of N partial differential equations N typically very large Problem of extremely large scale Winter School Inverse Problems, Geilo 31 Winter School Inverse Problems, Geilo 32
9 Electrical Impedance Tomography Interesting special case: piecewise constant conductivities Semiconductor Devices Related problem to impedance tomography appears for semiconductor devices: inverse dopant profiling Real interest is the subset where a takes a value different from a 1 (e.g. a tumour) Winter School Inverse Problems, Geilo 33 Winter School Inverse Problems, Geilo 34 Inverse Dopant Profiling Identify the device doping profile from measurements Current-Voltage map Analogous definitions of voltage and current, but more complicated mathematical model Mathematical Model Stationary Drift Diffusion Model: PDE system for potential V, electron density n and hole density p in Ω (subset of R 2 ) Doping Profile C(x) enters as source term Winter School Inverse Problems, Geilo 35 Winter School Inverse Problems, Geilo 36
10 Boundary Conditions Boundary of Ω : homogeneous Neumann boundary conditions on Γ N and Device Characteristics Measured on a contact Γ 0 on Γ D : Outflow current density on Dirichlet boundary Γ D (Ohmic Contacts) Winter School Inverse Problems, Geilo 37 Winter School Inverse Problems, Geilo 38 Numerical Tests Test for a P-N Diode Numerical Tests Different Voltage Sources Real Doping Profile Initial Guess Winter School Inverse Problems, Geilo 39 Winter School Inverse Problems, Geilo 40
11 Numerical Tests Reconstructions with first source Numerical Tests Reconstructions with second source Winter School Inverse Problems, Geilo 41 Winter School Inverse Problems, Geilo 42 Polymer Crystallization Polymer Crystallization Mesoscale model for polymer crystallization ξ = degree of crystallinity u,v = surface densities T = temperature G = growth rate N = nucleation rate Winter School Inverse Problems, Geilo 43 Winter School Inverse Problems, Geilo 44
12 Polymer Crystallization Source term Polymer Crystallization Traditional way of determining nucleation rate as function of temperature: - Make separate experiment for each value of T - Count (by eyes!!!) the final number of crystals (typically >> 10 6 ). - Divide by volume Extremely expensive, extremely timeconsuming! Winter School Inverse Problems, Geilo 45 Winter School Inverse Problems, Geilo 46 Polymer Crystallization Idea: determine nucleation rate as function of temperature by single nonisothermal experiment Polymer Crystallization Results Measured data: temperature T at the boundary of the sample Degree of crystallinity at final time Winter School Inverse Problems, Geilo 47 Winter School Inverse Problems, Geilo 48
13 Polymer Crystallization Results Polymer Crystallization Results Winter School Inverse Problems, Geilo 49 Winter School Inverse Problems, Geilo 50 Polymer Crystallization Reconstruction error vs noise level Polymer Crystallization Reconstruction error vs noise level Winter School Inverse Problems, Geilo 51 Winter School Inverse Problems, Geilo 52
14 Polymer Crystallization Iteration numbers Polymer Crystallization Iteration numbers Winter School Inverse Problems, Geilo 53 Winter School Inverse Problems, Geilo 54 Iterative Regularization For nonlinear inverse problems, iterative regularization is of even higher importance as for linear ones We need to perform iteration to minimize nonlinear Tikhonov functional anyway Iterative Regularization Gradient of the least-squares functional Perform simple gradient descent: Start from least-squares functional Landweber iteration Winter School Inverse Problems, Geilo 55 Winter School Inverse Problems, Geilo 56
15 Iterative Regularization Where is the regularization parameter??!! Iterative Regularization Regularization by appropriate early stopping, e.g. discrepancy principle Each step of Landweber iteration is well-posed, so we obtain iterative regularization method Winter School Inverse Problems, Geilo 57 Winter School Inverse Problems, Geilo 58 Iterative Regularization Analysis of Landweber iteration: Hanke-Neubauer- Scherzer 1994, Scherzer 1995 Usual semi-convergence properties as for linear problems, but restriction on the nonlinearity of the problem (replacing regularity of ) Iterative Regularization Landweber iteration is forward Euler time discretization (with time step τ of the flow) Asymptotical regularization (Tautenhahn 1994) Iteration Winter School Inverse Problems, Geilo 59 Winter School Inverse Problems, Geilo 60
16 Iterative Regularization We can consider other time discretizations: Backward Euler Iterative Regularization We can consider other time discretizations: Semi-Implicit Euler Iterated Tikhonov Regularization (Hanke- Groetsch 1999, Groetsch-Scherzer 1999) Levenberg-Marquardt (Hanke 1995) Arbitrary Runge-Kutta methods possibly, even inconsistent ones (Rieder, 2005) Winter School Inverse Problems, Geilo 61 Winter School Inverse Problems, Geilo 62 Iterative Regularization Levenberg-Marquardt can be rewritten to Iterative Regularization General approach for Newton-type methods: Apply linear regularization to Tikhonov regularization of the linear Newtonequation Iteratively regularized Gauss-Newton (Kaltenbacher-Neubauer-Scherzer 1996, Kaltenbacher 1997) Newton-CG (Hanke, 1998) Winter School Inverse Problems, Geilo 63 Winter School Inverse Problems, Geilo 64
17 Iterative Regularization Total Variation Denoising Newton-Landweber (Kaltenbacher 1999) Broyden (Kaltenbacher, 1997) Multigrid / Discretization (Kaltenbacher, ) In particular for parameter identification: Levenberg-Marquardt SQP (B-Mühlhuber 2002) Newton-Kaczmarcz (B-Kaltenbacher 2004) Winter School Inverse Problems, Geilo 65 Winter School Inverse Problems, Geilo 66 Total Variation Denoising TV-Images Basic problem in denoising: given noisy version g of image f, find optimal approximation u of f Basic problem in deblurring: given noisy version g of Kf, find optimal approximation u of f Winter School Inverse Problems, Geilo 67 Winter School Inverse Problems, Geilo 68
18 TV-Images What is optimal? Satisfy other criterion, e.g., minimization ROF-Model What is a suitable class of approximations? Discrepancy principle: allow images with with σ being the noise level: over a class of approximations Winter School Inverse Problems, Geilo 69 Winter School Inverse Problems, Geilo 70 ROF-Model Rudin-Osher-Fatemi 1989 (ROF): Minimize total variation ROF-Model Deblurring: Minimize total variation subject to subject to Winter School Inverse Problems, Geilo 71 Winter School Inverse Problems, Geilo 72
19 ROF-Model Acar-Vogel 1994: ROF is regularization method (wellposedness + convergence as σ to zero) Chambolle-Lions 1997: Solution unique, there exists Lagrange multiplier λ, problem equivalent to ROF-Model Acar-Vogel 1994: ROF is regularization method (wellposedness + convergence as σ to zero) Chambolle-Lions 1997: Solution unique, there exists Lagrange multiplier λ, problem equivalent to Winter School Inverse Problems, Geilo 73 Winter School Inverse Problems, Geilo 74 Total Variation Rigorous definition of total variation Total Variation BV includes discontinuous functions: 1D example Winter School Inverse Problems, Geilo 75 Winter School Inverse Problems, Geilo 76
20 Total Variation 1D example Total Variation Formal optimality for TV-Denoising Term on the right-hand side corresponds to mean curvature of level sets of u Winter School Inverse Problems, Geilo 77 Winter School Inverse Problems, Geilo 78 Duality Consider TV-regularization Duality Exchange inf and sup Solve inner inf problem for u Winter School Inverse Problems, Geilo 79 Winter School Inverse Problems, Geilo 80
21 Duality Remains sup problem for g Duality Introduce Rewrite equivalently Dual TV problem (Chambolle 2003) Dual space of BV Winter School Inverse Problems, Geilo 81 Winter School Inverse Problems, Geilo 82 Stair-Casing Lecture notes p.20 ROF-Model Meyer 2001 (and others before): ROF has a systematic error, since This means that jumps in the reconstructed image are smaller than jumps in the true image Winter School Inverse Problems, Geilo 83 Winter School Inverse Problems, Geilo 84
22 Noise Decomposition How to resolve the systematic error? Take some λ, and minimize TV-functional Noise Decomposition Minimize to obtain. Then take residual ( noise ) to obtain next iterate The second step may increase total variation! Winter School Inverse Problems, Geilo 85 Winter School Inverse Problems, Geilo 86 Iterated Total Variation Some Results Obtain new iterate by minimizing then decompose noise Need not change the code, just the data! Winter School Inverse Problems, Geilo 87 Winter School Inverse Problems, Geilo 88
23 Some Results Some Results Winter School Inverse Problems, Geilo 89 Winter School Inverse Problems, Geilo 90 Some Results Some Results Winter School Inverse Problems, Geilo 91 Winter School Inverse Problems, Geilo 92
24 Some Results Convergence Analysis So why does that work??? Expand fitting term, throw away constant terms: Winter School Inverse Problems, Geilo 93 Winter School Inverse Problems, Geilo 94 Convergence Analysis Optimality condition implies Convergence Analysis Hence, iterated TV is proximal point algorithm Write equivalent optimization problem Winter School Inverse Problems, Geilo 95 Winter School Inverse Problems, Geilo 96
25 Convergence Analysis The second term is called Bregman distance Convergence Analysis Example: quadratic case with yields Winter School Inverse Problems, Geilo 97 Iterated Tikhonov / Levenberg-Marquardt Method! Winter School Inverse Problems, Geilo 98 Inverse TV Flow Does the reconstruction depend on λ? To understand the dependence on the parameter, consider limit λ 0. Let λ N = T/N, set t k = k λ N. We compute the solution of the iterated TV flow with this specific parameter and set Inverse TV Flow We obtain a limiting flow, This flow is uniquely defined Formally, Linear interpolation for other times. Winter School Inverse Problems, Geilo 99 Winter School Inverse Problems, Geilo 100
26 Inverse TV Flow The inverse TV flow can be considered as an inverse scale space method (Scherzer-Weickert 99) Start with zero image, and gradually insert smaller and smaller scales. Noise should return in the end, this happened in experiments Inverse TV Flow The parameter λ can be interpreted as a time step for an implicit discretization of the inverse TV flow Should work well and be independent of the parameter, as long as λ is small enough. This can be observed in experiments, too. Winter School Inverse Problems, Geilo 101 Winter School Inverse Problems, Geilo 102 The algorithm can immediately be generalized to arbitrary linear and even nonlinear inverse problems: Find next iterate as minimizer of Winter School Inverse Problems, Geilo 103 Example: estimate diffusion coefficient q in from measurement of u Operator F: q u(q) Application: find material layers in groundwater modeling (piecewise constant) Winter School Inverse Problems, Geilo 104
27 Exact data Iterate 1 Exact data Iterate 2 Winter School Inverse Problems, Geilo 105 Winter School Inverse Problems, Geilo 106 Exact data Iterate 3 Exact data Iterate 4 Winter School Inverse Problems, Geilo 107 Winter School Inverse Problems, Geilo 108
28 Exact data Iterate 5 Exact data Iterate 6 Winter School Inverse Problems, Geilo 109 Winter School Inverse Problems, Geilo 110 Exact data Iterate 10 Exact data Iterate 20 Winter School Inverse Problems, Geilo 111 Winter School Inverse Problems, Geilo 112
29 Exact data Iterate 50 Exact data Iterate 100 Winter School Inverse Problems, Geilo 113 Winter School Inverse Problems, Geilo 114 Exact data Iterate 200 1% noise Iterate 1 Winter School Inverse Problems, Geilo 115 Winter School Inverse Problems, Geilo 116
30 1% noise Iterate 2 1% noise Iterate 3 Winter School Inverse Problems, Geilo 117 Winter School Inverse Problems, Geilo 118 1% noise Iterate 4 1% noise Iterate 5 Winter School Inverse Problems, Geilo 119 Winter School Inverse Problems, Geilo 120
31 1% noise Iterate 10 1% noise Iterate 50 Winter School Inverse Problems, Geilo 121 Winter School Inverse Problems, Geilo % noise Iterate 1 10% noise Iterate 2 Winter School Inverse Problems, Geilo 123 Winter School Inverse Problems, Geilo 124
32 10% noise Iterate 3 10% noise Iterate 4 Winter School Inverse Problems, Geilo 125 Winter School Inverse Problems, Geilo % noise Iterate 5 10% noise Iterate 10 Winter School Inverse Problems, Geilo 127 Winter School Inverse Problems, Geilo 128
33 Parameter Identification Implementation Aspects Landweber iteration: Discretize parameter Discretize state + state equation Compute gradient by adjoint method Try to use same kind of discretization for adjoint as for state equation ( Discretized Adjoint = Adjoint of Discretized Problem ) For multiple state equations compute gradients corresponding to each state equation immediately Winter School Inverse Problems, Geilo 129 Winter School Inverse Problems, Geilo 130 Implementation Aspects Landweber iteration with Black Box Solver: Discretize parameter Discretize state + state equation (determined by solver) If adjoint method not possible, try to compute gradients by finite differencing For multiple state equations compute gradients corresponding to each state equation immediately Implementation Aspects Quasi-Newton Methods (BFGS): Discretize parameter Discretize state + state equation (determined by solver) Compute gradients by adjoint method, use adjoint to update quasi-newton matrix Solve linear system for update Winter School Inverse Problems, Geilo 131 Winter School Inverse Problems, Geilo 132
34 Implementation Aspects Newton Methods (LM / IRGN / NCG): Discretize parameter Discretize state + state equation Never compute Newton-Matrix (NOT SPARSE)!! Implementation Aspects One shot methods: Discretize parameter + state + state equation + adjoint + adjoint equation Solve the full KKT system simultaneously, e.g. Try to use iterative method for computing the Newton step, only implement application of and its adjoint (two PDEs) Winter School Inverse Problems, Geilo 133 Winter School Inverse Problems, Geilo 134 Implementation Aspects One shot methods: KKT system is sparse and symmetric KKT system is not positive definite (no CG!) Use BiCGStab, MINRES, QMR or GMRES Apply appropriate preconditioning, note that there is a small parameter (α), acting like singular perturbation Preconditioning strategies: Battermann-Sachs 2000, Battermann-Heinckenschloss 2001, Ascher-Haber , B-Mühlhuber 2002, Griewank 2005 Kaczmarcz Methods For problems like impedance tomography, the operators are of the form Idea: perform multiplicative splitting (like in Gauss- Seidel), cyclic iteration over the N subproblems Winter School Inverse Problems, Geilo 135 Winter School Inverse Problems, Geilo 136
35 Kaczmarcz Methods Subproblem j: Landweber-Kaczmarcz (Natterer 1996, Kowar-Scherzer 2003) Kaczmarcz Methods Levenberg-Marquardt Kaczmarcz (B-Kaltenbacher 2004) State equation: Linear KKT: Winter School Inverse Problems, Geilo 137 Winter School Inverse Problems, Geilo 138 Kaczmarcz Methods Simple model problem: identify q in Numerical Solution Newton-Kaczmarz approach is equivalent to minimizing from measurements for in each step, subject to for different sources Winter School Inverse Problems, Geilo 139 Winter School Inverse Problems, Geilo 140
36 Numerical Solution Since we need another function to realize the H -1/2 -norm, corresponding KKT system is 4 x 4 Numerical Solution With boundary conditions and Winter School Inverse Problems, Geilo 141 Winter School Inverse Problems, Geilo 142 Numerical Results Numerical experiment with N=20 localized sources g j (5 on each side) Numerical Results Iterates 1, 2, 3, 4, 5, 6 Exact solution Constant initial value Winter School Inverse Problems, Geilo 143 Winter School Inverse Problems, Geilo 144
37 Numerical Results Iterates 10, 15, 20, 40, 60, 80 (1/2,3/4,1,2,3,4) Numerical Results Residuals and Error Winter School Inverse Problems, Geilo 145 Winter School Inverse Problems, Geilo 146 Numerical Results Numerical experiment with N=20 localized sources g j (5 on each side) Numerical Results Iterates 1, 6, 11, 16, 21, 41 Exact solution Constant initial value Winter School Inverse Problems, Geilo 147 Winter School Inverse Problems, Geilo 148
38 Numerical Results Iterates 80, 120, 160, 240, 320, 400 (4,6,8,12,16,20) Numerical Results Residuals and Error Winter School Inverse Problems, Geilo 149 Winter School Inverse Problems, Geilo 150
Iterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationFunctionalanalytic tools and nonlinear equations
Functionalanalytic tools and nonlinear equations Johann Baumeister Goethe University, Frankfurt, Germany Rio de Janeiro / October 2017 Outline Fréchet differentiability of the (PtS) mapping. Nonlinear
More informationRegularization in Banach Space
Regularization in Banach Space Barbara Kaltenbacher, Alpen-Adria-Universität Klagenfurt joint work with Uno Hämarik, University of Tartu Bernd Hofmann, Technical University of Chemnitz Urve Kangro, University
More informationConvergence rates of convex variational regularization
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 20 (2004) 1411 1421 INVERSE PROBLEMS PII: S0266-5611(04)77358-X Convergence rates of convex variational regularization Martin Burger 1 and Stanley Osher
More informationMarlis Hochbruck 1, Michael Hönig 1 and Alexander Ostermann 2
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher REGULARIZATION OF NONLINEAR ILL-POSED PROBLEMS BY EXPONENTIAL INTEGRATORS Marlis
More informationRegular article Nonlinear multigrid methods for total variation image denoising
Comput Visual Sci 7: 199 06 (004) Digital Object Identifier (DOI) 10.1007/s00791-004-0150-3 Computing and Visualization in Science Regular article Nonlinear multigrid methods for total variation image
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationThe Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator
The Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator G. Bal (gb2030@columbia.edu) 1 W. Naetar (wolf.naetar@univie.ac.at) 2 O. Scherzer (otmar.scherzer@univie.ac.at) 2,3
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationA multigrid method for large scale inverse problems
A multigrid method for large scale inverse problems Eldad Haber Dept. of Computer Science, Dept. of Earth and Ocean Science University of British Columbia haber@cs.ubc.ca July 4, 2003 E.Haber: Multigrid
More informationAdaptive and multilevel methods for parameter identification in partial differential equations
Adaptive and multilevel methods for parameter identification in partial differential equations Barbara Kaltenbacher, University of Stuttgart joint work with Hend Ben Ameur, Université de Tunis Anke Griesbaum,
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationarxiv: v1 [math.na] 21 Aug 2014 Barbara Kaltenbacher
ENHANCED CHOICE OF THE PARAMETERS IN AN ITERATIVELY REGULARIZED NEWTON- LANDWEBER ITERATION IN BANACH SPACE arxiv:48.526v [math.na] 2 Aug 24 Barbara Kaltenbacher Alpen-Adria-Universität Klagenfurt Universitätstrasse
More informationMarlis Hochbruck 1, Michael Hönig 1 and Alexander Ostermann 2
ESAIM: M2AN 43 (29) 79 72 DOI:.5/m2an/292 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org REGULARIZATION OF NONLINEAR ILL-POSED PROBLEMS BY EXPONENTIAL INTEGRATORS Marlis Hochbruck,
More informationIterative Methods for Ill-Posed Problems
Iterative Methods for Ill-Posed Problems Based on joint work with: Serena Morigi Fiorella Sgallari Andriy Shyshkov Salt Lake City, May, 2007 Outline: Inverse and ill-posed problems Tikhonov regularization
More informationAccelerated Newton-Landweber Iterations for Regularizing Nonlinear Inverse Problems
www.oeaw.ac.at Accelerated Newton-Landweber Iterations for Regularizing Nonlinear Inverse Problems H. Egger RICAM-Report 2005-01 www.ricam.oeaw.ac.at Accelerated Newton-Landweber Iterations for Regularizing
More informationNew Algorithms for Parallel MRI
New Algorithms for Parallel MRI S. Anzengruber 1, F. Bauer 2, A. Leitão 3 and R. Ramlau 1 1 RICAM, Austrian Academy of Sciences, Altenbergerstraße 69, 4040 Linz, Austria 2 Fuzzy Logic Laboratorium Linz-Hagenberg,
More informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationCoordinate Update Algorithm Short Course Operator Splitting
Coordinate Update Algorithm Short Course Operator Splitting Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 25 Operator splitting pipeline 1. Formulate a problem as 0 A(x) + B(x) with monotone operators
More informationCoordinate Update Algorithm Short Course Proximal Operators and Algorithms
Coordinate Update Algorithm Short Course Proximal Operators and Algorithms Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 36 Why proximal? Newton s method: for C 2 -smooth, unconstrained problems allow
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationDue Giorni di Algebra Lineare Numerica (2GALN) Febbraio 2016, Como. Iterative regularization in variable exponent Lebesgue spaces
Due Giorni di Algebra Lineare Numerica (2GALN) 16 17 Febbraio 2016, Como Iterative regularization in variable exponent Lebesgue spaces Claudio Estatico 1 Joint work with: Brigida Bonino 1, Fabio Di Benedetto
More informationIntroduction to Scientific Computing
(Lecture 5: Linear system of equations / Matrix Splitting) Bojana Rosić, Thilo Moshagen Institute of Scientific Computing Motivation Let us resolve the problem scheme by using Kirchhoff s laws: the algebraic
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationA fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring
A fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring Marco Donatelli Dept. of Science and High Tecnology U. Insubria (Italy) Joint work with M. Hanke
More informationFast algorithms for the inverse medium problem. George Biros University of Pennsylvania
Fast algorithms for the inverse medium problem George Biros University of Pennsylvania Acknowledgments S. Adavani, H. Sundar, S. Rahul (grad students) C. Davatzikos, D. Shen, H. Litt (heart project) Akcelic,
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 informationVariational Image Restoration
Variational Image Restoration Yuling Jiao yljiaostatistics@znufe.edu.cn School of and Statistics and Mathematics ZNUFE Dec 30, 2014 Outline 1 1 Classical Variational Restoration Models and Algorithms 1.1
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationAM 205: lecture 19. Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods Optimality Conditions: Equality Constrained Case As another example of equality
More informationParameter Identification
Lecture Notes Parameter Identification Winter School Inverse Problems 25 Martin Burger 1 Contents 1 Introduction 3 2 Examples of Parameter Identification Problems 5 2.1 Differentiation of Data...............................
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.
More informationIMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND
IMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND JIAN-FENG CAI, BIN DONG, STANLEY OSHER, AND ZUOWEI SHEN Abstract. The variational techniques (e.g., the total variation based method []) are
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More informationVariational Methods in Image Denoising
Variational Methods in Image Denoising Jamylle Carter Postdoctoral Fellow Mathematical Sciences Research Institute (MSRI) MSRI Workshop for Women in Mathematics: Introduction to Image Analysis 22 January
More informationA parallel method for large scale time domain electromagnetic inverse problems
A parallel method for large scale time domain electromagnetic inverse problems Eldad Haber July 15, 2005 Abstract In this work we consider the solution of 3D time domain electromagnetic inverse problems
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More informationNumerical Optimization Professor Horst Cerjak, Horst Bischof, Thomas Pock Mat Vis-Gra SS09
Numerical Optimization 1 Working Horse in Computer Vision Variational Methods Shape Analysis Machine Learning Markov Random Fields Geometry Common denominator: optimization problems 2 Overview of Methods
More informationLevenberg-Marquardt method in Banach spaces with general convex regularization terms
Levenberg-Marquardt method in Banach spaces with general convex regularization terms Qinian Jin Hongqi Yang Abstract We propose a Levenberg-Marquardt method with general uniformly convex regularization
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More information(4D) Variational Models Preserving Sharp Edges. Martin Burger. Institute for Computational and Applied Mathematics
(4D) Variational Models Preserving Sharp Edges Institute for Computational and Applied Mathematics Intensity (cnt) Mathematical Imaging Workgroup @WWU 2 0.65 0.60 DNA Akrosom Flagellum Glass 0.55 0.50
More informationAM 205: lecture 19. Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality, Newton s method for optimization Today: survey of optimization methods Quasi-Newton Methods General form of quasi-newton methods: x k+1 = x k α
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationSimulation based optimization
SimBOpt p.1/52 Simulation based optimization Feb 2005 Eldad Haber haber@mathcs.emory.edu Emory University SimBOpt p.2/52 Outline Introduction A few words about discretization The unconstrained framework
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Seminario di Calcolo delle Variazioni ed Equazioni
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationConvex Optimization and l 1 -minimization
Convex Optimization and l 1 -minimization Sangwoon Yun Computational Sciences Korea Institute for Advanced Study December 11, 2009 2009 NIMS Thematic Winter School Outline I. Convex Optimization II. l
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationA multilevel, level-set method for optimizing eigenvalues in shape design problems
A multilevel, level-set method for optimizing eigenvalues in shape design problems E. Haber July 22, 2003 Abstract In this paper we consider optimal design problems that involve shape optimization. The
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 informationNon-polynomial Least-squares fitting
Applied Math 205 Last time: piecewise polynomial interpolation, least-squares fitting Today: underdetermined least squares, nonlinear least squares Homework 1 (and subsequent homeworks) have several parts
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More information9.1 Preconditioned Krylov Subspace Methods
Chapter 9 PRECONDITIONING 9.1 Preconditioned Krylov Subspace Methods 9.2 Preconditioned Conjugate Gradient 9.3 Preconditioned Generalized Minimal Residual 9.4 Relaxation Method Preconditioners 9.5 Incomplete
More informationA uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance tomography
c de Gruyter 28 J. Inv. Ill-Posed Problems 16 28), 381 396 DOI 1.1515 / JIIP.28.21 A uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance
More informationNonlinear Flows for Displacement Correction and Applications in Tomography
Nonlinear Flows for Displacement Correction and Applications in Tomography Guozhi Dong 1 and Otmar Scherzer 1,2 1 Computational Science Center, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien,
More informationMumford-Shah Level-Set Ideas for Problems in Medical Imaging Paris Dauphine 2005
1/26 Mumford-Shah Level-Set Ideas for Problems in Medical Imaging Paris Dauphine 2005 Wolfgang Ring, Marc Droske, Ronny Ramlau Institut für Mathematik, Universität Graz General Mumford-Shah like 2/26 functionals
More informationComputational methods for large distributed parameter estimation problems with possible discontinuities
Computational methods for large distributed parameter estimation problems with possible discontinuities Uri Ascher Department of Computer Science, University of British Columbia, Vancouver, BC, V6T 1Z4,
More information1 Non-negative Matrix Factorization (NMF)
2018-06-21 1 Non-negative Matrix Factorization NMF) In the last lecture, we considered low rank approximations to data matrices. We started with the optimal rank k approximation to A R m n via the SVD,
More informationTuning of Fuzzy Systems as an Ill-Posed Problem
Tuning of Fuzzy Systems as an Ill-Posed Problem Martin Burger 1, Josef Haslinger 2, and Ulrich Bodenhofer 2 1 SFB F 13 Numerical and Symbolic Scientific Computing and Industrial Mathematics Institute,
More informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Colloquium of the Department of Mathematics Saarland
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationTikhonov Replacement Functionals for Iteratively Solving Nonlinear Operator Equations
Tikhonov Replacement Functionals for Iteratively Solving Nonlinear Operator Equations Ronny Ramlau Gerd Teschke April 13, 25 Abstract We shall be concerned with the construction of Tikhonov based iteration
More informationIterative Methods for Smooth Objective Functions
Optimization Iterative Methods for Smooth Objective Functions Quadratic Objective Functions Stationary Iterative Methods (first/second order) Steepest Descent Method Landweber/Projected Landweber Methods
More informationNumerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation
Numerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation Sören Bartels University of Freiburg, Germany Joint work with G. Buttazzo (U Pisa), L. Diening (U Bielefeld),
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationA Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization
A Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization Alexandra Smirnova Rosemary A Renaut March 27, 2008 Abstract The preconditioned iteratively regularized Gauss-Newton
More informationNovel integro-differential equations in image processing and its applications
Novel integro-differential equations in image processing and its applications Prashant Athavale a and Eitan Tadmor b a Institute of Pure and Applied Mathematics, University of California, Los Angeles,
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More informationA Four-Pixel Scheme for Singular Differential Equations
A Four-Pixel Scheme for Singular Differential Equations Martin Welk 1, Joachim Weickert 1, and Gabriele Steidl 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer Science, Bldg. 7 Saarland
More informationA Total Fractional-Order Variation Model for Image Restoration with Non-homogeneous Boundary Conditions and its Numerical Solution
A Total Fractional-Order Variation Model for Image Restoration with Non-homogeneous Boundary Conditions and its Numerical Solution Jianping Zhang and Ke Chen arxiv:509.04237v [cs.cv] 6 Sep 205 Abstract
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationA framework for the adaptive finite element solution of large inverse problems. I. Basic techniques
ICES Report 04-39 A framework for the adaptive finite element solution of large inverse problems. I. Basic techniques Wolfgang Bangerth Center for Subsurface Modeling, Institute for Computational Engineering
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
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 informationElectrical impedance tomography with resistor networks
IOP PUBLISHING Inverse Problems 24 (2008) 035013 (31pp) INVERSE PROBLEMS doi:10.1088/0266-5611/24/3/035013 Electrical impedance tomography with resistor networks Liliana Borcea 1, Vladimir Druskin 2 and
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
More informationNovel tomography techniques and parameter identification problems
Novel tomography techniques and parameter identification problems Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München Colloquium of the Institute of Biomathematics
More informationOn effective methods for implicit piecewise smooth surface recovery
On effective methods for implicit piecewise smooth surface recovery U. M. Ascher E. Haber H. Huang October 15, 2004 Abstract This paper considers the problem of reconstructing a piecewise smooth model
More informationPart 4: Active-set methods for linearly constrained optimization. Nick Gould (RAL)
Part 4: Active-set methods for linearly constrained optimization Nick Gould RAL fx subject to Ax b Part C course on continuoue optimization LINEARLY CONSTRAINED MINIMIZATION fx subject to Ax { } b where
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationAn image decomposition model using the total variation and the infinity Laplacian
An image decomposition model using the total variation and the inity Laplacian Christopher Elion a and Luminita A. Vese a a Department of Mathematics, University of California Los Angeles, 405 Hilgard
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationarxiv: v1 [math.na] 30 Jan 2018
Modern Regularization Methods for Inverse Problems Martin Benning and Martin Burger December 18, 2017 arxiv:1801.09922v1 [math.na] 30 Jan 2018 Abstract Regularization methods are a key tool in the solution
More informationLaplace-distributed increments, the Laplace prior, and edge-preserving regularization
J. Inverse Ill-Posed Probl.? (????), 1 15 DOI 1.1515/JIIP.????.??? de Gruyter???? Laplace-distributed increments, the Laplace prior, and edge-preserving regularization Johnathan M. Bardsley Abstract. For
More informationA NOTE ON THE NONLINEAR LANDWEBER ITERATION. Dedicated to Heinz W. Engl on the occasion of his 60th birthday
A NOTE ON THE NONLINEAR LANDWEBER ITERATION Martin Hanke Dedicated to Heinz W. Engl on the occasion of his 60th birthday Abstract. We reconsider the Landweber iteration for nonlinear ill-posed problems.
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
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 informationOn a family of gradient type projection methods for nonlinear ill-posed problems
On a family of gradient type projection methods for nonlinear ill-posed problems A. Leitão B. F. Svaiter September 28, 2016 Abstract We propose and analyze a family of successive projection methods whose
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationLINEARIZED FUNCTIONAL MINIMIZATION FOR INVERSE MODELING
XIX International Conference on Water Resources CMWR University of Illinois at Urbana-Champaign June 7-, LINEARIZED FUNCTIONAL MINIMIZATION FOR INVERSE MODELING Brendt Wohlberg, Daniel M. Tartakovsky,
More information