Mixed-Integer PDE-Constrained Optimization
|
|
- Arleen Knight
- 5 years ago
- Views:
Transcription
1 Mixed-Integer PDE-Constrained Optimization IMA Workshop on Uncertainty Quantification Pelin Cay, Bart van Bloemen Waanders, Drew Kouri, Anna Thuenen and Sven Leyffer Lehigh University, Universität Magdeburg, Argonne National Laboratory, and Sandia National Laboratories February 22, 2016
2 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 2 / 42
3 Mixed-Integer PDE-Constrained Optimization (MIPDECO) PDE-constrained MIP... u = u(t, x, y, z) infinite-dimensional! t is time index; minimize u,w x, y, z are spatial dimensions F(u, w) subject to C(u, w) = 0 u U, and w Z p (integers), u(t, x, y, z): PDE states, controls, & design parameters w discrete or integral variables MIPDECO Warning w = w(t, x, y, z) Z may be infinite-dimensional integers! 3 / 42
4 Mixed-Integer PDE-Constrained Optimization (MIPDECO) PDE-constrained MIP... u = u(t, x, y, z) infinite-dimensional! t is time index; minimize u,w x, y, z are spatial dimensions F(u, w) subject to C(u, w) = 0 u U, and w Z p (integers), u(t, x, y, z): PDE states, controls, & design parameters w discrete or integral variables MIPDECO Warning w = w(t, x, y, z) Z may be infinite-dimensional integers! It s a MIP, Jim, but not as we know it! 3 / 42
5 Grand-Challenge Applications of MIPDECO Topology optimization [Sigmund and Maute, 2013] Nuclear plant design: select core types & control flow rates [Committee, 2010] Well-selection for remediation of contaminated sites [Ozdogan, 2004] Design of next-generation solar cells [Reinke et al., 2011] Design of wind-farms [Zhang et al., 2013] Scheduling for disaster recovery: oil-spills [You and Leyffer, 2010] & wildfires [Donovan and Rideout, 2003] Design & control of gas networks, [De Wolf and Smeers, 2000, Martin et al., 2006, Zavala, 2014]... also as optimization under uncertainty 4 / 42
6 Uncertainty Quantification and MIPDECO Design of experiments, e.g. discrete sensor placement Akaike s Information Criterion: parameter & structure est. AIC: maximize log-likelihood & minimize nonzeros u minimize w {0,1} N l e k (u) T R 1 e k (u)+ k=1 i=1 w i s.t. Mw i u i Mw i where R is known co-variance 5 / 42
7 Source Inversion as MIP with PDE Constraints Simple Example: Locate number of sources to match observation ū minimize J = 1 (u ū) 2 dω least-squares fit u,w 2 Ω subject to u = w kl f kl in Ω Poisson equation k,l w kl S and w kl {0, 1} source budget k,l with Dirichlet boundary conditions u = 0 on Ω. E.g. Gaussian source term, σ > 0, centered at (x k, y l ) ( (xk, y l ) (x, y) 2 ) f kl (x, y) := exp, Motivated by porous-media flow application to determine number of boreholes, [Ozdogan, 2004, Fipki and Celi, 2008] σ 2 6 / 42
8 Source Inversion as MIP with PDE Constraints Consider 2D example with Ω = [0, 1] 2 and discretize PDE: 5-point finite-difference stencil; uniform mesh h = 1/N Denote u i,j u(ih, jh) approximation at grid points minimize u,w J h = h2 2 N (u i,j ū i,j ) 2 i,j=0 subject to 4u i,j u i,j 1 u i,j+1 u i 1,j u i+1,j h 2 = u 0,j = u N,j = u i,0 = u i,n = 0 N w kl S and w kl {0, 1} k,l=1 finite-dimensional (convex) MIQP N w kl f kl (ih, jh) k,l=1 7 / 42
9 Source Inversion as MIP with PDE Constraints Potential source locations (blue dots) on mesh Create target ū using red square sources 8 / 42
10 Source Inversion as MIP with PDE Constraints 1 Contours of TARGET ubar(x,y) 0.02 Relative error in States ubar(x,y) - u(x,y) & Source Location Target (3 sources), reconstructed sources, & error on mesh 9 / 42
11 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 10 / 42
12 Mixed-Integer PDE-Constrained Optimization (MIPDECO) minimize u,w F(u, w) subject to C(u, w) = 0 u U, and w Z p (integers), u(t, x, y, z): PDE states, controls, & design parameters w discrete or integral variables Towards a problem characterization Type of PDE: different classes of PDEs e.g. elliptic, parabolic, hyperbolic, nonlinear,... Class of Integers: binary, general integers, etc Type of Objective: functional form of objective Type of Constraints: characterize c/s other than PDE Discretization: discretization method & CUTEr classification 11 / 42
13 Mesh-Independent & Mesh-Dependent Integers Definition (Mesh-Independent & Mesh-Dependent Integers) 1 The integer variables are mesh-independent, iff number of integer variables is independent of the mesh. 2 The integer variables are mesh-dependent, iff the number of integer variables depends on the mesh. Mesh-Independent Mesh-Dependent Manageable tree Theory possible Exploding tree size Theory??? 12 / 42
14 Theoretical Challenges of MIPDECO Functional Analysis (mesh-dependent integers) Denis Ridzal: What function space is w(x, y) {0, 1}? Consistently approximate w(x, y) {0, 1} as h 0? Conjecture: {w(x, y) {0, 1}} L 2 (Ω)... e.g. binary support of Cantor set not integrable Likely need additional regularity assumptions Coupling between Discretization & Integers Discretization scheme (e.g. upwinding for wave equation) depends on direction of flow (integers). Application: gas network models with flow reversals... open postdoc position at Argonne! 13 / 42
15 Computational Challenges of MIPDECO Approaches for humongous branch-and-bound trees... e.g. 3D topology optimization with 10 9 binary variables Warm-starts for PDE-constrained optimization (nodes) Guarantees for nonconvex (nonlinear) PDE constraints... factorable programming approach hopeless for 10 9 vars! + * ^ x1 log x 2 3 x2... f (x 1, x 2 ) = x 1 log(x 2 ) + x / 42
16 MIPDECO: Two Cultures Collide Observation PDE-optimization & MIP developed separately different assumptions, methodologies, and computational kernels! PDE-Optimization Obtain good solutions efficiently Nonlinear optimization: Newton s method Iterative Krylov solvers Run on bleeding-edge HPC Mixed-Integer Programming Deliver certificate of optimality Combinatorial optimization: branch-and-cut Factors & rank-one updates Limited HPC developments Potential for Disaster, or Opportunity for Innovation! 15 / 42
17 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 16 / 42
18 Source Inversion as MIP with PDE Constraints Find number and location of sources to match observation ū minimize J = 1 (u(w) ū) 2 dω least-squares fit u,w 2 Ω subject to u = w kl f kl in Ω Poisson equation k,l w kl S and w kl {0, 1} source budget k,l MIP with convex quadratic objective on Ω = [0, 1] 2 5-point finite-difference stencil; uniform mesh h = 1/N Denote u i,j u(ih, jh) approximation at grid points 17 / 42
19 Cool MIPDECO Trick: Eliminating the PDE Discretized PDE constraint (Poisson equation) 4u i,j u i,j 1 u i,j+1 u i 1,j u i+1,j h 2 = k,l w kl f kl (ih, jh), i, j Au = w kl f kl, where w kl {0, 1} only appear on RHS! Elimination of PDE and states u(x, y, z) Au = w kl f kl u = A 1 w kl f kl = k,l k,l k,l w kl A 1 f kl Solve n 2 2 n PDEs: u (kl) := A 1 f kl Eliminate u = k,l w klu (kl) from Source Inversion 18 / 42
20 Cool MIPDECO Trick: Eliminating the PDE Eliminating u = k,l w klu (kl) in MINLP gives: minimize w subject to N J h = h2 w kl u (kl) ij ū i,j 2 i,j=0 k,l N w kl S and w kl {0, 1} k,l=1 2 Eliminates the states u (N 2 variables) Eliminates the PDE constraint (N 2 constraints)... generalizes to other PDEs (with integer controls on RHS) Simplified model is quadratic knapsack problem 19 / 42
21 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 20 / 42
22 Numerical Experience with Source Inversion Find number and location of sources to match observation ū minimize J = 1 (u(w) ū) 2 dω least-squares fit u,w 2 Ω subject to u = w kl f kl in Ω Poisson equation k,l w kl S and w kl {0, 1} source budget k,l MIP with convex quadratic objective Computational Experiments: 1 Test NLP-plus-rounding heuristic versus MINLP 2 Effect of mesh-dependent vs. mesh-independent integers Mesh-independent: pick sources from 36 potential locations Mesh-dependent: all nodes are potential locations 3 Effect of state-elimination trick 21 / 42
23 1st Example Mixed-Integer PDE-Constrained Optimization Potential source locations (blue dots) on mesh Create target ū using red square sources 22 / 42
24 Approach 1: NLP-Solve, Knapsack Rounding, and MIP Knapsack Rounding 1 Solve continuous relaxation using NLP solver 2 Round largest S locations, w i, to one & set all others to zero 23 / 42
25 Approach 1: NLP-Solve, Knapsack Rounding, and MIP Knapsack Rounding 1 Solve continuous relaxation using NLP solver 2 Round largest S locations, w i, to one & set all others to zero Error in States ubar(x,y) - u(x,y) & Source Location Error in States ubar(x,y) - u(x,y) & Source Location Knapsack-rounded NLP (left) and MINLP (right) MINLP solution better: NLP-err = > = MIP-err 23 / 42
26 Mesh-Independent Source Inversion: MINLP Solvers Number of Nodes and CPU time for Increasing Mesh Sizes Nodes CPU Time Mesh-Size BonminOA MINLP Minotaur Mesh-Size BonminOA MINLP Minotaur Number of Nodes independent of mesh size! MINLP & Minotaur: filtersqp runs out of memory for N 32 BonminOA takes roughly 100 iterations... quadratic objective 24 / 42
27 Mesh-Dependent (all) Source Inversion: MINLP Solvers Number of Nodes and CPU time for Increasing Mesh Sizes Nodes CPU Time Mesh-Size BonminOA BonminBB MINLP Minotaur Mesh-Size BonminOA BonminBB MINLP Minotaur Number of nodes explodes with mesh size! OA <BREAK> after 130,000 seconds 25 / 42
28 Elimination of States & PDEs: Source Inversion CPU Time for Increasing Mesh Sizes: Simplified vs. Original Model CPU Time Mesh-Size MINLP MINLP-Simple Eliminating PDEs is two orders of magnitude faster! 26 / 42
29 Elimination of States & PDEs: Source Inversion CPU Time for Increasing Mesh Sizes: Simplified vs. Original Model Presolve Time Simplified Model Total Simplified Full PDE Model using NLP solve for PDE (inefficient) Presolve is cheap... simplified model solves much faster! 27 / 42
30 First Conclusions: Source Inversion Numerical Results Solve mesh-independent problems with coarse discretization Mesh-dependent instances cannot be solved Outer Approximation (Bon-OA) inefficient for these instances Trick # 1: elimination of states and PDE constraint Nonlinear solvers run into storage issues 28 / 42
31 First Conclusions: Source Inversion Numerical Results Solve mesh-independent problems with coarse discretization Mesh-dependent instances cannot be solved Outer Approximation (Bon-OA) inefficient for these instances Trick # 1: elimination of states and PDE constraint Nonlinear solvers run into storage issues... not surprising: MIPDECO trees grow like tribbles! 28 / 42
32 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 29 / 42
33 Control Regularization: Not All Norms Are Equal Poisson with Distributed Control [OPTPDE, 2014] & [Tröltzsch, 1984] minimize u u d 2 u,w L 2 (Ω) + e Γ u ds + α w 2 L x Γ subject to u + u = w + e Ω in Ω u = 0 on boundary Γ n w(t) {0, 1} L 1 or L 2 regularization term for control w(t) {0, 1}? Good Norms for MIPs MIP ers prefer polyhedral norms... promote integrality Old MIP trick: w 2 (t) = w(t) for w(t) {0, 1} L 1 -norm same as L 2 -norm on binary variables! 30 / 42
34 Not All Norms Are Equal Consider Distributed Control for increasing mesh-size CPU for L 2 Regularization Mesh Minotaur B-BB B-Hyb B-OA 8x x Time 32x32 Time Time Time Time 31 / 42
35 Not All Norms Are Equal Consider Distributed Control for increasing mesh-size CPU for L 2 Regularization Mesh Minotaur B-BB B-Hyb B-OA 8x x Time 32x32 Time Time Time Time CPU for L 1 Regularization Mesh Minotaur B-BB B-Hyb B-OA 8x x x L 1 regularization is equivalent to L 2, but faster Many fewer nodes in tree-searches solve up to / 42
36 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 32 / 42
37 Problem 2: Actuator Placement and Operation [Falk Hante] Goal: Control temperature with actuators Select sequence of control inputs (actuators) Choose continuous control (heat/cool) at locations Match prescribed temperature profile... de-mist bathroom mirror with hair-drier Potential Actuator Locations l = 1,..., L 33 / 42
38 Problem 2: Actuator Placement and Operation Find optimal sequence of actuators, w l (t), and controls, v l (t): where minimize u,v,w u(t f, ) 2 Ω + 2 u 2 T Ω v 2 T L subject to u t κ u = w l (t) {0, 1}, l=1 v l (t)f l in T Ω L w l (t) W, t T l=1 Lw l (t) v l (t) Uw l (t), l = 1,..., L, t T f l (x, y) = 1 ( (x, y) (xl, y l ) 2 ) ) exp 2πσ 2σ point-source for actuators at (x l, y l )... movies! 34 / 42
39 Outline 1 Introduction Problem Definition and Applications Source Inversion as MIP with PDE Constraints Problem Classification and Challenges 2 Early Theoretical & Numerical Results Eliminating the PDE & State Variables Numerical Experience with Source Inversion Control Regularization: Not All Norms Are Equal Heat Equation: Actuator Design 3 Rounding-Based Heuristic for Cloaking 4 Conclusions 35 / 42
40 Topology Design of Cloaking Devices/Scatterers Design of cloaking device on domain Ω Cloak subdomain Ω 0 (red dashes) by preventing (complex) wave from entering domain Design scatterer in subdomain ˆΩ...w(x, y) {0, 1} PDE: 2D Helmholtz (over C) with Robin boundary conditions Incident wave is exp(ik 0 y) for wavelength k 0 = 6π where i = 1 Romulan Warbird Scatterer 36 / 42
41 Topology Design of Cloaking Devices/Scatterers Control: w = w(x, y) in ˆΩ States: u = u(x, y) in Ω Target: u 0 = u 0 (x, y) in Ω 0 minimize J(u) = 1 u,v,w 2 u + u 0 2 2,Ω 0 subject to u k0 2(1 + qw)u = k2 0 qwu 0 in Ω u n ik 0u = 0 on Ω w {0, 1} in ˆΩ. Discretization: finite-differences with l = 3 nodes per scatter element, w(x, y). 37 / 42
42 Strip Rounding Heuristic Cannot solve on reasonable mesh/domain with any MINLP solver. Algorithm: Strip Rounding Heuristic Solve continuous relaxation & initialize i = 1 for i=1,...,n do Round a strip w(x i, y j ) for all j Resolve relaxation with w(x k, y) fixed for all k i end Round fractional w(x, y) following direction of wave 38 / 42
43 Strip Rounding Heuristic Cannot solve on reasonable mesh/domain with any MINLP solver. Algorithm: Strip Rounding Heuristic Solve continuous relaxation & initialize i = 1 for i=1,...,n do Round a strip w(x i, y j ) for all j Resolve relaxation with w(x k, y) fixed for all k i end Round fractional w(x, y) following direction of wave 38 / 42
44 Strip Rounding Heuristic Cannot solve on reasonable mesh/domain with any MINLP solver. Algorithm: Strip Rounding Heuristic Solve continuous relaxation & initialize i = 1 for i=1,...,n do Round a strip w(x i, y j ) for all j Resolve relaxation with w(x k, y) fixed for all k i end Round fractional w(x, y) following direction of wave 38 / 42
45 Results for Strip Rounding Scatterer, w(x, y) States u(x, y)... resolve PDE on finer mesh for fixed controls 39 / 42
46 ... Solution Not Physical! Coarse States Resolved States... not clear we re getting the correct physics! 40 / 42
47 Conclusions Mixed-Integer PDE-Constrained Optimization (MIPDECO) Class of challenging problems with important applications Subsurface flow: oil recovery or environmental remediation Design and operation of gas-/power-networks Classification: mesh-dependent vs. mesh-independent On-going work: Building library of test problems... formulation matters: interplay of binary and continuous Elimination of PDE and state variables u(t, x, y, z) Discretized PDEs huge MINLPs... push solvers to limit Outlook and Extensions Consider multi-level in space (network) and time Move toward truly multi-level approach similar to PDEs Interested in new UQ applications involving MIP & PDEs / 42
48 Our five-year mission To boldly go where no optimizer has gone before to explore strange new PDEs & MIPs! 42 / 42
49 Committee, T. (2010). Advanced fueld pellet materials and fuel rod design for water cooled reactors. Technical report, International Atomic Energy Agency. De Wolf, D. and Smeers, Y. (2000). The gas transmission problem solved by an extension of the simplex algorithm. Management Science, 46: Donovan, G. and Rideout, D. (2003). An integer programming model to optimize resource allocation for wildfire contrainment. Forest Science, 61(2). Fipki, S. and Celi, A. (2008). The use of multilateral well designs for improved recovery in heavy oil reservoirs. In IADV/SPE Conference and Exhibition, Orlanda, Florida. SPE. Martin, A., Möller, M., and Moritz, S. (2006). Mixed integer models for the stationary case of gas network optimization. Mathematical Programming, 105: OPTPDE (2014). OPTPDE a collection of problems in PDE-constrained optimization. Ozdogan, U. (2004). Optimization of well placement under time-dependent uncertainty. Master s thesis, Stanford University. 42 / 42
50 Reinke, C. M., la Mata Luque, T. M. D., Su, M. F., Sinclair, M. B., and El-Kady, I. (2011). Group-theory approach to tailored electromagnetic properties of metamaterials: An inverse-problem solution. Physical Review E, 83(6): Sigmund, O. and Maute, K. (2013). Topology optimization approaches: A comparative review. Structural and Multidisciplinary Optimization, 48(6): Tröltzsch, F. (1984). The generalized bang-bang-principle and the numerical solution of a parabolic boundary-control problem with constraints on the control and the state. Zeitschrift für Angewandte Mathematik und Mechanik, 64(12): You, F. and Leyffer, S. (2010). Oil spill response planning with MINLP. SIAG/OPT Views-and-News, 21(2):1 8. Zavala, V. M. (2014). Stochastic optimal control model for natural gas networks. Computers & Chemical Engineering, 64: Zhang, P., Romero, D., Beck, J., and Amon, C. (2013). Integration of AI and OR Techniques in Contraint Programming for Combinatorial Optimization Problems, chapter Solving Wind Farm Layout Optimization with Mixed Integer Programming and Constraint Programming. Springer Verlag. 42 / 42
Software for Integer and Nonlinear Optimization
Software for Integer and Nonlinear Optimization Sven Leyffer, leyffer@mcs.anl.gov Mathematics & Computer Science Division Argonne National Laboratory Roger Fletcher & Jeff Linderoth Advanced Methods and
More informationDevelopment of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks
Development of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks E. Wanufelle 1 S. Leyffer 2 A. Sartenaer 1 Ph. Toint 1 1 FUNDP, University of Namur 2
More informationA Branch-and-Refine Method for Nonconvex Mixed-Integer Optimization
A Branch-and-Refine Method for Nonconvex Mixed-Integer Optimization Sven Leyffer 2 Annick Sartenaer 1 Emilie Wanufelle 1 1 University of Namur, Belgium 2 Argonne National Laboratory, USA IMA Workshop,
More informationUtilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification
Utilizing Adjoint-Based Techniques to Improve the Accuracy and Reliability in Uncertainty Quantification Tim Wildey Sandia National Laboratories Center for Computing Research (CCR) Collaborators: E. Cyr,
More informationNumerical Optimization. Review: Unconstrained Optimization
Numerical Optimization Finding the best feasible solution Edward P. Gatzke Department of Chemical Engineering University of South Carolina Ed Gatzke (USC CHE ) Numerical Optimization ECHE 589, Spring 2011
More information4y Springer NONLINEAR INTEGER PROGRAMMING
NONLINEAR INTEGER PROGRAMMING DUAN LI Department of Systems Engineering and Engineering Management The Chinese University of Hong Kong Shatin, N. T. Hong Kong XIAOLING SUN Department of Mathematics Shanghai
More informationDevelopment of the new MINLP Solver Decogo using SCIP - Status Report
Development of the new MINLP Solver Decogo using SCIP - Status Report Pavlo Muts with Norman Breitfeld, Vitali Gintner, Ivo Nowak SCIP Workshop 2018, Aachen Table of contents 1. Introduction 2. Automatic
More informationFeasibility Pump for Mixed Integer Nonlinear Programs 1
Feasibility Pump for Mixed Integer Nonlinear Programs 1 Presenter: 1 by Pierre Bonami, Gerard Cornuejols, Andrea Lodi and Francois Margot Mixed Integer Linear or Nonlinear Programs (MILP/MINLP) Optimize
More informationSome Recent Advances in Mixed-Integer Nonlinear Programming
Some Recent Advances in Mixed-Integer Nonlinear Programming Andreas Wächter IBM T.J. Watson Research Center Yorktown Heights, New York andreasw@us.ibm.com SIAM Conference on Optimization 2008 Boston, MA
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationSolving Mixed-Integer Nonlinear Programs
Solving Mixed-Integer Nonlinear Programs (with SCIP) Ambros M. Gleixner Zuse Institute Berlin MATHEON Berlin Mathematical School 5th Porto Meeting on Mathematics for Industry, April 10 11, 2014, Porto
More informationMILP reformulation of the multi-echelon stochastic inventory system with uncertain demands
MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands Axel Nyberg Åbo Aademi University Ignacio E. Grossmann Dept. of Chemical Engineering, Carnegie Mellon University,
More informationLarge-scale optimization and decomposition methods: outline. Column Generation and Cutting Plane methods: a unified view
Large-scale optimization and decomposition methods: outline I Solution approaches for large-scaled problems: I Delayed column generation I Cutting plane methods (delayed constraint generation) 7 I Problems
More informationGreedy Kernel Techniques with Applications to Machine Learning
Greedy Kernel Techniques with Applications to Machine Learning Robert Schaback Jochen Werner Göttingen University Institut für Numerische und Angewandte Mathematik http://www.num.math.uni-goettingen.de/schaback
More informationMINLP: Theory, Algorithms, Applications: Lecture 3, Basics of Algorothms
MINLP: Theory, Algorithms, Applications: Lecture 3, Basics of Algorothms Jeff Linderoth Industrial and Systems Engineering University of Wisconsin-Madison Jonas Schweiger Friedrich-Alexander-Universität
More information23. Cutting planes and branch & bound
CS/ECE/ISyE 524 Introduction to Optimization Spring 207 8 23. Cutting planes and branch & bound ˆ Algorithms for solving MIPs ˆ Cutting plane methods ˆ Branch and bound methods Laurent Lessard (www.laurentlessard.com)
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective Vehicle routing problem 1 / 33 Cutting stock
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationMVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms
MVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms Ann-Brith Strömberg 2017 04 07 Lecture 8 Linear and integer optimization with applications
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationA Stochastic-Oriented NLP Relaxation for Integer Programming
A Stochastic-Oriented NLP Relaxation for Integer Programming John Birge University of Chicago (With Mihai Anitescu (ANL/U of C), Cosmin Petra (ANL)) Motivation: The control of energy systems, particularly
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective 1 / 24 Cutting stock problem 2 / 24 Problem description
More informationIndicator Constraints in Mixed-Integer Programming
Indicator Constraints in Mixed-Integer Programming Andrea Lodi University of Bologna, Italy - andrea.lodi@unibo.it Amaya Nogales-Gómez, Universidad de Sevilla, Spain Pietro Belotti, FICO, UK Matteo Fischetti,
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationHot-Starting NLP Solvers
Hot-Starting NLP Solvers Andreas Wächter Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu 204 Mixed Integer Programming Workshop Ohio
More informationImprovements to Benders' decomposition: systematic classification and performance comparison in a Transmission Expansion Planning problem
Improvements to Benders' decomposition: systematic classification and performance comparison in a Transmission Expansion Planning problem Sara Lumbreras & Andrés Ramos July 2013 Agenda Motivation improvement
More informationColumn Generation. i = 1,, 255;
Column Generation The idea of the column generation can be motivated by the trim-loss problem: We receive an order to cut 50 pieces of.5-meter (pipe) segments, 250 pieces of 2-meter segments, and 200 pieces
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationScalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty
Scalable Algorithms for Optimal Control of Systems Governed by PDEs Under Uncertainty Alen Alexanderian 1, Omar Ghattas 2, Noémi Petra 3, Georg Stadler 4 1 Department of Mathematics North Carolina State
More informationSolving LP and MIP Models with Piecewise Linear Objective Functions
Solving LP and MIP Models with Piecewise Linear Obective Functions Zonghao Gu Gurobi Optimization Inc. Columbus, July 23, 2014 Overview } Introduction } Piecewise linear (PWL) function Convex and convex
More informationMultiperiod Blend Scheduling Problem
ExxonMobil Multiperiod Blend Scheduling Problem Juan Pablo Ruiz Ignacio E. Grossmann Department of Chemical Engineering Center for Advanced Process Decision-making University Pittsburgh, PA 15213 1 Motivation
More informationBasic notions of Mixed Integer Non-Linear Programming
Basic notions of Mixed Integer Non-Linear Programming Claudia D Ambrosio CNRS & LIX, École Polytechnique 5th Porto Meeting on Mathematics for Industry, April 10, 2014 C. D Ambrosio (CNRS) April 10, 2014
More informationTopology Optimization for the Microstructure Design of Plasmonic Composites
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Topology Optimization for the Microstructure Design of Plasmonic Composites Masaki Otomori 1,
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationOutline. 1 Introduction. 2 Modeling and Applications. 3 Algorithms Convex. 4 Algorithms Nonconvex. 5 The Future of MINLP? 2 / 126
Outline 1 Introduction 2 3 Algorithms Convex 4 Algorithms Nonconvex 5 The Future of MINLP? 2 / 126 Larry Leading Off 3 / 126 Introduction Outline 1 Introduction Software Tools and Online Resources Basic
More informationMachine learning, ALAMO, and constrained regression
Machine learning, ALAMO, and constrained regression Nick Sahinidis Acknowledgments: Alison Cozad, David Miller, Zach Wilson MACHINE LEARNING PROBLEM Build a model of output variables as a function of input
More informationUncertainty Quantification in Discrete Fracture Network Models
Uncertainty Quantification in Discrete Fracture Network Models Claudio Canuto Dipartimento di Scienze Matematiche, Politecnico di Torino claudio.canuto@polito.it Joint work with Stefano Berrone, Sandra
More informationCS Algorithms and Complexity
CS 50 - Algorithms and Complexity Linear Programming, the Simplex Method, and Hard Problems Sean Anderson 2/15/18 Portland State University Table of contents 1. The Simplex Method 2. The Graph Problem
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationSOLVING ELLIPTIC PDES
university-logo SOLVING ELLIPTIC PDES School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 POISSON S EQUATION Equation and Boundary Conditions Solving the Model Problem 3 THE LINEAR ALGEBRA PROBLEM
More informationLin-Kernighan Heuristic. Simulated Annealing
DM63 HEURISTICS FOR COMBINATORIAL OPTIMIZATION Lecture 6 Lin-Kernighan Heuristic. Simulated Annealing Marco Chiarandini Outline 1. Competition 2. Variable Depth Search 3. Simulated Annealing DM63 Heuristics
More informationOpen-source finite element solver for domain decomposition problems
1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University
More informationPoisson Solvers. William McLean. April 21, Return to Math3301/Math5315 Common Material.
Poisson Solvers William McLean April 21, 2004 Return to Math3301/Math5315 Common Material 1 Introduction Many problems in applied mathematics lead to a partial differential equation of the form a 2 u +
More informationNONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS. Max Gunzburger
NONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS Max Gunzburger Department of Scientific Computing Florida State University North Carolina State University, March 10, 2011
More informationExtended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications
Extended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications François Vanderbeck University of Bordeaux INRIA Bordeaux-Sud-Ouest part : Defining Extended Formulations
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationMULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW
MULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW VAN EMDEN HENSON CENTER FOR APPLIED SCIENTIFIC COMPUTING LAWRENCE LIVERMORE NATIONAL LABORATORY Abstract Since their early application to elliptic
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationLINEAR MODELS FOR CLASSIFICATION. J. Elder CSE 6390/PSYC 6225 Computational Modeling of Visual Perception
LINEAR MODELS FOR CLASSIFICATION Classification: Problem Statement 2 In regression, we are modeling the relationship between a continuous input variable x and a continuous target variable t. In classification,
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationIE418 Integer Programming
IE418: Integer Programming Department of Industrial and Systems Engineering Lehigh University 2nd February 2005 Boring Stuff Extra Linux Class: 8AM 11AM, Wednesday February 9. Room??? Accounts and Passwords
More informationMixed Integer Non Linear Programming
Mixed Integer Non Linear Programming Claudia D Ambrosio CNRS Research Scientist CNRS & LIX, École Polytechnique MPRO PMA 2016-2017 Outline What is a MINLP? Dealing with nonconvexities Global Optimization
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationFirst-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012
First-order overdetermined systems for elliptic problems John Strain Mathematics Department UC Berkeley July 2012 1 OVERVIEW Convert elliptic problems to first-order overdetermined form Control error via
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationMidterm exam CS 189/289, Fall 2015
Midterm exam CS 189/289, Fall 2015 You have 80 minutes for the exam. Total 100 points: 1. True/False: 36 points (18 questions, 2 points each). 2. Multiple-choice questions: 24 points (8 questions, 3 points
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationColumn Generation for Extended Formulations
1 / 28 Column Generation for Extended Formulations Ruslan Sadykov 1 François Vanderbeck 2,1 1 INRIA Bordeaux Sud-Ouest, France 2 University Bordeaux I, France ISMP 2012 Berlin, August 23 2 / 28 Contents
More informationCode: 101MAT4 101MT4B. Today s topics Finite-difference method in 2D. Poisson equation Wave equation
Code: MAT MTB Today s topics Finite-difference method in D Poisson equation Wave equation Finite-difference method for elliptic PDEs in D Recall that u is the short version of u x + u y Dirichlet BVP:
More informationExponential families also behave nicely under conditioning. Specifically, suppose we write η = (η 1, η 2 ) R k R p k so that
1 More examples 1.1 Exponential families under conditioning Exponential families also behave nicely under conditioning. Specifically, suppose we write η = η 1, η 2 R k R p k so that dp η dm 0 = e ηt 1
More informationComputer Vision Group Prof. Daniel Cremers. 10a. Markov Chain Monte Carlo
Group Prof. Daniel Cremers 10a. Markov Chain Monte Carlo Markov Chain Monte Carlo In high-dimensional spaces, rejection sampling and importance sampling are very inefficient An alternative is Markov Chain
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationA simple FEM solver and its data parallelism
A simple FEM solver and its data parallelism Gundolf Haase Institute for Mathematics and Scientific Computing University of Graz, Austria Chile, Jan. 2015 Partial differential equation Considered Problem
More informationApplications and algorithms for mixed integer nonlinear programming
Applications and algorithms for mixed integer nonlinear programming Sven Leyffer 1, Jeff Linderoth 2, James Luedtke 2, Andrew Miller 3, Todd Munson 1 1 Mathematics and Computer Science Division, Argonne
More informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
More informationarxiv: v1 [math.na] 6 Nov 2017
Efficient boundary corrected Strang splitting Lukas Einkemmer Martina Moccaldi Alexander Ostermann arxiv:1711.02193v1 [math.na] 6 Nov 2017 Version of November 6, 2017 Abstract Strang splitting is a well
More informationMotivating the Covariance Matrix
Motivating the Covariance Matrix Raúl Rojas Computer Science Department Freie Universität Berlin January 2009 Abstract This note reviews some interesting properties of the covariance matrix and its role
More informationTruss topology design with integer variables made easy
Mathematical Programming manuscript No. (will be inserted by the editor) Michal Kočvara Truss topology design with integer variables made easy Dedicated to Herbert Hörnlein on the occasion of his 65th
More informationReconnect 04 Introduction to Integer Programming
Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, Reconnect 04 Introduction to Integer Programming Cynthia Phillips, Sandia National Laboratories Integer programming
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationInteger Linear Programming Modeling
DM554/DM545 Linear and Lecture 9 Integer Linear Programming Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. Assignment Problem Knapsack Problem
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationHigh Performance Nonlinear Solvers
What is a nonlinear system? High Performance Nonlinear Solvers Michael McCourt Division Argonne National Laboratory IIT Meshfree Seminar September 19, 2011 Every nonlinear system of equations can be described
More informationMultilevel Monte Carlo for two phase flow and transport in random heterogeneous porous media
Eidgenössische Technische Hochschule Zürich Ecole polytechnique fédérale de Zurich Politecnico federale di Zurigo Swiss Federal Institute of Technology Zurich Multilevel Monte Carlo for two phase flow
More informationSource Reconstruction for 3D Bioluminescence Tomography with Sparse regularization
1/33 Source Reconstruction for 3D Bioluminescence Tomography with Sparse regularization Xiaoqun Zhang xqzhang@sjtu.edu.cn Department of Mathematics/Institute of Natural Sciences, Shanghai Jiao Tong University
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More informationDecentralized Control of Stochastic Systems
Decentralized Control of Stochastic Systems Sanjay Lall Stanford University CDC-ECC Workshop, December 11, 2005 2 S. Lall, Stanford 2005.12.11.02 Decentralized Control G 1 G 2 G 3 G 4 G 5 y 1 u 1 y 2 u
More informationLagrangean Decomposition for Mean-Variance Combinatorial Optimization
Lagrangean Decomposition for Mean-Variance Combinatorial Optimization Frank Baumann, Christoph Buchheim, and Anna Ilyina Fakultät für Mathematik, Technische Universität Dortmund, Germany {frank.baumann,christoph.buchheim,anna.ilyina}@tu-dortmund.de
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationInterior-Point versus Simplex methods for Integer Programming Branch-and-Bound
Interior-Point versus Simplex methods for Integer Programming Branch-and-Bound Samir Elhedhli elhedhli@uwaterloo.ca Department of Management Sciences, University of Waterloo, Canada Page of 4 McMaster
More informationDesign of optimal RF pulses for NMR as a discrete-valued control problem
Design of optimal RF pulses for NMR as a discrete-valued control problem Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Carla Tameling (Göttingen) and Benedikt Wirth
More informationRegression with Numerical Optimization. Logistic
CSG220 Machine Learning Fall 2008 Regression with Numerical Optimization. Logistic regression Regression with Numerical Optimization. Logistic regression based on a document by Andrew Ng October 3, 204
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationMixed Integer Programming Solvers: from Where to Where. Andrea Lodi University of Bologna, Italy
Mixed Integer Programming Solvers: from Where to Where Andrea Lodi University of Bologna, Italy andrea.lodi@unibo.it November 30, 2011 @ Explanatory Workshop on Locational Analysis, Sevilla A. Lodi, MIP
More informationSparse Covariance Selection using Semidefinite Programming
Sparse Covariance Selection using Semidefinite Programming A. d Aspremont ORFE, Princeton University Joint work with O. Banerjee, L. El Ghaoui & G. Natsoulis, U.C. Berkeley & Iconix Pharmaceuticals Support
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationMultiobjective Mixed-Integer Stackelberg Games
Solving the Multiobjective Mixed-Integer SCOTT DENEGRE TED RALPHS ISE Department COR@L Lab Lehigh University tkralphs@lehigh.edu EURO XXI, Reykjavic, Iceland July 3, 2006 Outline Solving the 1 General
More informationPriority Program 1253
Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations Klaus Deckelnick and Michael Hinze A note on the approximation of elliptic control problems with bang-bang
More informationConic optimization under combinatorial sparsity constraints
Conic optimization under combinatorial sparsity constraints Christoph Buchheim and Emiliano Traversi Abstract We present a heuristic approach for conic optimization problems containing sparsity constraints.
More informationParallel Simulation of Subsurface Fluid Flow
Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationReformulation and Sampling to Solve a Stochastic Network Interdiction Problem
Network Interdiction Stochastic Network Interdiction and to Solve a Stochastic Network Interdiction Problem Udom Janjarassuk Jeff Linderoth ISE Department COR@L Lab Lehigh University jtl3@lehigh.edu informs
More informationSchur Complement Technique for Advection-Diffusion Equation using Matching Structured Finite Volumes
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 8, Number 1, pp. 51 62 (213) http://campus.mst.edu/adsa Schur Complement Technique for Advection-Diffusion Equation using Matching Structured
More information