Mechanizing Optimization and Statistics
|
|
- Sabrina Chandler
- 6 years ago
- Views:
Transcription
1 Mechanizing Optimization and Statistics Ashish Agarwal Yale University IBM Programming Languages Day Watson Research Center July 29, 2010 Ashish Agarwal () 1 / 34
2 Acknowledgments Optimization: Ignacio Grossmann (Carnegie Mellon) Nick Sawaya and Vikas Goel (Exxon Mobil) Indexing: Bob Harper (Carnegie Mellon) Statistics: Sooraj Bhat and Alex Gray (GeorgiaTech) Rich Vuduc (GeorgiaTech, linear algebra and autotuning) Ashish Agarwal () 2 / 34
3 Linear programs (LP) Dantzig (1982) x x 1 x x x 1 + x Ashish Agarwal () 3 / 34
4 Linear programs (LP) Dantzig (1982) x x 1 max x 1 x 2 x x x 1 + x Ashish Agarwal () 3 / 34
5 Linear programs (LP) Dantzig (1982) x x 1 max x 1 x 2 x x x 1 + x Cannot represent multiple polyhedra. Ashish Agarwal () 3 / 34
6 Declaring discrete choice with disjunction Disjunctive programs (DP), Balas (1974), Jeroslow and Lowe (1984), Raman and Grossmann (1994) x R 1 R x 1 x 1 1 x 2 1 } x 1 + x 2 5 {{ } R 1 [ ] 5 x1 8 4 x 2 7 }{{} R 2 Ashish Agarwal () 4 / 34
7 Declaring discrete choice with disjunction Disjunctive programs (DP), Balas (1974), Jeroslow and Lowe (1984), Raman and Grossmann (1994) x R 1 R x 1 x 1 1 x 2 1 x 1 + x 2 5 } {{ } R 1 [ 5 x1 8 4 x 2 7 ] } {{ } R 2 Language of DP extends LP with disjunction Ashish Agarwal () 4 / 34
8 Declaring discrete choice with disjunction Disjunctive programs (DP), Balas (1974), Jeroslow and Lowe (1984), Raman and Grossmann (1994) x R 1 R x 1 x 1 1 x 2 1 x 1 + x 2 5 } {{ } R 1 [ 5 x1 8 4 x 2 7 ] } {{ } R 2 Language of DP extends LP with disjunction Few algorithms for solving DPs directly. Ashish Agarwal () 4 / 34
9 Declaring discrete choice with integers Mixed-integer linear programs (MILP) Basic idea: multiply terms by y {0, 1} 0 y 1 x 3.0y + 2.0(1 y) if y = 1, then x 3.0 if y = 0, then x 2.0 Ashish Agarwal () 5 / 34
10 Declaring discrete choice with integers Mixed-integer linear programs (MILP) Basic idea: multiply terms by y {0, 1} 0 y 1 x 3.0y + 2.0(1 y) if y = 1, then x 3.0 if y = 0, then x 2.0 Language of MILP extends LP with the integer type Ashish Agarwal () 5 / 34
11 Declaring discrete choice with integers Mixed-integer linear programs (MILP) Basic idea: multiply terms by y {0, 1} 0 y 1 x 3.0y + 2.0(1 y) if y = 1, then x 3.0 if y = 0, then x 2.0 Language of MILP extends LP with the integer type Goal: Express as DP (intuitive) Convert to MILP (accepted by solvers) Ashish Agarwal () 5 / 34
12 Multiple Transformation Techniques Available Choice affects computational efficiency x 2 best big-m reformulation x 2 convex-hull reformulation x x 1 Ashish Agarwal () 6 / 34
13 System overview disjunctive constraint [A 1 x b 1 ] [A 2 x b 2 ] x2 7 6 R R x Boolean propositions big-m convex-hull indicator constr x2 x x x y = 0 A 1 x b 1 y = 1 A 2 x b 2 Pure integer inequalities Ashish Agarwal () 7 / 34
14 Previous definition Convex-hull reformulation of is [A 1 x b 1 ] [A 2 x b 2 ] A 1 x 1 b 1 y 1 A 2 x 2 b 2 y 2 y 1 + y 2 = 1 x = x 1 + x 2 Ashish Agarwal () 8 / 34
15 Previous definition Convex-hull reformulation of is [A 1 x b 1 ] [A 2 x b 2 ] A 1 x 1 b 1 y 1 A 2 x 2 b 2 y 2 y 1 + y 2 = 1 x = x 1 + x 2 Insufficient for automation Real programs not declared in canonical matrix form How are variables introduced? Disjuncts should be bounded, how is this checked? How are variable bounds tracked? Disaggregated variables should have same bounds as those they replace (except range must include zero) Ashish Agarwal () 8 / 34
16 A syntactic foundation for mathematical programs Agarwal, Bhat, Gray, Grossmann (2010) ρ ::= [r L,r U ] [r L, ) (,r U ] real r L,r U r L, ) (,r U int {true} {false} bool e ::= x r true false note e 1 or e 2 e 1 ande 2 e e 1 + e 2 e 1 e 2 e 1 e 2 c ::= T F istrue e e 1 = e 2 e 1 e 2 c 1 c 2 c 1 c 2 x :ρ c p ::= max x1 :ρ 1,...,x m:ρ m {e c} Υ ::= Υ,x :ρ Ashish Agarwal () 9 / 34
17 Convex-hull transformation Agarwal, Bhat, Gray, Grossmann (2010) Υ c cvx c { } prop Υ c j c j { Υ c j x j Υ c A c B cvx Υ j {A,B} } { y j { x j / x } } c j c j j {A,B} j {A,B} ( x A : ρ ( x B : ρ y A : 0,1 y B : 0,1 x = x A + x B) ( y A + y B = 1 ) (c A c B ) c 1,...,xj m j Compile disjuncts and context Add bounding constraints Υ ctxt Substitute disaggregated variables in each disjunct Multiply constant terms by respective y ) Ashish Agarwal () 10 / 34
18 Example: single disjunctive constraint Input var x:<10.0, 100.0> var w:<2.0, 50.0> min x + w subject_to (x <= w) disj (x >= w + 4.0) Output var x:<10.0, 100.0> var w:<2.0, 50.0> min x + w subject_to exists y1:[0, 1] exists y2:[0, 1] exists x1:<0.0, 100.0> exists x2:<0.0, 100.0> exists w1:<0.0, 50.0> exists w2:<0.0, 50.0> w = w1 + w2, x = x1 + x2, y1 + y2 = 1, 10.0 * y1 <= x1, x1 <= * y1, 2.0 * y1 <= w1, w1 <= 50.0 * y1, x1 <= w1, 10.0 * y2 <= x2, x2 <= * y2, 2.0 * y2 <= w2, w2 <= 50.0 * y2, x2 >= w * y2 Ashish Agarwal () 11 / 34
19 Example: single disjunctive constraint Input var x:<10.0, 100.0> var w:<2.0, 50.0> min x + w subject_to (x <= w) disj (x >= w + 4.0) Output generated in MPS and AMPL formats Implemented as a DSL embedded in OCaml Output var x:<10.0, 100.0> var w:<2.0, 50.0> min x + w subject_to exists y1:[0, 1] exists y2:[0, 1] exists x1:<0.0, 100.0> exists x2:<0.0, 100.0> exists w1:<0.0, 50.0> exists w2:<0.0, 50.0> w = w1 + w2, x = x1 + x2, y1 + y2 = 1, 10.0 * y1 <= x1, x1 <= * y1, 2.0 * y1 <= w1, w1 <= 50.0 * y1, x1 <= w1, 10.0 * y2 <= x2, x2 <= * y2, 2.0 * y2 <= w2, w2 <= 50.0 * y2, x2 >= w * y2 Ashish Agarwal () 11 / 34
20 Switched flow process on off m α m β hi low α M max β M min F out Ashish Agarwal () 12 / 34
21 Switched flow process, example constraint Pump α has three kinds of mode transition dynamics: i N\{n}, istrue Z α (on,off,i) ĉ α (i) = 0.0 ˆr α (i) = R e (i) istrue Z α (off,on,i) R e (i) 2.0 ĉ α (i) = 50.0 ˆr α (i) = R e (i) istrue YY α (i) ĉ α (i) = 0.0 ˆr α (i) = 0.0 Booleans and disjunction enable the natural modeling of such logical relations between constraints. Ashish Agarwal () 13 / 34
22 Switched flow process, example constraint Pump α has three kinds of mode transition dynamics: i N\{n}, istrue Z α (on,off,i) ĉ α (i) = 0.0 ˆr α (i) = R e (i) istrue Z α (off,on,i) R e (i) 2.0 ĉ α (i) = 50.0 ˆr α (i) = R e (i) istrue YY α (i) ĉ α (i) = 0.0 ˆr α (i) = which interact with each other: i N\{n}, a {α,β}, istrue YY a (i) q Q a Z a (q,q,i) Booleans and disjunction enable the natural modeling of such logical relations between constraints. Ashish Agarwal () 13 / 34
23 Switched flow process, comparison to ILOG Concert Method #vars (#binary) #constr. (#IC) time (sec) flow-concert 1061 (874) 1080 (718) flow-ic 477 (291) 1001 (438) flow-bm 477 (291) flow-ch 1194 (631) Ashish Agarwal () 14 / 34
24 Switched flow process, comparison to ILOG Concert Method #vars (#binary) #constr. (#IC) time (sec) flow-concert 1061 (874) 1080 (718) flow-ic 477 (291) 1001 (438) flow-bm 477 (291) flow-ch 1194 (631) All 3 of our methods improve on state-of-the-art. Ashish Agarwal () 14 / 34
25 Strip packing Ashish Agarwal () 15 / 34
26 Strip packing, formulation The most natural formulation uses disjunction. min no overlapping stay in bounds length s.t. length x i + L i i N [x i + L i x j ] [x j + L j x i ] [y i H i y j ] [y j H j y i ] i,j N,i < j { 0 xi L max L i i N H i y i W i N (x i,y i ) is the position of the top-left corner of rectangle i. Ashish Agarwal () 16 / 34
27 Strip packing, comparison to expert Method #vars (#binary) #constr. (#IC) time (sec) pack12-ic 289 (264) 342 (264) 1.83 pack12-bm 289 (264) pack12-ch 1345 (264) pack12-bm-expert 289 (264) pack12-ch-expert 1345 (264) pack21-ic 883 (840) 1071 (840) pack21-bm 883 (840) pack21-ch 4243 (840) pack21-bm-expert 883 (840) pack21-ch-expert 4243 (840) 5271 > Ashish Agarwal () 17 / 34
28 Strip packing, comparison to expert Method #vars (#binary) #constr. (#IC) time (sec) pack12-ic 289 (264) 342 (264) 1.83 pack12-bm 289 (264) pack12-ch 1345 (264) pack12-bm-expert 289 (264) pack12-ch-expert 1345 (264) pack21-ic 883 (840) 1071 (840) pack21-bm 883 (840) pack21-ch 4243 (840) pack21-bm-expert 883 (840) pack21-ch-expert 4243 (840) 5271 > Our mechanizations perform just as well as expert encodings. Ashish Agarwal () 17 / 34
29 Indexing is Ubiquitous We solved a problem with 150,000 equations and 25,000 variables. How were so many equations and variables declared? Ashish Agarwal () 18 / 34
30 Indexing is Ubiquitous We solved a problem with 150,000 equations and 25,000 variables. How were so many equations and variables declared? Sometimes, use matrix notation. Often, use indexing. Ashish Agarwal () 18 / 34
31 Indexing is Ubiquitous We solved a problem with 150,000 equations and 25,000 variables. How were so many equations and variables declared? Sometimes, use matrix notation. Often, use indexing. Indexed operators: Families of equations: n i=1 x i i [1... n] x i+1 = x i + y i Indexing is a meta-programming feature Index variables i distinct from mathematical variables x Ashish Agarwal () 18 / 34
32 Complex Index Sets Arise in Real Problems Job shop scheduling: j J s S j j Pre j,s t j,s t j,s Mappings from sets to set of all sets: S Dependent types: S j depends on value of j Ashish Agarwal () 19 / 34
33 Indexing Language: Syntax Agarwal (2006) Index Expressions Index Sets (Types) ε ::= i k (ε 1,...,ε m ) ε.k ε ε 1 + ε 2 ε 1 ε 2 ε 1 ε 2 case ε of {k j ε j } m j=1 σ ::= [ε L..ε U ] i 1 : σ 1 i m : σ m case ε of {k j σ j } m j=1 λi σ σ ε σ :: κ Kinds κ ::= IndexSet i : σ κ Ashish Agarwal () 20 / 34
34 Example Index Sets set JOBS = { a, b, c } set STAGES = fn i. case i of a => { s1, s2 } b => { s1, s3, s4 } c => { s3, s4 } set JOBS_STAGES = i:jobs * STAGES[i] Explicitly: {( a, s1 ), ( a, s2 ), ( b, s1 ), ( b, s3 ), ( b, s4 ), ( c, s3 ), ( c, s4 )} Ashish Agarwal () 21 / 34
35 Memory Reduction Load this program: Other software expand this to: i [1... n] x i+1 = x i + y i x 2 = x 1 + y 1 x 3 = x 2 + y 2 x 4 = x 3 + y 3 x 5 = x 4 + y 4... Ashish Agarwal () 22 / 34
36 Memory Reduction Load this program: Other software expand this to: i [1... n] x i+1 = x i + y i x 2 = x 1 + y 1 x 3 = x 2 + y 2 x 4 = x 3 + y 3 x 5 = x 4 + y 4... We retain indexing structure: Memory requirements reduced from O(n) to O(1). Ashish Agarwal () 22 / 34
37 Computational Improvements Input to our software: i:[1..10] w x i Our software s output: i:[1..10] 10.0 y i w i w i 90.0 y i [ 5.0 yi x i,d ] x i,d 75.0 y i w i x i,i y i d:[1..10] Ashish Agarwal () 23 / 34
38 Computational Improvements Input to our software: i:[1..10] w x i Our software s output: i:[1..10] 10.0 y i w i w i 90.0 y i [ 5.0 yi x i,d ] x i,d 75.0 y i w i x i,i y i d:[1..10] Reformulation time reduced from O(n) to O(1). Ashish Agarwal () 23 / 34
39 Indexing Language: Type System and Semantics Agarwal (2006) Usual Way: Syntax Type System Semantics Existence is prior to meaning. Alternative Way: Syntax Semantics Type System Meaning is prior to existence. Ashish Agarwal () 24 / 34
40 Indexing Language: Type System and Semantics Agarwal (2006) Usual Way: Syntax Type System Semantics Existence is prior to meaning. Alternative Way: Syntax Semantics Type System Meaning is prior to existence. Admits more programs. Possible only because all types are finitary. Ashish Agarwal () 24 / 34
41 What Are Random Variables? Wasserman (2004) says: A random variable is a mapping X : Ω R that assigns a real number X(ω) to each outcome ω. Ashish Agarwal () 25 / 34
42 What Are Random Variables? Wasserman (2004) says: However: A random variable is a mapping X : Ω R that assigns a real number X(ω) to each outcome ω. Treated as real: P(X 5) Not random: We write X Bernoulli(p) to mean that X is exactly distributed as f (x) = p x (1 p) 1 x for x {0,1} Ashish Agarwal () 25 / 34
43 What Are Random Variables? Not variables: Cannot substitute occurrences of X for anything. e.g. In P(X 5), certainly cannot replace X with its distribution. Dependence matters. e.g. Two random variables X and Y, both distributed as Bernoulli(0.5), each 0 or 1 with probability 0.5. What is P(X + Y = 2)? Perhaps 0.25? But not if Y = 1 X. Ashish Agarwal () 26 / 34
44 What Are Random Variables? Not variables: Cannot substitute occurrences of X for anything. e.g. In P(X 5), certainly cannot replace X with its distribution. Dependence matters. e.g. Two random variables X and Y, both distributed as Bernoulli(0.5), each 0 or 1 with probability 0.5. What is P(X + Y = 2)? Perhaps 0.25? But not if Y = 1 X. Random variables are neither random nor variable. Ashish Agarwal () 26 / 34
45 Previous Work Giry (1981), Jones and Plotkin (1989) Probability distributions are a monad. Kozen (1981) Formalized semantics. Ramsey and Pfeffer (2002) Efficient expectations, but discrete distributions only. Park, Pfenning, and Thrun (2004) Continuous distributions also, but support only sampling. Erwig and Kollmansberger (2006) Provide Haskell library, but discrete distributions only, computational efficiency not optimized. Ashish Agarwal () 27 / 34
46 Previous Work Giry (1981), Jones and Plotkin (1989) Probability distributions are a monad. Kozen (1981) Formalized semantics. Ramsey and Pfeffer (2002) Efficient expectations, but discrete distributions only. Park, Pfenning, and Thrun (2004) Continuous distributions also, but support only sampling. Erwig and Kollmansberger (2006) Provide Haskell library, but discrete distributions only, computational efficiency not optimized. Our goal: Unify these results in a single system. Ashish Agarwal () 27 / 34
47 Syntax: Probability Language Bhat, Agarwal, Gray, Vuduc (2010) T ::= Bool Int Real T 1 T 2 Prob T E ::= x true false r E 1 + E 2 E 1 E 2 (E 1,E 2 ) fst E snd E if E 1 then E 2 else E 3 E 1 = E 2 E 1 E 2 uniform return E let x E 1 in E 2 Ashish Agarwal () 28 / 34
48 Language: Type System Example typing rule: Γ E 1 : Prob T 1 Γ,x : T 1 E 2 : Prob T 2 Γ let x E 1 in E 2 : Prob T 2 Ashish Agarwal () 29 / 34
49 Language: Type System Example typing rule: Γ E 1 : Prob T 1 Γ,x : T 1 E 2 : Prob T 2 Γ let x E 1 in E 2 : Prob T 2 Pass: var U uniform in return (U 0.7) Fail: var U uniform in (U 0.7) Ashish Agarwal () 29 / 34
50 Gaussian Model Height (cm) Ashish Agarwal () 30 / 34
51 Mixture of Gaussians Model Height (cm) Ashish Agarwal () 31 / 34
52 Trying alternative statistical models Formulation: X i Normal(θ,1) ˆθ = arg max f (x θ) θ Ashish Agarwal () 32 / 34
53 Trying alternative statistical models Formulation: X i Normal(θ,1) ˆθ = arg max f (x θ) θ Solution: n ˆθ = 1 n i=1 x i Ashish Agarwal () 32 / 34
54 Trying alternative statistical models Formulation: X i Normal(θ,1) ˆθ = arg max f (x θ) θ Formulation: Z i Bernoulli(0.5) X i Normal((1 Z i )θ 0 + Z i θ 1,1) ˆθ = arg max f (x θ) θ Solution: n ˆθ = 1 n i=1 x i Ashish Agarwal () 32 / 34
55 Trying alternative statistical models Formulation: X i Normal(θ,1) ˆθ = arg max f (x θ) θ Formulation: Z i Bernoulli(0.5) X i Normal((1 Z i )θ 0 + Z i θ 1,1) ˆθ = arg max f (x θ) θ Solution: ˆθ = 1 n n i=1 x i Solution: (ˆθ 0, ˆθ 1 ) := rand(); while (...) for i = 1 to n do γ i := φ(x i ; ˆθ 1,1)/(φ(x i ; ˆθ 0,1)+φ(x i ; ˆθ 1,1)); ˆθ 0 := P n i=1 (1-γ i)*x i / P n i=1 (1-γ i); ˆθ 1 := P n i=1 γ i*x i / P n i=1 γ i; return (ˆθ 0, ˆθ 1 ); Ashish Agarwal () 32 / 34
56 Interactive algorithm assistant Bhat, Agarwal, Gray, Vuduc (2010) Features Status enter problems apply schemas undo/redo combinators can solve several textbook examples of MLE, incl. via EM autotuning + more sophisticated code generation is planned Ashish Agarwal () 33 / 34
57 Conclusions Automated bigm and convex-hull methods Beginnings of a formalization of probability and statistics Library of transformations Formalization of indexing provides: advances on previous MP languages: e.g. GAMS, AMPL, OPL fundamental improvements in time and space performance possible Ashish Agarwal () 34 / 34
58 Conclusions Automated bigm and convex-hull methods Beginnings of a formalization of probability and statistics Library of transformations Formalization of indexing provides: advances on previous MP languages: e.g. GAMS, AMPL, OPL fundamental improvements in time and space performance possible challenges remain: e.g. conversion of to indexed CNF e i:σ i :σ f :(i:σ σ ) i:σ not supported in current theory. {f (i)/i }e Ashish Agarwal () 34 / 34
Critical Reading of Optimization Methods for Logical Inference [1]
Critical Reading of Optimization Methods for Logical Inference [1] Undergraduate Research Internship Department of Management Sciences Fall 2007 Supervisor: Dr. Miguel Anjos UNIVERSITY OF WATERLOO Rajesh
More informationA note on : A Superior Representation Method for Piecewise Linear Functions by Li, Lu, Huang and Hu
A note on : A Superior Representation Method for Piecewise Linear Functions by Li, Lu, Huang and Hu Juan Pablo Vielma, Shabbir Ahmed and George Nemhauser H. Milton Stewart School of Industrial and Systems
More informationA note on : A Superior Representation Method for Piecewise Linear Functions
A note on : A Superior Representation Method for Piecewise Linear Functions Juan Pablo Vielma Business Analytics and Mathematical Sciences Department, IBM T. J. Watson Research Center, Yorktown Heights,
More informationA brief introduction to Logic. (slides from
A brief introduction to Logic (slides from http://www.decision-procedures.org/) 1 A Brief Introduction to Logic - Outline Propositional Logic :Syntax Propositional Logic :Semantics Satisfiability and validity
More informationA Monadic Probabilistic Language
A Monadic Probabilistic Language Sungwoo Park, Frank Pfenning, Sebastian Thrun School of Computer Science Carnegie Mellon University Pittsburgh, PA 15213 USA {gla,fp,thrun}@cs.cmu.edu ABSTRACT Motivated
More informationOptimization of Linear Coupled Component Automata
Optimization of Linear Coupled Component Automata Ashish Agarwal Ignacio E. Grossmann Department of Chemical Engineering, Carnegie Mellon University, Pittsburgh, PA 15213, USA November 28, 2006 Abstract
More informationProofs of randomized algorithms in Coq
Proofs of randomized algorithms in Coq Philippe Audebaud LIP - ENS Lyon 46, allée d Italie - 69437 Lyon - France Christine Paulin-Mohring INRIA Futurs, ProVal - Parc Orsay Université - F-91893 LRI, Univ
More informationA Probabilistic Language Based on Sampling Functions
A Probabilistic Language Based on Sampling Functions SUNGWOO PARK Pohang University of Science and Technology FRANK PFENNING Carnegie Mellon University and SEBASTIAN THRUN Stanford University As probabilistic
More informationOn mathematical programming with indicator constraints
On mathematical programming with indicator constraints Andrea Lodi joint work with P. Bonami & A. Tramontani (IBM), S. Wiese (Unibo) University of Bologna, Italy École Polytechnique de Montréal, Québec,
More informationImproved Big-M Reformulation for Generalized Disjunctive Programs
Improved Big-M Reformulation for Generalized Disjunctive Programs Francisco Trespalacios and Ignacio E. Grossmann Department of Chemical Engineering Carnegie Mellon University, Pittsburgh, PA 15213 Author
More informationA hierarchy of relaxations for nonlinear convex generalized disjunctive programming
A hierarchy of relaxations for nonlinear convex generalized disjunctive programming Juan P. Ruiz, Ignacio E. Grossmann 1 Carnegie Mellon University - Department of Chemical Engineering Pittsburgh, PA 15213
More informationWrite your own Theorem Prover
Write your own Theorem Prover Phil Scott 27 October 2016 Phil Scott Write your own Theorem Prover 27 October 2016 1 / 31 Introduction We ll work through a toy LCF style theorem prover for classical propositional
More informationECE473 Lecture 15: Propositional Logic
ECE473 Lecture 15: Propositional Logic Jeffrey Mark Siskind School of Electrical and Computer Engineering Spring 2018 Siskind (Purdue ECE) ECE473 Lecture 15: Propositional Logic Spring 2018 1 / 23 What
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 informationEmbedded probabilistic programming
Embedded probabilistic programming Oleg Kiselyov FNMOC oleg@pobox.com Chung-chieh Shan Rutgers University ccshan@cs.rutgers.edu 17 July 2009 Probabilistic inference Model (what) Pr(Reality) Pr(Obs j Reality)
More information02 Propositional Logic
SE 2F03 Fall 2005 02 Propositional Logic Instructor: W. M. Farmer Revised: 25 September 2005 1 What is Propositional Logic? Propositional logic is the study of the truth or falsehood of propositions or
More informationA categorical view of computational effects
Emily Riehl Johns Hopkins University A categorical view of computational effects C mp se::conference 1. Functions, composition, and categories 2. Categories for computational effects (monads) 3. Categories
More informationLecture 2: Syntax. January 24, 2018
Lecture 2: Syntax January 24, 2018 We now review the basic definitions of first-order logic in more detail. Recall that a language consists of a collection of symbols {P i }, each of which has some specified
More informationValidating QBF Invalidity in HOL4
Interactive Theorem Proving (ITP) 14 July, 2010 Quantified Boolean Formulae Quantified Boolean Formulae Motivation System Overview Related Work QBF = propositional logic + quantifiers over Boolean variables
More informationNICTA Advanced Course. Theorem Proving Principles, Techniques, Applications
NICTA Advanced Course Theorem Proving Principles, Techniques, Applications λ 1 CONTENT Intro & motivation, getting started with Isabelle Foundations & Principles Lambda Calculus Higher Order Logic, natural
More informationTyping λ-terms. Types. Typed λ-terms. Base Types. The Typing Relation. Advanced Formal Methods. Lecture 3: Simply Typed Lambda calculus
Course 2D1453, 200607 Advanced Formal Methods Lecture 3: Simply Typed Lambda calculus Mads Dam KTH/CSC Some material from B. Pierce: TAPL + some from G. Klein, NICTA Typing λterms The uptyped λcalculus
More informationTecniche di Verifica. Introduction to Propositional Logic
Tecniche di Verifica Introduction to Propositional Logic 1 Logic A formal logic is defined by its syntax and semantics. Syntax An alphabet is a set of symbols. A finite sequence of these symbols is called
More informationEDA045F: Program Analysis LECTURE 10: TYPES 1. Christoph Reichenbach
EDA045F: Program Analysis LECTURE 10: TYPES 1 Christoph Reichenbach In the last lecture... Performance Counters Challenges in Dynamic Performance Analysis Taint Analysis Binary Instrumentation 2 / 44 Types
More informationTruth-Functional Logic
Truth-Functional Logic Syntax Every atomic sentence (A, B, C, ) is a sentence and are sentences With ϕ a sentence, the negation ϕ is a sentence With ϕ and ψ sentences, the conjunction ϕ ψ is a sentence
More informationLikelihood, MLE & EM for Gaussian Mixture Clustering. Nick Duffield Texas A&M University
Likelihood, MLE & EM for Gaussian Mixture Clustering Nick Duffield Texas A&M University Probability vs. Likelihood Probability: predict unknown outcomes based on known parameters: P(x q) Likelihood: estimate
More informationPrinciples of Program Analysis: A Sampler of Approaches
Principles of Program Analysis: A Sampler of Approaches Transparencies based on Chapter 1 of the book: Flemming Nielson, Hanne Riis Nielson and Chris Hankin: Principles of Program Analysis Springer Verlag
More informationPrice: $25 (incl. T-Shirt, morning tea and lunch) Visit:
Three days of interesting talks & workshops from industry experts across Australia Explore new computing topics Network with students & employers in Brisbane Price: $25 (incl. T-Shirt, morning tea and
More information= x ( x 2 4 ) [factor out x] (2)
Connexions module: m10717 1 Part IIc: Propositional Logic 1 Propositional Equivalences What are the roots of x 3 4x? Well, in high-school algebra you learned how to deal with such numeric formulas: x 3
More informationPrinciples of Program Analysis: Control Flow Analysis
Principles of Program Analysis: Control Flow Analysis Transparencies based on Chapter 3 of the book: Flemming Nielson, Hanne Riis Nielson and Chris Hankin: Principles of Program Analysis. Springer Verlag
More informationIntroduction to Artificial Intelligence Propositional Logic & SAT Solving. UIUC CS 440 / ECE 448 Professor: Eyal Amir Spring Semester 2010
Introduction to Artificial Intelligence Propositional Logic & SAT Solving UIUC CS 440 / ECE 448 Professor: Eyal Amir Spring Semester 2010 Today Representation in Propositional Logic Semantics & Deduction
More informationOperational Semantics
Operational Semantics Semantics and applications to verification Xavier Rival École Normale Supérieure Xavier Rival Operational Semantics 1 / 50 Program of this first lecture Operational semantics Mathematical
More informationComputational Integer Programming. Lecture 2: Modeling and Formulation. Dr. Ted Ralphs
Computational Integer Programming Lecture 2: Modeling and Formulation Dr. Ted Ralphs Computational MILP Lecture 2 1 Reading for This Lecture N&W Sections I.1.1-I.1.6 Wolsey Chapter 1 CCZ Chapter 2 Computational
More informationLimitations of OCAML records
Limitations of OCAML records The record types must be declared before they are used; a label e can belong to only one record type (otherwise fun x x.e) would have several incompatible types; we cannot
More informationDiscrete Mathematics Fall 2018 Midterm Exam Prof. Callahan. Section: NetID: Multiple Choice Question (30 questions in total, 4 points each)
Discrete Mathematics Fall 2018 Midterm Exam Prof. Callahan Section: NetID: Name: Multiple Choice Question (30 questions in total, 4 points each) 1 Consider the following propositions: f: The student got
More informationPredicate Logic. Xinyu Feng 09/26/2011. University of Science and Technology of China (USTC)
University of Science and Technology of China (USTC) 09/26/2011 Overview Predicate logic over integer expressions: a language of logical assertions, for example x. x + 0 = x Why discuss predicate logic?
More informationAn Introduction to Modal Logic III
An Introduction to Modal Logic III Soundness of Normal Modal Logics Marco Cerami Palacký University in Olomouc Department of Computer Science Olomouc, Czech Republic Olomouc, October 24 th 2013 Marco Cerami
More informationInteger vs. constraint programming. IP vs. CP: Language
Discrete Math for Bioinformatics WS 0/, by A. Bockmayr/K. Reinert,. Januar 0, 0:6 00 Integer vs. constraint programming Practical Problem Solving Model building: Language Model solving: Algorithms IP vs.
More informationThe TLA + proof system
The TLA + proof system Stephan Merz Kaustuv Chaudhuri, Damien Doligez, Leslie Lamport INRIA Nancy & INRIA-MSR Joint Centre, France Amir Pnueli Memorial Symposium New York University, May 8, 2010 Stephan
More informationProjection, Inference, and Consistency
Projection, Inference, and Consistency John Hooker Carnegie Mellon University IJCAI 2016, New York City A high-level presentation. Don t worry about the details. 2 Projection as a Unifying Concept Projection
More informationA Probabilistic Language based upon Sampling Functions
A Probabilistic Language based upon Sampling Functions Sungwoo Park Frank Pfenning Sebastian Thrun Computer Science Department Computer Science Department Carnegie Mellon University Stanford University
More informationExpectation Maximization (EM) Algorithm. Each has it s own probability of seeing H on any one flip. Let. p 1 = P ( H on Coin 1 )
Expectation Maximization (EM Algorithm Motivating Example: Have two coins: Coin 1 and Coin 2 Each has it s own probability of seeing H on any one flip. Let p 1 = P ( H on Coin 1 p 2 = P ( H on Coin 2 Select
More informationLecture 3: Semantics of Propositional Logic
Lecture 3: Semantics of Propositional Logic 1 Semantics of Propositional Logic Every language has two aspects: syntax and semantics. While syntax deals with the form or structure of the language, it is
More informationBuilding Algorithmic Polytopes. David Bremner with D. Avis, H.R. Tiwary, and O. Watanabe December 16, 2014
Building Algorithmic Polytopes David Bremner with D. Avis, H.R. Tiwary, and O. Watanabe December 16, 2014 Outline Matchings in graphs Given G = (V, E), M E is called a matching {e M v e} 1 v V A matching
More informationNotes on Inductive Sets and Induction
Notes on Inductive Sets and Induction Finite Automata Theory and Formal Languages TMV027/DIT21 Ana Bove, March 15th 2018 Contents 1 Induction over the Natural Numbers 2 1.1 Mathematical (Simple) Induction........................
More informationCritical Reading of Optimization Methods for Logical Inference [1]
Critical Reading of Optimization Methods for Logical Inference [1] Undergraduate Research Internship Department of Management Sciences Winter 2008 Supervisor: Dr. Miguel Anjos UNIVERSITY OF WATERLOO Rajesh
More information19. Logic constraints, integer variables
CS/ECE/ISyE 524 Introduction to Optimization Spring 2016 17 19. Logic constraints, integer variables If-then constraints Generalized assignment problems Logic constraints Modeling a restricted set of values
More informationPropositional and Predicate Logic - VII
Propositional and Predicate Logic - VII Petr Gregor KTIML MFF UK WS 2015/2016 Petr Gregor (KTIML MFF UK) Propositional and Predicate Logic - VII WS 2015/2016 1 / 11 Theory Validity in a theory A theory
More informationCOMP219: Artificial Intelligence. Lecture 19: Logic for KR
COMP219: Artificial Intelligence Lecture 19: Logic for KR 1 Overview Last time Expert Systems and Ontologies Today Logic as a knowledge representation scheme Propositional Logic Syntax Semantics Proof
More informationLogic: Propositional Logic (Part I)
Logic: Propositional Logic (Part I) Alessandro Artale Free University of Bozen-Bolzano Faculty of Computer Science http://www.inf.unibz.it/ artale Descrete Mathematics and Logic BSc course Thanks to Prof.
More informationPropositional Logic. Testing, Quality Assurance, and Maintenance Winter Prof. Arie Gurfinkel
Propositional Logic Testing, Quality Assurance, and Maintenance Winter 2018 Prof. Arie Gurfinkel References Chpater 1 of Logic for Computer Scientists http://www.springerlink.com/content/978-0-8176-4762-9/
More informationInteger and Constraint Programming for Batch Annealing Process Planning
Integer and Constraint Programming for Batch Annealing Process Planning Willem-Jan van Hoeve and Sridhar Tayur Tepper School of Business, Carnegie Mellon University 5000 Forbes Avenue, Pittsburgh, PA 15213,
More informationA Game Semantics for a Non-Classical Logic
Can BAŞKENT INRIA, Nancy can@canbaskent.net www.canbaskent.net October 16, 2013 Outlook of the Talk Classical (but Extended) Game Theoretical Semantics for Negation Classical Game Theoretical Semantics
More informationSmall (Explicit) Extended Formulation for Knapsack Cover Inequalities from Monotone Circuits
Small (Explicit) Extended Formulation for Knapsack Cover Inequalities from Monotone Circuits Abbas Bazzi Samuel Fiorini Sangxia Huang Ola Svensson & Samuel Fiorini Tony Huynh Stefan Weltge Extended formulations
More informationStabilization in Column Generation: numerical study
1 / 26 Stabilization in Column Generation: numerical study Artur Pessoa 3 Ruslan Sadykov 1,2 Eduardo Uchoa 3 François Vanderbeck 2,1 1 INRIA Bordeaux, France 2 Univ. Bordeaux I, France 3 Universidade Federal
More informationSAT Solvers: Theory and Practice
Summer School on Verification Technology, Systems & Applications, September 17, 2008 p. 1/98 SAT Solvers: Theory and Practice Clark Barrett barrett@cs.nyu.edu New York University Summer School on Verification
More informationDecision Procedures for Satisfiability and Validity in Propositional Logic
Decision Procedures for Satisfiability and Validity in Propositional Logic Meghdad Ghari Institute for Research in Fundamental Sciences (IPM) School of Mathematics-Isfahan Branch Logic Group http://math.ipm.ac.ir/isfahan/logic-group.htm
More informationCDM Propositional Logic
CDM Propositional Logic 1 Propositional Statements Klaus Sutner Carnegie Mellon University 10-prop-logic 2017/12/15 23:21 Semantics Normal Forms A Conundrum 3 Propositional Reasoning 4 We would like our
More informationLecture Notes on SAT Solvers & DPLL
15-414: Bug Catching: Automated Program Verification Lecture Notes on SAT Solvers & DPLL Matt Fredrikson André Platzer Carnegie Mellon University Lecture 10 1 Introduction In this lecture we will switch
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 informationSAT, NP, NP-Completeness
CS 473: Algorithms, Spring 2018 SAT, NP, NP-Completeness Lecture 22 April 13, 2018 Most slides are courtesy Prof. Chekuri Ruta (UIUC) CS473 1 Spring 2018 1 / 57 Part I Reductions Continued Ruta (UIUC)
More informationComputing the MLE and the EM Algorithm
ECE 830 Fall 0 Statistical Signal Processing instructor: R. Nowak Computing the MLE and the EM Algorithm If X p(x θ), θ Θ, then the MLE is the solution to the equations logp(x θ) θ 0. Sometimes these equations
More informationComp487/587 - Boolean Formulas
Comp487/587 - Boolean Formulas 1 Logic and SAT 1.1 What is a Boolean Formula Logic is a way through which we can analyze and reason about simple or complicated events. In particular, we are interested
More informationcse541 LOGIC FOR COMPUTER SCIENCE
cse541 LOGIC FOR COMPUTER SCIENCE Professor Anita Wasilewska Spring 2015 LECTURE 2 Chapter 2 Introduction to Classical Propositional Logic PART 1: Classical Propositional Model Assumptions PART 2: Syntax
More informationA probabilistic lambda calculus - Some preliminary investigations
A probabilistic lambda calculus - Some preliminary investigations Ugo Dal Lago, Margherita Zorzi Università di Bologna, Università di Verona June, 9-11, 2010, Torino Introduction: Λ P We present some results
More informationBetween proof theory and model theory Three traditions in logic: Syntactic (formal deduction)
Overview Between proof theory and model theory Three traditions in logic: Syntactic (formal deduction) Jeremy Avigad Department of Philosophy Carnegie Mellon University avigad@cmu.edu http://andrew.cmu.edu/
More informationScalable and Accurate Verification of Data Flow Systems. Cesare Tinelli The University of Iowa
Scalable and Accurate Verification of Data Flow Systems Cesare Tinelli The University of Iowa Overview AFOSR Supported Research Collaborations NYU (project partner) Chalmers University (research collaborator)
More informationCS 4700: Foundations of Artificial Intelligence
CS 4700: Foundations of Artificial Intelligence Bart Selman selman@cs.cornell.edu Module: Knowledge, Reasoning, and Planning Part 2 Logical Agents R&N: Chapter 7 1 Illustrative example: Wumpus World (Somewhat
More informationChange in the State of the Art of (Mixed) Integer Programming
Change in the State of the Art of (Mixed) Integer Programming 1977 Vancouver Advanced Research Institute 24 papers 16 reports 1 paper computational, 4 small instances Report on computational aspects: only
More informationLOGIC PROPOSITIONAL REASONING
LOGIC PROPOSITIONAL REASONING WS 2017/2018 (342.208) Armin Biere Martina Seidl biere@jku.at martina.seidl@jku.at Institute for Formal Models and Verification Johannes Kepler Universität Linz Version 2018.1
More informationCSE 1400 Applied Discrete Mathematics Definitions
CSE 1400 Applied Discrete Mathematics Definitions Department of Computer Sciences College of Engineering Florida Tech Fall 2011 Arithmetic 1 Alphabets, Strings, Languages, & Words 2 Number Systems 3 Machine
More informationValid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints
Valid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints Nilay Noyan Andrzej Ruszczyński March 21, 2006 Abstract Stochastic dominance relations
More informationBilevel Integer Optimization: Theory and Algorithms
: Theory and Algorithms Ted Ralphs 1 Joint work with Sahar Tahernajad 1, Scott DeNegre 3, Menal Güzelsoy 2, Anahita Hassanzadeh 4 1 COR@L Lab, Department of Industrial and Systems Engineering, Lehigh University
More informationCSCE 222 Discrete Structures for Computing. Propositional Logic. Dr. Hyunyoung Lee. !!!!!! Based on slides by Andreas Klappenecker
CSCE 222 Discrete Structures for Computing Propositional Logic Dr. Hyunyoung Lee Based on slides by Andreas Klappenecker 1 Propositions A proposition is a declarative sentence that is either true or false
More informationMVE165/MMG630, Applied Optimization Lecture 6 Integer linear programming: models and applications; complexity. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming: models and applications; complexity Ann-Brith Strömberg 2011 04 01 Modelling with integer variables (Ch. 13.1) Variables Linear programming (LP) uses continuous
More informationLogic, Optimization and Data Analytics
Logic, Optimization and Data Analytics John Hooker Carnegie Mellon University United Technologies Research Center, Cork, Ireland August 2015 Thesis Logic and optimization have an underlying unity. Ideas
More informationINTRODUCTION TO PREDICATE LOGIC HUTH AND RYAN 2.1, 2.2, 2.4
INTRODUCTION TO PREDICATE LOGIC HUTH AND RYAN 2.1, 2.2, 2.4 Neil D. Jones DIKU 2005 Some slides today new, some based on logic 2004 (Nils Andersen), some based on kernebegreber (NJ 2005) PREDICATE LOGIC:
More informationMixed Integer Linear Programming Formulations for Probabilistic Constraints
Mixed Integer Linear Programming Formulations for Probabilistic Constraints J. P. Vielma a,, S. Ahmed b, G. Nemhauser b a Department of Industrial Engineering, University of Pittsburgh 1048 Benedum Hall,
More informationMidterm sample questions
Midterm sample questions CS 585, Brendan O Connor and David Belanger October 12, 2014 1 Topics on the midterm Language concepts Translation issues: word order, multiword translations Human evaluation Parts
More informationFuzzy Propositional Logic for the Knowledge Representation
Fuzzy Propositional Logic for the Knowledge Representation Alexander Savinov Institute of Mathematics Academy of Sciences Academiei 5 277028 Kishinev Moldova (CIS) Phone: (373+2) 73-81-30 EMAIL: 23LSII@MATH.MOLDOVA.SU
More informationPredicate Logic. Xinyu Feng 11/20/2013. University of Science and Technology of China (USTC)
University of Science and Technology of China (USTC) 11/20/2013 Overview Predicate logic over integer expressions: a language of logical assertions, for example x. x + 0 = x Why discuss predicate logic?
More informationOptimization Methods in Logic
Optimization Methods in Logic John Hooker Carnegie Mellon University February 2003, Revised December 2008 1 Numerical Semantics for Logic Optimization can make at least two contributions to boolean logic.
More informationOn the Complexity of the Reflected Logic of Proofs
On the Complexity of the Reflected Logic of Proofs Nikolai V. Krupski Department of Math. Logic and the Theory of Algorithms, Faculty of Mechanics and Mathematics, Moscow State University, Moscow 119899,
More informationLecture 4: Proposition, Connectives and Truth Tables
Discrete Mathematics (II) Spring 2017 Lecture 4: Proposition, Connectives and Truth Tables Lecturer: Yi Li 1 Overview In last lecture, we give a brief introduction to mathematical logic and then redefine
More informationCS1021. Why logic? Logic about inference or argument. Start from assumptions or axioms. Make deductions according to rules of reasoning.
3: Logic Why logic? Logic about inference or argument Start from assumptions or axioms Make deductions according to rules of reasoning Logic 3-1 Why logic? (continued) If I don t buy a lottery ticket on
More informationInteger Programming ISE 418. Lecture 13. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 13 Dr. Ted Ralphs ISE 418 Lecture 13 1 Reading for This Lecture Nemhauser and Wolsey Sections II.1.1-II.1.3, II.1.6 Wolsey Chapter 8 CCZ Chapters 5 and 6 Valid Inequalities
More informationGlobal Optimization of Non-convex Generalized Disjunctive Programs: A Review on Relaxations and Solution Methods
Noname manuscript No. (will be inserted by the editor) Global Optimization of Non-convex Generalized Disjunctive Programs: A Review on Relaxations and Solution Methods Juan P. Ruiz Ignacio E. Grossmann
More informationReview. Principles of Programming Languages. Equality. The Diamond Property. The Church-Rosser Theorem. Corollaries. CSE 230: Winter 2007
CSE 230: Winter 2007 Principles of Programming Languages Lecture 12: The λ-calculus Ranjit Jhala UC San Diego Review The lambda calculus is a calculus of functions: e := x λx. e e 1 e 2 Several evaluation
More informationSection Notes 9. IP: Cutting Planes. Applied Math 121. Week of April 12, 2010
Section Notes 9 IP: Cutting Planes Applied Math 121 Week of April 12, 2010 Goals for the week understand what a strong formulations is. be familiar with the cutting planes algorithm and the types of cuts
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 informationLogic Overview, I. and T T T T F F F T F F F F
Logic Overview, I DEFINITIONS A statement (proposition) is a declarative sentence that can be assigned a truth value T or F, but not both. Statements are denoted by letters p, q, r, s,... The 5 basic logical
More informationAn Introduction to SAT Solving
An Introduction to SAT Solving Applied Logic for Computer Science UWO December 3, 2017 Applied Logic for Computer Science An Introduction to SAT Solving UWO December 3, 2017 1 / 46 Plan 1 The Boolean satisfiability
More informationInteger Programming ISE 418. Lecture 12. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 12 Dr. Ted Ralphs ISE 418 Lecture 12 1 Reading for This Lecture Nemhauser and Wolsey Sections II.2.1 Wolsey Chapter 9 ISE 418 Lecture 12 2 Generating Stronger Valid
More information03 Review of First-Order Logic
CAS 734 Winter 2014 03 Review of First-Order Logic William M. Farmer Department of Computing and Software McMaster University 18 January 2014 What is First-Order Logic? First-order logic is the study of
More informationHierarchy of Relaxations for Convex. and Extension to Nonconvex Problems
Hierarchy of Relaations for Conve Nonlinear Generalized Disjunctive Programs and Etension to Nonconve Problems Ignacio Grossmann Sangbum Lee, Juan Pablo Ruiz, Nic Sawaya Center for Advanced Process Decision-maing
More informationCOMP219: Artificial Intelligence. Lecture 19: Logic for KR
COMP219: Artificial Intelligence Lecture 19: Logic for KR 1 Overview Last time Expert Systems and Ontologies Today Logic as a knowledge representation scheme Propositional Logic Syntax Semantics Proof
More informationFundamentals of Software Engineering
Fundamentals of Software Engineering First-Order Logic Ina Schaefer Institute for Software Systems Engineering TU Braunschweig, Germany Slides by Wolfgang Ahrendt, Richard Bubel, Reiner Hähnle (Chalmers
More informationIntroduction to Solving Combinatorial Problems with SAT
Introduction to Solving Combinatorial Problems with SAT Javier Larrosa December 19, 2014 Overview of the session Review of Propositional Logic The Conjunctive Normal Form (CNF) Modeling and solving combinatorial
More informationThe Curry-Howard Isomorphism
The Curry-Howard Isomorphism Software Formal Verification Maria João Frade Departmento de Informática Universidade do Minho 2008/2009 Maria João Frade (DI-UM) The Curry-Howard Isomorphism MFES 2008/09
More informationLinearity and Passivity
Linearity and Passivity David A. 1 School of Computing University of Tasmania GPO Box 252-100 Hobart 7001 Australia Abstract A simple symmetric logic is proposed which captures both the notions of Linearity
More informationCSE507. Satisfiability Modulo Theories. Computer-Aided Reasoning for Software. Emina Torlak
Computer-Aided Reasoning for Software CSE507 Satisfiability Modulo Theories courses.cs.washington.edu/courses/cse507/18sp/ Emina Torlak emina@cs.washington.edu Today Last lecture Practical applications
More information