The worst case approach to UQ
|
|
- Garey Lang
- 6 years ago
- Views:
Transcription
1 The worst case approach to UQ Houman Owhadi The gods to-day stand friendly, that we may, Lovers of peace, lead on our days to age! But, since the affairs of men rests still uncertain, Let s reason with the worst that may befall. Julius Caesar, Act 5, Scene 1 William Shakespeare ( )
2 You want to certify that Problem and
3 You want to certify that Problem and You only know
4 Compute Worst and best case optimal bounds P[G(X) a] given available information.
5 Optimal Uncertainty Quantification H. Owhadi, C. Scovel, T. J. Sullivan, M. McKerns, and M. Ortiz. Optimal Uncertainty Quantification. SIAM Review, 55(2): , Robust Optimization A. Ben-Tal, L. El Ghaoui, and A.Nemirovski. Robust optimization. Princeton Series in Applied Mathematics. Princeton University Press, Princeton, NJ, 2009 D. Bertsimas, D. B. Brown, and C. Caramanis. Theory and applications of robust optimization. SIAM Rev., 53(3): , 2011 A. Ben-Tal and A. Nemirovski. Robust convex optimization. Math. Oper. Res., 23(4): , 1998 I. Elishakoff and M. Ohsaki. Optimization and Anti-Optimization of Structures Under Uncertainty. World Scientific, London, Global Sensitivity Analysis Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. (2008). Global Sensitivity Analysis: The Primer. John Wiley & Sons.
6 Set based design in the aerospace industry Bernstein JI (1998) Design methods in the aerospace industry: looking for evidence of set-based practices. Master s thesis. Massachusetts Institute of Technology, Cambridge, MA. Set based design/analysis David J. Singer, PhD., Captain Norbert Doerry, PhD., and Michael E. Buckley, What is Set-Based Design?, Presented at ASNE DAY 2009, National Harbor, MD., April 8-9, Also published in ASNE Naval Engineers Journal, 2009 Vol 121 No 4, pp B. Rustem and Howe M. Algorithms for Worst-Case Design and Applications to Risk Management. Princeton University Press, Princeton, 2002.
7
8 P. L. Chebyshev A. A. Markov M. G. Krein
9 Answer
10 Answer Markov s inequality
11 History of classical inequalities S. Karlin and W. J. Studden. Tchebycheff Systems: With Applications in Analysis and Statistics. Pure and Applied Mathematics, Vol. XV. Interscience Publishers John Wiley & Sons, New York-London-Sydney, Classical Markov-Krein theorem and classical works of Krein, Markov and Chebyshev M. G. Krein. The ideas of P. L. Cebysev and A. A. Markov in the theory of limiting values of integrals and their further development. In E. B. Dynkin, editor, Eleven papers on Analysis, Probability, and Topology, American Mathematical Society Translations, Series 2, Volume 12, pages American Mathematical Society, New York, Theory of majorization A. W. Marshall and I. Olkin. Inequalities: Theory of Majorization and its Applications, volume 143 of Mathematics in Science and Engineering. Academic Press Inc. [Harcourt Brace Jovanovich Publishers], New York, 1979
12 Connections between Chebyshev inequalities and optimization theory H. J. Godwin. On generalizations of Tchebychef s inequality. J. Amer. Statist. Assoc., 50(271): , K. Isii. On a method for generalizations of Tchebycheff s inequality. Ann. Inst. Statist. Math. Tokyo, 10(2):65 88, K. Isii. On sharpness of Tchebycheff-type inequalities. Ann. Inst. Statist. Math., 14(1): , 1962/1963. K. Isii. The extrema of probability determined by generalized moments. I. Bounded random variables. Ann. Inst. Statist. Math., 12(2): ; errata, 280, I. Olkin and J. W. Pratt. A multivariate Tchebycheff inequality. Ann. Math. Statist, 29(1): , A. W. Marshall and I. Olkin. Multivariate Chebyshev inequalities. Ann. Math. Statist., 31(4): , A. W. Marshall and I. Olkin. Inequalities: Theory of Majorization and its Applications, volume 143 of Mathematics in Science and Engineering. Academic Press Inc. [Harcourt Brace Jovanovich Publishers], New York, 1979
13 Connection between Chebyshev inequalities and optimization theory E. B. Dynkin. Sufficient statistics and extreme points. Ann. Prob., 6(5): , A. F. Karr. Extreme points of certain sets of probability measures, with applications. Math. Oper. Res., 8(1):74 85, H. Joe. Majorization, randomness and dependence for multivariate distributions. Ann. Probab., 15(3): , J. E. Smith. Generalized Chebychev inequalities: theory and applications in decision analysis. Oper. Res., 43(5): , D. Bertsimas and I. Popescu. Optimal inequalities in probability theory: a convex optimization approach. SIAM J. Optim., 15(3): (electronic), L. Vandenberghe, S. Boyd, and K. Comanor. Generalized Chebyshev bounds via semidefinite programming. SIAM Rev., 49(1):52 64 (electronic), 2007
14 Stochastic linear programming and Stochastic Optimization G. B. Dantzig. Linear programming under uncertainty. Management Sci., 1: , A. Madansky. Bounds on the expectation of a convex function of a multivariate random variable. The Annals of Mathematical Statistics, pages , 1959 A. Madansky. Inequalities for stochastic linear programming problems. Management science, 6(2): , J. Zjackovja. On minimax solutions of stochastic linear programming problems. Casopis Pest. Mat., 91: , C. C. Huang, W. T. Ziemba, and A. Ben-Tal. Bounds on the expectation of a convex function of a random variable: With applications to stochastic programming. Operations Research, 25(2): , P. Kall. Stochastric programming with recourse: upper bounds and moment problems: a review. Mathematical research, 45:86 103, 1988 Y. Ermoliev, A. Gaivoronski, and C. Nedeva. Stochastic optimization problems with incomplete information on distribution functions. SIAM Journal on Control and Optimization, 23(5): , 1985 J. R. Birge and R. J.-B. Wets. Designing approximation schemes for stochastic optimization problems, in particular for stochastic programs with recourse. Math. Prog. Stud., 27:54 102, 1986
15 Chance constrained/distributionally robust optimization A. A. Gaivoronski. A numerical method for solving stochastic programming problems with moment constraints on a distribution function. Annals of Operations Research, 31(1): , R. I. Bot N. Lorenz, and G. Wanka. Duality for linear chance-constrained optimization problems. J. Korean Math. Soc., 47(1):17 28, J. Goh and M. Sim. Distributionally robust optimization and its tractable approximations. Oper. Res., 58(4, part 1): , 2010 S. Zymler, D. Kuhn, and B. Rustem. Distributionally robust joint chance constraints with second-order moment information. Math. Program., 137(1-2, Ser.A): , L. Xu, B. Yu, and W. Liu. The distributionally robust optimization reformulation for stochastic complementarity problems. Abstr. Appl. Anal., pages 7, W. Wiesemann, D. Kuhn, and M. Sim. Distributionally robust convex optimization. Oper. Res., 62(6): , G. A. Hanasusanto, V. Roitch, D. Kuhn, and W. Wiesemann. A distributionally robust perspective on uncertainty quantification and chance constrained programming. Mathematical Programming, 151(1):35 62, 2015 Value at Risk Artzner, P.; Delbaen, F.; Eber, J. M.; Heath, D. (1999). Coherent Measures of Risk. Mathematical Finance 9 (3): 203. W. Chen, M. Sim, J. Sun, and C.-P. Teo. From CVaR to uncertainty set: implications in joint chance-constrained optimization. Oper. Res., 58(2): , 2010.
16 Optimal Uncertainty Quantification H. Owhadi, C. Scovel, T. J. Sullivan, M. McKerns, and M. Ortiz. Optimal Uncertainty Quantification. SIAM Review, 55(2): , T. J. Sullivan, M. McKerns, D. Meyer, F. Theil, H. Owhadi, and M. Ortiz. Optimal uncertainty quantification for legacy data observations of Lipschitz functions. ESAIM Math. Model. Numer. Anal., 47(6): , S. Han, U. Topcu, M. Tao, H. Owhadi, and R. Murray. Convex optimal uncertainty quantification: Algorithms and a case study in energy storage placement for power grids. In American Control Conference (ACC), 2013, pages IEEE, 2013 S. Han, M. Tao, U. Topcu, H. Owhadi, and R. M. Murray. Convex optimal uncertainty quantification. SIAM Journal on Optimization, 25(23): , J Chen, MD Flood, R Sowers. Measuring the Unmeasurable: An Application of Uncertainty Quantification to Financial Portfolios, OFR WP, 2015 L Ming, W Chenglin. An improved algorithm for convex optimal uncertainty quantification with polytopic canonical form. Control Conference (CCC), 2015 H. Owhadi, C. Scovel and T. Sullivan. Brittleness of Bayesian Inference under Finite Information in a Continuous World. Electronic Journal of Statistics, vol 9, pp 1-79, arxiv: H. Owhadi and Clint Scovel. Extreme points of a ball about a measure with finite support (2015). arxiv:
17 Our proof relies on Winkler (1988, Extreme points of moment sets) Follows from an extension of Choquet theory (Phelps 2001, lectures on Choquet s theorem) by Von Weizsacker & Winkler (1979, Integral representation in the set of solutions of a generalized moment problem) Kendall (1962, Simplexes & Vector lattices) G. Winkler. On the integral representation in convex noncompact sets of tight measures. Mathematische Zeitschrift, 158(1):71 77, 1978 G. Winkler. Extreme points of moment sets. Math. Oper. Res., 13(4): , H. von Weizsacker and G. Winkler. Integral representation in the set of solutions of a generalized moment problem. Math. Ann., 246(1):23 32, 1979/80. D. G. Kendall. Simplexes and vector lattices. J. London Math. Soc., 37(1): , 1962.
18 Theorem H. Owhadi, C. Scovel, T. J. Sullivan, M. McKerns, and M. Ortiz. Optimal Uncertainty Quantification. SIAM Review, 55(2): , 2013.
19 Further Reduction of optimization variables
20 Another example: Optimal concentration inequality H. Owhadi, C. Scovel, T. J. Sullivan, M. McKerns, and M. Ortiz. Optimal Uncertainty Quantification. SIAM Review, 55(2): , McDiarmid inequality s
21 Theorem Reduction of optimization variables
22 Theorem m =2 Explicit Solution m=2 C = {(1, 1)} h C (s) =a (1 s 1 )D 1 (1 s 2 )D 2
23 Theorem m =2 Explicit Solution m=2 Corollary
24 Each piece of information is a constraint on an optimization problem. Optimization concepts (binding, active) transfer to UQ concepts Non binding constraint f Binding but non active constraint Active constraint A Extremizer/ Worst case scenario μ
25 Optimal Hoeffding= Optimal McDiarmid for m=2
26 Theorem m =3 Explicit Solution m=3
27 Seismic Safety Assessment of a Truss Structure Ground Acceleration F min( Yield Strain - Axial Strain ) We want to certify that
28 Power Spectrum Mean PS Observed PS Frequency Amplitude
29 Power Spectrum Mean PS Observed PS Frequency Amplitude
30 Power Spectrum Mean PS Observed PS Frequency Amplitude
31 White noise Filtered White Noise Model N. Lama, J. Wilsona, and G. Hutchinsona. Generation of synthetic earthquake accelograms using seismological modeling: a review. Journal of Earthquake Engineering, 4(3): , Filter Amplitude Mean PS Observed PS Frequency Ground acceleration
32 Vulnerability Curves (vs earthquake magnitude)
33 Identification of the weakest elements H. Owhadi, C. Scovel, T. J. Sullivan, M. McKerns, and M. Ortiz. Optimal Uncertainty Quantification. SIAM Review, 55(2): , 2013.
34 Caltech Small Particle Hypervelocity Impact Range Plate thickness Plate Obliquity Projectile velocity G We want to certify that Perforation area Problem
35 What do we know? Plate thickness Plate Obliquity Projectile velocity Thickness, obliquity, velocity: independent random variables Mean perforation area: in between 5.5 and 7.5 mm^2 Bounds on the sensitivity of the response function w.r. to each variable
36 We only know Worst case bound Reduction calculus
37 What if we know the response function? Deterministic surrogate model for the perforation area (in mm^2)
38 Optimal bound on the probability of non perforation Application of the reduction calculus The measure of probability can be reduced to the tensorization of 2 Dirac masses on thickness, obliquity and velocity
39 The optimization variables can be reduced to the tensorization of 2 Dirac masses on thickness, obliquity and velocity Support Points at iteration 0
40 Numerical optimization Support Points at iteration 150
41 Numerical optimization Support Points at iteration 200
42 Velocity and obliquity marginals each collapse to a single Dirac mass. The plate thickness marginal collapses to have support on the extremes of its range. Iteration 1000 Probability non-perforation maximized by distribution supported on minimal, not maximal, impact obliquity. Dirac on velocity at a non extreme value.
43 Important observations Extremizers are singular They identify key players i.e. vulnerabilities of the physical system Extremizers are attractors
44 Initialization with 3 support points per marginal Support Points at iteration 0
45 Initialization with 3 support points per marginal Support Points at iteration 500
46 Initialization with 3 support points per marginal Support Points at iteration 1000
47 Initialization with 3 support points per marginal Support Points at iteration 2155
48 Initialization with 5 support points per marginal Support Points at iteration 0
49 Initialization with 5 support points per marginal Support Points at iteration 1000
50 Initialization with 5 support points per marginal Support Points at iteration 3000
51 Initialization with 5 support points per marginal Support Points at iteration 7100
52 Unknown response function G + Legacy data Objective Constraints on input variables Constraint on the mean perf. area Modified Lipschitz continuity constraints on response function
53 Legacy Data 32 data points (steel-on-aluminium shots A48 A81) from summer 2010 at Caltech s SPHIR facility: These constrain the value of G at 32 points T. J. Sullivan, M. McKerns, D. Meyer, F. Theil, H. Owhadi, and M. Ortiz. Optimal uncertainty quantification for legacy data observations of Lipschitz functions. ESAIM Math. Model. Numer. Anal., 47(6): , 2013.
54 The numerical results demonstrate agreement with the Markov bound Only 2 data points out of 32 carry information about the optimal bound
55 Legacy Data 32 data points (steel-on-aluminium shots A48 A81) from summer 2010 at Caltech s SPHIR facility: Only A54 and A67 carry information The other 30 data points carry no information about least upper bound and could have be ignored. T. J. Sullivan, M. McKerns, D. Meyer, F. Theil, H. Owhadi, and M. Ortiz. Optimal uncertainty quantification for legacy data observations of Lipschitz functions. ESAIM Math. Model. Numer. Anal., 47(6): , 2013.
56 What if we have model uncertainty? What do we want? What do we know?
57 PSAAP numerical model Plate Obliquity=0 F Plate thickness Perforation area Projectile velocity 30 mm 2 F 0mm 2 4.5km/s v 7km/s
58 59 data/experimental points Plate Obliquity=0 G Plate thickness Perforation area Projectile velocity 30 mm 2 F 0mm 2 4.5km/s v 7km/s
59 Confidence sausage around the model Perforation Area F (h, v) C y
60 Admissible set Confidence sausage What we compute
61
62
63 30 mm 2 F 0mm 2 4.5km/s v 7km/s
64 The extremizers led to the identification of a bug in an old model
65 Caltech PSAAP Center UQ analysis L. J. Lucas, H. Owhadi, and M. Ortiz. Rigorous verification, validation, uncertainty quantification and certification through concentration-of-measure inequalities. Comput. Methods Appl. Mech. Engrg., 197(51-52): , M. M. McKerns, L. Strand, T. J. Sullivan, A. Fang, and M. A. G. Aivazis. Building a framework for predictive science. In Proceedings of the 10th Python in Science Conference (SciPy 2011), P.-H. T. Kamga, B. Li, M. McKerns, L. H. Nguyen, M. Ortiz, H. Owhadi, and T. J. Sullivan. Optimal uncertainty quantification with model uncertainty and legacy data. Journal of the Mechanics and Physics of Solids, 72:1 19, 2014 A. A. Kidane, A. Lashgari, B. Li, M. McKerns, M. Ortiz, H. Owhadi, G. Ravichandran, M. Stalzer, and T. J. Sullivan. Rigorous model-based uncertainty quantification with application to terminal ballistics. Part I: Systems with controllable inputs and small scatter. Journal of the Mechanics and Physics of Solids, 60(5): , T. J. Sullivan, M. McKerns, D. Meyer, F. Theil, H. Owhadi, and M. Ortiz. Optimal uncertainty quantification for legacy data observations of Lipschitz functions. ESAIM Math. Model. Numer. Anal., 47(6): , H. Owhadi, C. Scovel, T. J. Sullivan, M. McKerns, and M. Ortiz. Optimal Uncertainty Quantification. SIAM Review, 55(2): , 2013.
66 Reduced numerical optimization problems solved using mystic: a highly-configurable optimization framework pathos: a distributed parallel graph execution framework providing a highlevel programmatic interface to heterogeneous computing Mike McKerns
67 Important observations In presence of incomplete information on the distribution of input variables the dependence of the least upper bound on the accuracy of the model is very weak We need to extract as much information as possible from the sample/experimental data on the underlying distributions How do we reason with the worst in presence of data sampled from an unknown distribution?
68 Quantity of Interest You know μ A You observe Problem: θ(d)
69 Player I Player II Max Min Chooses θ Mean squared error Confidence error
70 Game theory and statistical decision theory John Von Neumann Abraham Wald J. Von Neumann. Zur Theorie der Gesellschaftsspiele. Math. Ann., 100(1): , 1928 J. Von Neumann and O. Morgenstern. Theory of Games and Economic Behavior. Princeton University Press, Princeton, New Jersey, A. Wald. Contributions to the theory of statistical estimation and testing hypotheses. Ann. Math. Statist., 10(4): , A. Wald. Statistical decision functions which minimize the maximum risk. Ann. of Math. (2), 46: , A. Wald. An essentially complete class of admissible decision functions. Ann. Math. Statistics, 18: , A. Wald. Statistical decision functions. Ann. Math. Statistics, 20: , 1949.
71 Deterministic zero sum game Player II Player I Player I s payoff Player I & II both have a blue and a red marble At the same time, they show each other a marble How should I & II play the game?
72 Pure strategy solution II should play blue and loose 1 in the worst case I should play red and loose 2 in the worst case
73 Mixed strategy (repeated game) solution II should play red with probability 3/8 and win 1/8 on average I should play red with probability 3/8 and loose 1/8 on average
74 3-2 J. Von Neumann -2 1 Player I s payoff
75 Max Min Pure strategy solution for Player II Optimal bound on the statistical error max min θ E(μ, θ) μ A Optimal statistical estimators max E(μ, θ) μ A
76 Max Min Mixed strategy (repeated game) solution for Player I Mixed strategy (repeated game) solution for Player II
77 Theorem Can we have equality?
78 Theorem The best mixed strategy for I and II = worst prior for II The best estimator is not random if the loss function is strictly convex A. Dvoretzky, A. Wald, and J. Wolfowitz. Elimination of randomization in certain statistical decision procedures and zero-sum two-person games. Ann. Math. Statist., 22(1):1 21, 1951.
79 Complete class theorem Estimator Non cooperative Minmax loss/error Non Bayesian Over-estimate risk Risk cooperative Bayesian loss/error Bayesian Under-estimate risk Prior
80 Further generalization of Statistical decision theory L. J. Savage. The theory of statistical decision. Journal of the American Statistical Association, 46:55 67, L. Le Cam. An extension of Wald s theory of statistical decision functions. Ann. Math. Statist., 26:69 81, 1955 L. D. Brown. Minimaxity, more or less. In Statistical Decision Theory and Related Topics V, pages Springer, L. D. Brown. An essay on statistical decision theory. Journal of the American Statistical Association, 95(452): , I. Gilboa and D. Schmeidler. Maxmin expected utility with non-unique prior. Journal of Mathematical Economics, 18(2): , 1989 A. Shapiro and A. Kleywegt. Minimax analysis of stochastic problems. Optim. Methods Softw., 17(3): , M. Sniedovich. The art and science of modeling decision-making under severe uncertainty. Decis. Mak. Manuf. Serv., 1(1-2): , 2007 M. Sniedovich. A classical decision theoretic perspective on worst-case analysis. Appl. Math., 56(5): , 2011.
81 Impact in econometrics and social sciences O. Morgenstern. Abraham Wald, Econometrica: Journal of the Econometric Society, pages , G. Tintner. Abraham Wald s contributions to econometrics. Ann. Math. Statistics, 23:21 28, R. Leonard. Von Neumann, Morgenstern, and the Creation of Game Theory: From Chess to Social Science, Cambridge University Press, If we want to make decision theory practical for UQ we need to introduce computational complexity constraints H. Owhadi and C. Scovel. Towards Machine Wald. Handbook for Uncertainty Quantication, arxiv: How do we do that? Is there a natural relation between game theory, computational complexity and numerical approximations?
82 A simple approximation problem Approximate solution x of Ax = b A: Knownn n symmetric positive definite matrix b: Unknown element of R n Based on the information that Φx = y Φ: Knownm n rank m matrix (m <n) y: KnownelementofR m
83
84 A Max Min A
85
86 Deterministic zero sum game Player I s payoff How should I and II play the game?
87 Pure strategy (classical numerical analysis) solution II should play blue and loose 1 in the worst case I should play red and loose 2 in the worst case
88 Mixed strategy (repeated game) solution II should play red with probability 3/8 and win 1/8 on average I should play red with probability 3/8 and loose 1/8 on average
89 3-2 J. Von Neumann -2 1 Player I s payoff
90 Game theoretic formulation Ax = b Max Min Abraham Wald Continuous game but as in decision theory under compactness it can be approximated by a finite game
91 Best strategy: lift minimax to measures Ax = b Max Min The best strategy for I is to play at random Player II s best strategy live in the Bayesian class of estimators
92 Player II s mixed strategy Player II s bet Ax = b AX = ξ ξ N (0,Q)
93 Player II s mixed strategy Ax = b AX = ξ Theorem ξ N (0,Q) Owhadi 2015, Multi-grid with rough coefficients and Multiresolution PDE decomposition from Hierarchical Information Games, arxiv: , SIAM Review (to appear)
94 Main Question Can we turn the process of discovery of a scalable numerical method into a UQ problem and, to some degree, solve it as such in an automated fashion? Can we use a computer, not only to implement a numerical method but also to find the method itself?
95 Example: Find a method for solving (1) as fast as possible to a given accuracy (1) div(a u) =g, x Ω, u =0, x Ω, Ω R d Ω is piec. Lip. a unif. ell. a i,j L (Ω) log 10 (a)
96 Multigrid Methods Multigrid: [Fedorenko, 1961, Brandt, 1973, Hackbusch, 1978] Multiresolution/Wavelet based methods [Brewster and Beylkin, 1995, Beylkin and Coult, 1998, Averbuch et al., 1998] Linear complexity with smooth coefficients Problem Severely affected by lack of smoothness
97 Robust/Algebraic multigrid [Mandel et al., 1999,Wan-Chan-Smith, 1999, Xu and Zikatanov, 2004, Xu and Zhu, 2008], [Ruge-Stüben, 1987] [Panayot ] Stabilized Hierarchical bases, Multilevel preconditioners [Vassilevski - Wang, 1997, 1998] [Panayot - Vassilevski, 1997] [Chow - Vassilevski, 2003] [Aksoylu- Holst, 2010] Some degree of robustness but problem remains open with rough coefficients Why? Interpolation operators are unknown Don t know how to bridge scales with rough coefficients!
98 Low Rank Matrix Decomposition methods Fast Multipole Method: [Greengard and Rokhlin, 1987] Hierarchical Matrix Method: [Hackbusch et al., 2002] [Bebendorf, 2008]: N ln 2d+8 N complexity To achieve grid-size accuracy in L 2 -norm
99 Common theme between these methods Their process of discovery is based on intuition, brilliant insight, and guesswork Can we turn this process of discovery into an algorithm?
100 Answer: YES Compute fast Play adversarial Information game Identify game Compute with partial information Find optimal strategy Owhadi 2015, Multi-grid with rough coefficients and Multiresolution PDE decomposition from Hierarchical Information Games, arxiv: , SIAM Review (to appear) Resulting method: This is a theorem N ln 3d N complexity To achieve grid-size accuracy in H 1 -norm Subsequent solves: N ln d+1 N complexity
101 ( div(a u) =g in Ω, Resulting method: u =0on Ω, H 1 0 (Ω) =W (1) a W (2) a a W (k) a < ψ, χ > a := R Ω ( ψ)t a χ =0for(ψ, χ) W (i) W (j), i 6= j Theorem For v W (k) C 1 2 k kvk a k div(a v)k L 2 (Ω) C 2 2 k kvk 2 a :=< v,v> a = R Ω ( v)t a v Looks like an eigenspace decomposition
102 u = w (1) + w (2) + + w (k) + w (k) =F.E.sol.ofPDEinW (k) Can be computed independently u 0.14 w (1) w (2) w (3) = w (4) w (5) w (6) + + Multiresolution decomposition of solution space
103 u = w (1) + w (2) + + w (k) + w (k) =F.E.sol.ofPDEinW (k) Can be computed independently B (k) :Stiffness matrix of PDE in W (k) Theorem λ max(b (k) ) λ min (B (k) ) C Just relax in W (k) to find w (k) Quacks like an eigenspace decomposition
104 Application to time dependent problems [Owhadi-Zhang 2016, From gamblets to near FFT-complexity solvers for wave and parabolic PDEs with rough coefficients] μ(x) 2 t u div(a u) =g(x, t) μ(x) t u div(a u) =g(x, t) Hyperbolic and parabolic PDEs with rough coefficients can be solved in O(N ln 3d N)(nearFFT)complexity u 0.14 w (1) w (2) w (3) = w (4) w (5) w (6) + + Swims like an eigenspace decomposition
105 V: F.E.spaceofH 1 0 (Ω) ofdim.n Theorem The decomposition V = W (1) a W (2) a a W (k) Canbeperformedandstoredin O(N ln 3d N)operations Doesn t have the complexity of an eigenspace decomposition
106 ψ (1) i χ (2) i χ (3) i χ (4) i χ (5) i χ (6) i Basis functions look like and behave like wavelets: Localized and can be used to compress the operator and locally analyze the solution space
107 H 1 0 (Ω) u div(a ) Reduced operator H 1 (Ω) g Inverse Problem R m u m g m Numerical implementation requires R m computation with partial information. φ 1,...,φ m L 2 (Ω) u m =( φ Ω 1 u,..., Ω φ mu) u m R m Missing information u H 1 0 (Ω)
108 Discovery process ( div(a u) =g in Ω, u =0on Ω, Identify underlying information game Measurement functions: φ 1,...,φ m L 2 (Ω) Player I Player II Chooses g L 2 (Ω) Sees Ω uφ 1,..., Ω uφ m kgk L 2 (Ω) 1 Chooses u L 2 (Ω) Max Min u u a kfk 2 a := R Ω ( f)t a f
109 Deterministic zero sum game Player II Player I Player I s payoff Player I & II both have a blue and a red marble At the same time, they show each other a marble How should I & II play the (repeated) game?
110 Optimal strategies are mixed strategies Optimal way to play is at random q Player II 1 q Game theory Player I p 3-2 John Von Neumann 1 p -2 1 John Nash =3pq +(1 p)(1 q) 2p(1 q) 2q(1 p) =1 3q + p(8q 3) = 1 8 for q = 3 8
111 Player A Chooses g L 2 (Ω) kgk L 2 (Ω) 1 Player B Sees R Ω uφ 1,..., R Ω uφ m Chooses u L 2 (Ω) u u a Continuous game but as in decision theory under compactness it can be approximated by a finite game Abraham Wald The best strategy for A is to play at random Player B s best strategy live in the Bayesian class of estimators
112 Player II s class of mixed strategies Pretend that player I is choosing g at random g L 2 (Ω) ξ: Randomfield ( div(a u) =g in Ω, u =0on Ω, ( div(a v) =ξ in Ω, v =0on Ω, Player II s u (x) :=E bet v(x) R Ω v(y)φ i(y) dy = R Ω u(y)φ i(y) dy, i Player II s optimal strategy? min max problem over distribution of ξ
113 Computational efficiency ξ N (0, Γ) Theorem Gamblets Elementary gambles form deterministic basis functions for player B s bet u (x) = P m i=1 ψ i(x) R Ω u(y)φ i(y) dy ψ i : Elementary gambles/bets ψ i ψ i (x) :=E ξ N (0,Γ) hv(x) RΩ i v(y)φ j(y) dy = δ i,j,j {1,...,m}
114 What are these gamblets? Depend on Γ: Covariance function of ξ (Player B s decision) (φ i ) m i=1 : Measurements functions (rules of the game) Example x 1 Γ(x, y) =δ(x y) φ i (x) =δ(x x i ) Ω x i x m a = I d a i,j L (Ω) ψ i : Polyharmonic splines [Harder-Desmarais, 1972][Duchon 1976, 1977,1978] ψ i : Rough Polyharmonic splines [Owhadi-Zhang-Berlyand 2013]
115 What is Player II s best strategy? What is Player II s best choice for Γ(x, y) =E ξ(x)ξ(y)? Γ = L R Ω ξ(x)f(x) dx N (0, kfk2 a) kfk 2 a := R Ω ( f)t a f L = div(a ) Why? See algebraic generalization
116 The recovery is optimal (Galerkin projection) Theorem If Γ = L then u (x) is the F.E. solution of (1) in span{l 1 φ i i =1,...,m} ku u k a =inf ψ span{l 1 φ i :i {1,...,m}} ku ψk a L = div(a ) ( div(a u) =g, x Ω, (1) u =0, x Ω,
117 Optimal variational properties Theorem P m i=1 w iψ i minimizes kψk a over all ψ such that R Ω φ jψ = w j for j =1,...,m Variational characterization ψ i : Unique minimizer of Theorem ( Minimize kψk a Subject to ψ H0 1(Ω) and R Ω φ jψ = δ i,j, j =1,...,m
118 Selection of measurement functions Example Indicator functions of a Partition of Ω of resolution H φ i =1 τi τ i Ω τ j diam(τ i ) H Theorem ku u k a H λ min (a) kgk L 2 (Ω)
119 Elementary gamble ψ i Your best bet on the value of u given the information that R τ i u =1and R τ j u =0forj 6= i τ i ( div(a u) =g, x Ω, (1) u =0, x Ω, τ j Ω
120 Exponential decay of gamblets Ω τ i r Theorem RΩ (B(τ i,r)) c ( ψ i ) T a ψ i e r lh kψ i k 2 a ψ i ψ i 4 x-axis slice 10 x-axis slice log ψ i
121 Localization of the computation of gamblets ψ loc,r i ( Minimize : Minimizer of Subject to kψk a ψ H 1 0 (S r )and R S r φ j ψ = δ i,j τ i S r r Ω No loss of accuracy if localization H ln 1 H u,loc (x) = P m for τ j S r i=1 ψloc,r i (x) R Ω u(y)φ i(y) dy Theorem If r CH ln 1 H ku u,loc 1 k a Hkgk L λmin (a) 2 (Ω)
122 Formulation of the hierarchical game
123 Hierarchy of nested Measurement functions φ (k) i 1,...,i k φ (k) i Example with k {1,...,q} = P j c i,jφ (k+1) i,j φ (1) i 1 φ (2) i 1,j 1 φ (2) i 1,j 2 φ (2) i 1,j 3 φ (2) i 1,j 4 φ (3) i 1,j 2,k 1 φ (3) i 1,j 2,k 2 φ (3) i 1,j 2,k 3 φ (3) i 1,j 2,k 4 : Indicator functions of a hierarchical nested partition of Ω of resolution H k =2 k φ (k) i τ (1) 2 τ (2) 2,3 τ (3) 2,3,1 φ (1) 2 =1 τ (1) 2 φ (2) 2,3 =1 τ (2) 2,3 φ (3) 2,3,1 =1 τ (3) 2,3,1
124 In the discrete setting simply aggregate elements (as in algebraic multigrid) i Π 1,2 i j Π 1,2 i Π 2 I τ (1) i τ (2) j Π 1,3 i Π 2,3 j I 1 I 2 I 3 φ (1) i φ (2) i φ (3) i φ (4) i φ (5) i φ (6) i
125 Formulation of the hierarchy of games Player I Chooses g L 2 (Ω) kgk L2 (Ω) 1 ( div(a u) =g in Ω, u =0on Ω, Player II Sees { uφ (k),i I Ω i k } Must predict u and { uφ (k+1),j I Ω j k+1 }
126 Player II s best strategy ( div(a u) =g in Ω, u =0on Ω, ξ N (0, L) ( div(a v) =ξ in Ω, v =0on Ω, Player II s bets u (k) (x) :=E v(x) R Ω v(y)φ(k) i (y) dy = R Ω u(y)φ(k) i (y) dy, i I k The sequence of approximations forms a martingale under the mixed strategy emerging from the game F k = σ( R Ω vφ(k) i,i I k ) v (k) (x) :=E v(x) F k Theorem F k F k+1 v (k) (x) :=E v (k+1) (x) F k
127 Player II s best strategy ( div(a u) =g in Ω, u =0on Ω, ξ N (0, L) ( div(a v) =ξ in Ω, v =0on Ω, Player II s bets u (k) (x) :=E v(x) R Ω v(y)φ(k) i (y) dy = R Ω u(y)φ(k) i (y) dy, i I k u (1) u(2) u (3) u (4) u (5) u (6)
128 Gamblets Elementary gambles form a hierarchy of deterministic basis functions for player II s hierarchy of bets Theorem u (k) (x) = P i ψ(k) i ψ (k) i R (x) Ω u(y)φ(k) (y) dy i : Elementary gambles/bets at resolution H k =2 k h ψ (k) i (x) :=E v(x) R i Ω v(y)φ(k) j (y) dy = δ i,j,j I k ψ (1) i ψ (2) i ψ (3) i ψ (4) i ψ (5) i ψ (6) i
129 Gamblets are nested V (k) := span{ψ (k),i I i k } ψ (1) i 1 Theorem V (k) V (k+1) ψ (2) i 1,j 1 ψ (2) i 1,j 2 ψ (2) i 1,j 3 ψ (2) i 1,j 4 ψ (3) i 1,j 2,k 1 ψ (3) i 1,j 2,k 2 ψ (3) i 1,j 2,k 3 ψ (3) i 1,j 2,k 4 ψ (k) i (x) = P j I k+1 R (k) i,j ψ(k+1) j (x)
130 Interpolation/Prolongation operator R (k) i,j = E R Ω v(y)φ(k+1) j (y) dy R Ω v(y)φ(k) l (y) dy = δ i,l,l I k R (k) i,j τ (k) i R Your best bet on the value of τ (k+1) j given the information that R τ (k) i u =1and R τ l u =0forl 6= i R (k) i,j u τ (k+1) j
131 At this stage you can finish with classical multigrid But we want multiresolution decomposition
132 Elementary gamble χ (k) i R τ (k) i Your best bet on the value of u given the information that u =1, R τ (k) i u = 1 and R τ (k) j τ (k) i u =0forj 6= i τ (k) i τ (k) j Ω
133 χ (k) i = ψ (k) i ψ (k) i ψ (1) i 1 i =(i 1,...,i k 1,i k ) i =(i 1,...,i k 1,i k 1) ψ (2) i 1,j 1 ψ (2) i 1,j 2 ψ (2) i 1,j 3 ψ (2) i 1,j
134 χ (k) i = ψ (k) i ψ (k) i ψ (1) i χ (2) i χ (3) i χ (4) i χ (5) i χ (6) i
135 Multiresolution decomposition of the solution space V (k) := span{ψ (k) i,i I k } W (k) := span{χ (k),i} i W (k+1) : Orthogonal complement of V (k) in V (k+1) with respect to < ψ, χ > a := R Ω ( ψ)t a χ Theorem H 1 0 (Ω) =V (1) a W (2) a a W (k) a
136 Multiresolution decomposition of the solution Theorem u (k+1) u (k) = F.E. sol. of PDE in W (k+1) u = 0.14 u (1) 0.03 u (2) u (1) u (3) u (2) + + u (4) u (3) u (5) u (4) u (6) u (5) Subband solutions u (k+1) u (k) can be computed independently
137 A (k) i,j Uniformly bounded condition numbers := ψ (k) i, ψ (k) j a B (k) i,j := χ (k), χ (k) i j a Theorem λ max (B (k) ) λ min (B (k) ) C 4.5 log10 ( λ max(a (k) ) λ min (A (k) ) ) log 10 ( λ max(b (k) ) λ min (B (k) ) ) Just relax! In v W (k) to get u (k) u (k 1)
138 V = W (1) a W (2) a a W (k) Ranges of eigenvalues in V and W (k) (k =1,...,5) in log scale infψ V kψk 2 a kψk 2 L 2, sup ψ V kψk 2 a kψk 2 L 2 infψ W (k) kψk 2 a kψk 2 L 2, sup ψ W (k) kψk 2 a kψk 2 L 2
139 c (1) c (2) i j c (3) j 4 x x c (6) j c (4) j u = P i c(1) i c (5) -2 j c (6) ψ (1) i kψ (1) i k a + P q k=2 0 Pj c(k) j χ (k) j kχ (k) j k a Coefficients of the solution in the gamblet basis j
140 Operator Compression Gamblets behave like wavelets but they are adapted to the PDE and can compress its solution space u Compression ratio = 105 Energy norm relative error = 0.07 Gamblet compression Throw 99% of the coefficients
141 Fast gamblet transform O(N ln 3d N)complexity Nesting A (k) =(R (k,k+1) ) T A (k+1) R (k,k+1) Level(k) gamblets and stiffness matrices can be computed from level(k+1) gamblets and stiffness matrices Well conditioned linear systems Underlying linear systems have uniformly bounded condition numbers ψ (k) i = ψ (k+1) (i,1) + P j C(k+1),χ i,j χ (k+1) j C (k+1),χ =(B (k+1) ) 1 Z (k+1) Localization Z (k+1) j,i := (e (k+1) j e (k+1) j ) T A (k+1) e (k+1) (i,1) The nested computation can be localized without compromising accuracy or condition numbers
142 Theorem Localizing (ψ (k) i ) i Ik and (χ (k) i ) i to subdomains of size CH k ln 1 H k Cond. No (B (k),loc ) C CH k (ln 1 H k +ln 1 ² ) u u (1),loc P q 1 k=1 (u(k+1),loc u (k),loc ) a ² Theorem The number of operations to compute gamblets and achieve accuracy ² is O N ln 3d max( 1 ²,N1/d ) (and O N ln d (N 1/d )ln 1 ² for subsequent solves) Complexity Gamblet Transform Linear Solve O(N ln 3d N) O(N ln d+1 N)
143 ϕ i,a h,m h ψ (q) χ (q) i,b (q) i,a (q) u (q) u (q 1) Parallel operating diagram both in space and in frequency ψ (3) i ψ (q 1) i,a (q 1) χ (q 1) i,b (q 1) u (q 1) u (q 2) ψ (3) χ (3) i,b (3) u (1) i,a (3) u (3) u (2) ψ (2) i,a (2) χ (2) i,b (2) u (2) u (1) ψ (1) i,a (1) χ (3) i u (3) u (2) ψ (2) i χ (2) i 0.03 u (2) u (1) ψ (1) i ψ (1) i 0.14 u (1)
144 Harmonic Coordinates HMM Kozlov, 1979 Numerical Homogenization Babuska, Caloz, Osborn, 1994 Allaire Brizzi 2005; Owhadi, Zhang 2005 MsFEM [Hou, Wu: 1997]; [Efendiev, Hou, Wu: 1999] [Fish - Wagiman, 1993] [Gloria 2010] Nolen, Papanicolaou, Pironneau, 2008 Engquist, E, Abdulle, Runborg, Schwab, et Al Projection based method Flux norm Berlyand, Owhadi 2010; Symes 2012 Arbogast, 2011: Mixed MsFEM [Chu-Graham-Hou-2010] (limited inclusions) [Efendiev-Galvis-Wu-2010] (limited inclusions or mask) [Babuska-Lipton 2010] (local boundary eigenvectors) [Owhadi-Zhang 2011] (localized transfer property) [Malqvist-Peterseim 2012] Volume averaged interpolation Localization
145 Statistical approach to numerical approximation
146 Statistical approach to numerical approximation
147 Some high level remarks
148 What is the worst? u : Unknown element of A A Φ : A R u Φ(u) Quantity of Interest What is Φ(u )? inf u A Φ(u) Φ(u ) sup u A Φ(u) Robust Optimization worst case
149 A Max Min A Game theoretic worst case
150 Robust Optimization worst case Confidence error Failure is not an option. You want to always be right. Game theoretic worst case Interpretation depends on the choice of loss function. You want to be right with high probability. Quadratic error You want to be right on average. Well suited for numerical computation where you need to keep computing with partial information (e.g. invert a 1,000,000 by 1,000,000 matrix)
151 Complete class theorem Estimator Non cooperative Minmax loss/error Non Bayesian Over-estimate risk Risk cooperative Bayesian loss/error Bayesian Under-estimate risk Prior Can we approximate the optimal prior?
152 Numerical robustness of Bayesian inference Can we numerically approximate the prior when closed form expressions are not available for posterior values?
153
154 Robustness of Bayesian conditioning in continuous spaces Brittleness of Bayesian Inference under Finite Information in a Continuous World. H. Owhadi, C. Scovel and T. Sullivan. Electronic Journal of Statistics, vol 9, pp 1-79, arxiv: Brittleness of Bayesian inference and new Selberg formulas. H. Owhadi and C. Scovel. Communications in Mathematical Sciences (2015). arxiv: On the Brittleness of Bayesian Inference. H. Owhadi, C. Scovel and T. Sullivan. SIAM Review, 57(4), , 2015, arxiv: Qualitative Robustness in Bayesian Inference (2015). H. Owhadi and C. Scovel. arxiv: Positive Classical Bernstein Von Mises Wasserman, Lavine, Wolpert (1993) P Gustafson & L Wasserman (1995) Castillo and Rousseau (2013) Castillo and Nickl (2013) Stuart & Al (2010+).. Negative Freedman (1963, 1965) P Gustafson & L Wasserman (1995) Diaconis & Freedman 1998 Johnstone 2010 Leahu 2011 Belot 2013
155 Brittleness of Bayesian Inference under Finite Information in a Continuous World. H. Owhadi, C. Scovel and T. Sullivan. Electronic Journal of Statistics, vol 9, pp 1-79, arxiv: children are given one pound of play-doh. On average, how much mass can they put above a while, on average, keeping the seesaw balanced around m? Paul is given one pound of play-doh. What can you say about how much mass he is putting above a if all you have is the belief that he is keeping the seesaw balanced around m?
156 What is the least upper bound on If all you know is? Answer
157 Theorem
158
159
160 Reduction calculus with measures over measures M(X ) M(A) Ψ Ψ 1 Q Q Polish space M(Q) Theorem Brittleness of Bayesian Inference under Finite Information in a Continuous World. H. Owhadi, C. Scovel and T. Sullivan. Electronic Journal of Statistics, vol 9, pp 1-79, arxiv: sup Q Q sup π Ψ 1 Q he q Q E μ π Φ(μ) = sup μ Ψ 1 (q) Φ(μ) i
161 What is the worst with random data? A
162 A Frequentist/Concentration of measure worst case
163 Theorem N. Fournier and A. Guillin, On the rate of convergence in Wasserstein distance of the empirical measure, Probability Theory and Related Fields, (2014), pp A
164 Reduction calculus of the ball about the empirical distribution D. Wozabal. A framework for optimization under ambiguity. Annals of Operations Research, 193(1):21 47, P. M. Esfahani and D. Kuhn. Data-driven distributionally robust optimization using the Wasserstein metric: performance guarantees and tractable reformulations. arxiv: , Extreme points of a ball about a measure with finite support (2015). H. Owhadi and Clint Scovel. arxiv: The extreme points of the Prokhorov, Monge-Wasserstein and Kantorovich metric balls about a measure whose support has at most n points, consist of measures whose supports have at most n+2 points.
165 Question Game/Decision Theory + Information Based Complexity Turn the process of discovery of scalable numerical solvers into an algorithm Worst case calculus?
166 The truncated moment problem Ψ ³ E X μ [X], E X μ [X 2 ],...,E X μ [X k ] Study of the geometry of M k := Ψ M([0, 1]) P. L. Chebyshev A. A. Markov M. G. Krein
167 Ψ M k := Ψ M([0, 1]) ³ E X μ [X], E X μ [X 2 ],...,E X μ [X k ] Infinite dim. Finite dim. Finite dim.
168 Infinite dim. Finite dim. Ψ t 1 t 2 t j t N
169 Index Theorem Upper Lower t 1 t 2 t j t N t 1 t 2 t j t N
170 Selberg Identities S n (α, β, γ) = Q n 1 j=0 Γ(α+jγ)Γ(β+jγ)Γ(1+(j+1)γ) Γ(α+β+(n+j 1)γ)Γ(1+γ) S n (α, β, γ) := R [0,1] n Q n j=1 tα 1 j (1 t j ) β 1 (t) 2γ dt. (t) := Q j<k (t k t j ) Brittleness of Bayesian inference and new Selberg formulas. H. Owhadi and C. Scovel. Communications in Mathematical Sciences (2015). arxiv:
171 Forrester and Warnaar 2008 The importance of the Selberg integral Used to prove outstanding conjectures in Random matrix theory and cases of the Macdonald conjectures Central role in random matrix theory, Calogero- Sutherland quantum many-body systems, Knizhnik- Zamolodchikov equations, and multivariable orthogonal polynomial theory
172 Index Theorem t 1 t 2 t t j t N
173 New Reproducing Kernel Hilbert Spaces and Selberg Integral formulas related to the Markov-Krein representations of moment spaces. Ψ ³ E X μ [X], E X μ [X 2 ],...,E X μ [X k ] Σt RI 1 Qm m j=1 t2 j (1 t j) 2 4 m(t)dt = S m(5,1,2) S m (3,3,2) 2 Σt RI 1 Qm m j=1 t2 j 4 m(t)dt = m 2 S m 1(5, 3, 2) m (t) := Q j<k (t k t j ) I := [0, 1] Σφ (t) := P m j=1 φ(t j), t I m S n (α, β, γ) = Q n 1 j=0 Γ(α+jγ)Γ(β+jγ)Γ(1+(j+1)γ) Γ(α+β+(n+j 1)γ)Γ(1+γ)
174 e j (t) := X Theorem i 1 < <i j t i1 t ij Π n 0 : n-th degree polynomials which vanish on the boundary of [0, 1] M n R n : set of q =(q 1,...,q n ) R n such that there exists a probability measure μ on [0, 1] with E μ [X i ]=q i with i {1,...,n}. Bi-orthogonal systems of Selberg Integral formulas Consider the basis of Π 2m 1 0 consisting of the associated Legendre polynomials Q j,j = 2,..,2m 1 of order 2 translated to the unit interval I. For k =2,..,2m 1define a jk := (j + k + k2 )Γ(j +2)Γ(j) Γ(j + k +2)Γ(j k +1), k j 2m 1 Then for j = k mod 2, j, k =2,..,2m 1, we have Z m 1 Y hk (t)σq j (t) t 2 j (k +2)! 0 4 m 1(t)dt = Vol(M 2m 1 )(2m 1)!(m 1)! I (8k +4)(k 2)! δ jk. m 1 j 0 =1 h k (t) := 2m 1 X j=k ( 1) j+1 a jk e 2m 1 j (t, t).
175 Collaborators Clint Scovel (Caltech), Tim Sullivan (Warwick), Mike McKerns (Caltech), Michael Ortiz (Caltech), Lei Zhang (Jiaotong), Leonid Berlyand (PSU), Lan Huong Nguyen, Paul Herve Tamogoue Kamga. Research supported by Air Force Office of Scientific Research DARPA EQUiPS Program (Enabling Quantification of Uncertainty in Physical Systems) U.S. Department of Energy Office of Science, Office of Advanced Scientific Computing Research, through the Exascale Co-Design Center for Materials in Extreme Environments National Nuclear Security Administration
Universal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis.
Universal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis Houman Owhadi Joint work with Clint Scovel IPAM Apr 3, 2017 DARPA EQUiPS / AFOSR
More informationComputational Information Games A minitutorial Part II Houman Owhadi ICERM June 5, 2017
Computational Information Games A minitutorial Part II Houman Owhadi ICERM June 5, 2017 DARPA EQUiPS / AFOSR award no FA9550-16-1-0054 (Computational Information Games) Question Can we design a linear
More informationBayesian Brittleness
Bayesian Brittleness Houman Owhadi Bayesian Brittleness. H. Owhadi, C. Scovel, T. Sullivan. arxiv:1304.6772 Brittleness of Bayesian inference and new Selberg formulas. H. Owhadi, C. Scovel. arxiv:1304.7046
More informationComputational Information Games Houman Owhadi
Computatonal Informaton Games Houman Owhad TexAMP 215 Man Queston Can we turn the process of dscovery of a scalable numercal method nto a UQ problem and, to some degree, solve t as such n an automated
More informationUncertainty quantification via codimension one domain partitioning and a new concentration inequality
Uncertainty quantification via codimension one domain partitioning and a new concentration inequality Tim Sullivan, Ufuk Topcu, Mike McKerns & Houman Owhadi tjs@..., utopcu@..., mmckerns@...& owhadi@caltech.edu
More informationDistributionally Robust Convex Optimization
Submitted to Operations Research manuscript OPRE-2013-02-060 Authors are encouraged to submit new papers to INFORMS journals by means of a style file template, which includes the journal title. However,
More informationDistributionally Robust Convex Optimization
Distributionally Robust Convex Optimization Wolfram Wiesemann 1, Daniel Kuhn 1, and Melvyn Sim 2 1 Department of Computing, Imperial College London, United Kingdom 2 Department of Decision Sciences, National
More informationData-Driven Optimization under Distributional Uncertainty
Data-Driven Optimization under Distributional Uncertainty Postdoctoral Researcher Electrical and Systems Engineering University of Pennsylvania A network of physical objects - Devices - Vehicles - Buildings
More informationarxiv: v1 [math.oc] 27 Nov 2013
Submitted, SIAM Journal on Optimization (SIOPT) http://www.cds.caltech.edu/~murray/papers/han+13-siopt.html Convex Optimal Uncertainty Quantification arxiv:1311.7130v1 [math.oc] 27 Nov 2013 Shuo Han 1,
More informationBrittleness and Robustness of Bayesian Inference
Brittleness and Robustness of Bayesian Inference Tim Sullivan with Houman Owhadi and Clint Scovel (Caltech) Mathematics Institute, University of Warwick Predictive Modelling Seminar University of Warwick,
More informationHandout 8: Dealing with Data Uncertainty
MFE 5100: Optimization 2015 16 First Term Handout 8: Dealing with Data Uncertainty Instructor: Anthony Man Cho So December 1, 2015 1 Introduction Conic linear programming CLP, and in particular, semidefinite
More informationRandom Convex Approximations of Ambiguous Chance Constrained Programs
Random Convex Approximations of Ambiguous Chance Constrained Programs Shih-Hao Tseng Eilyan Bitar Ao Tang Abstract We investigate an approach to the approximation of ambiguous chance constrained programs
More informationRobust Optimization for Risk Control in Enterprise-wide Optimization
Robust Optimization for Risk Control in Enterprise-wide Optimization Juan Pablo Vielma Department of Industrial Engineering University of Pittsburgh EWO Seminar, 011 Pittsburgh, PA Uncertainty in Optimization
More informationDistributionally robust optimization techniques in batch bayesian optimisation
Distributionally robust optimization techniques in batch bayesian optimisation Nikitas Rontsis June 13, 2016 1 Introduction This report is concerned with performing batch bayesian optimization of an unknown
More informationDistributionally Robust Stochastic Optimization with Wasserstein Distance
Distributionally Robust Stochastic Optimization with Wasserstein Distance Rui Gao DOS Seminar, Oct 2016 Joint work with Anton Kleywegt School of Industrial and Systems Engineering Georgia Tech What is
More informationHandout 6: Some Applications of Conic Linear Programming
ENGG 550: Foundations of Optimization 08 9 First Term Handout 6: Some Applications of Conic Linear Programming Instructor: Anthony Man Cho So November, 08 Introduction Conic linear programming CLP, and
More informationQuantifying Stochastic Model Errors via Robust Optimization
Quantifying Stochastic Model Errors via Robust Optimization IPAM Workshop on Uncertainty Quantification for Multiscale Stochastic Systems and Applications Jan 19, 2016 Henry Lam Industrial & Operations
More informationDistributionally Robust Optimization with Infinitely Constrained Ambiguity Sets
manuscript Distributionally Robust Optimization with Infinitely Constrained Ambiguity Sets Zhi Chen Imperial College Business School, Imperial College London zhi.chen@imperial.ac.uk Melvyn Sim Department
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization
Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012
More informationDistributionally robust simple integer recourse
Distributionally robust simple integer recourse Weijun Xie 1 and Shabbir Ahmed 2 1 Department of Industrial and Systems Engineering, Virginia Tech, Blacksburg, VA 24061 2 School of Industrial & Systems
More informationWorst-case design of structures using stopping rules in k-adaptive random sampling approach
10 th World Congress on Structural and Multidisciplinary Optimization May 19-4, 013, Orlando, Florida, USA Worst-case design of structures using stopping rules in k-adaptive random sampling approach Makoto
More informationarxiv: v3 [math.pr] 23 May 2012
Optimal Uncertainty Quantification H. Owhadi, T. J. Sullivan, M. McKerns, M. Ortiz California Institute of Technology owhadi@caltech.edu arxiv:1009.0679v3 [math.pr] 23 May 2012 C. Scovel Los Alamos National
More informationDistributionally Robust Optimization with ROME (part 1)
Distributionally Robust Optimization with ROME (part 1) Joel Goh Melvyn Sim Department of Decision Sciences NUS Business School, Singapore 18 Jun 2009 NUS Business School Guest Lecture J. Goh, M. Sim (NUS)
More informationV&V MURI Overview Caltech, October 2008
V&V MURI Overview Caltech, October 2008 Pablo A. Parrilo Laboratory for Information and Decision Systems Massachusetts Institute of Technology Goals!! Specification, design, and certification!! Coherent
More informationM E M O R A N D U M. Faculty Senate approved November 1, 2018
M E M O R A N D U M Faculty Senate approved November 1, 2018 TO: FROM: Deans and Chairs Becky Bitter, Sr. Assistant Registrar DATE: October 23, 2018 SUBJECT: Minor Change Bulletin No. 5 The courses listed
More informationBrittleness and Robustness of Bayesian Inference
Brittleness and Robustness of Bayesian Inference Tim Sullivan 1 with Houman Owhadi 2 and Clint Scovel 2 1 Free University of Berlin / Zuse Institute Berlin, Germany 2 California Institute of Technology,
More informationModel-Based Rigorous Uncertainty Quantification in Complex Systems
Model-Based Rigorous Uncertainty Quantification in Complex Systems California Institute of Technology MECOM 2010 Buenos Aires, November 17, 2010 Work done in collaboration with: Marc Adams, Addis Kadani,
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationAmbiguity in portfolio optimization
May/June 2006 Introduction: Risk and Ambiguity Frank Knight Risk, Uncertainty and Profit (1920) Risk: the decision-maker can assign mathematical probabilities to random phenomena Uncertainty: randomness
More informationConvergence Rates of Kernel Quadrature Rules
Convergence Rates of Kernel Quadrature Rules Francis Bach INRIA - Ecole Normale Supérieure, Paris, France ÉCOLE NORMALE SUPÉRIEURE NIPS workshop on probabilistic integration - Dec. 2015 Outline Introduction
More informationCS295: Convex Optimization. Xiaohui Xie Department of Computer Science University of California, Irvine
CS295: Convex Optimization Xiaohui Xie Department of Computer Science University of California, Irvine Course information Prerequisites: multivariate calculus and linear algebra Textbook: Convex Optimization
More informationApproximation algorithms for nonnegative polynomial optimization problems over unit spheres
Front. Math. China 2017, 12(6): 1409 1426 https://doi.org/10.1007/s11464-017-0644-1 Approximation algorithms for nonnegative polynomial optimization problems over unit spheres Xinzhen ZHANG 1, Guanglu
More informationDistributionally Robust Discrete Optimization with Entropic Value-at-Risk
Distributionally Robust Discrete Optimization with Entropic Value-at-Risk Daniel Zhuoyu Long Department of SEEM, The Chinese University of Hong Kong, zylong@se.cuhk.edu.hk Jin Qi NUS Business School, National
More informationA Robust von Neumann Minimax Theorem for Zero-Sum Games under Bounded Payoff Uncertainty
A Robust von Neumann Minimax Theorem for Zero-Sum Games under Bounded Payoff Uncertainty V. Jeyakumar, G.Y. Li and G. M. Lee Revised Version: January 20, 2011 Abstract The celebrated von Neumann minimax
More informationConditional and Dynamic Preferences
Conditional and Dynamic Preferences How can we Understand Risk in a Dynamic Setting? Samuel Drapeau Joint work with Hans Föllmer Humboldt University Berlin Jena - March 17th 2009 Samuel Drapeau Joint work
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationQuadratic Two-Stage Stochastic Optimization with Coherent Measures of Risk
Noname manuscript No. (will be inserted by the editor) Quadratic Two-Stage Stochastic Optimization with Coherent Measures of Risk Jie Sun Li-Zhi Liao Brian Rodrigues Received: date / Accepted: date Abstract
More informationCOMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO
COMPARISON PRINCIPLES FOR CONSTRAINED SUBHARMONICS PH.D. COURSE - SPRING 2019 UNIVERSITÀ DI MILANO KEVIN R. PAYNE 1. Introduction Constant coefficient differential inequalities and inclusions, constraint
More informationRobust control and applications in economic theory
Robust control and applications in economic theory In honour of Professor Emeritus Grigoris Kalogeropoulos on the occasion of his retirement A. N. Yannacopoulos Department of Statistics AUEB 24 May 2013
More informationRobust linear optimization under general norms
Operations Research Letters 3 (004) 50 56 Operations Research Letters www.elsevier.com/locate/dsw Robust linear optimization under general norms Dimitris Bertsimas a; ;, Dessislava Pachamanova b, Melvyn
More informationVariational Inequalities. Anna Nagurney Isenberg School of Management University of Massachusetts Amherst, MA 01003
Variational Inequalities Anna Nagurney Isenberg School of Management University of Massachusetts Amherst, MA 01003 c 2002 Background Equilibrium is a central concept in numerous disciplines including economics,
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1110.0918ec e-companion ONLY AVAILABLE IN ELECTRONIC FORM informs 2011 INFORMS Electronic Companion Mixed Zero-One Linear Programs Under Objective Uncertainty: A Completely
More informationLecture: Examples of LP, SOCP and SDP
1/34 Lecture: Examples of LP, SOCP and SDP Zaiwen Wen Beijing International Center For Mathematical Research Peking University http://bicmr.pku.edu.cn/~wenzw/bigdata2018.html wenzw@pku.edu.cn Acknowledgement:
More informationRobust Farkas Lemma for Uncertain Linear Systems with Applications
Robust Farkas Lemma for Uncertain Linear Systems with Applications V. Jeyakumar and G. Li Revised Version: July 8, 2010 Abstract We present a robust Farkas lemma, which provides a new generalization of
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
More informationOn Distributionally Robust Chance Constrained Program with Wasserstein Distance
On Distributionally Robust Chance Constrained Program with Wasserstein Distance Weijun Xie 1 1 Department of Industrial and Systems Engineering, Virginia Tech, Blacksburg, VA 4061 June 15, 018 Abstract
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationarxiv: v1 [cs.dm] 27 Jan 2014
Randomized Minmax Regret for Combinatorial Optimization Under Uncertainty Andrew Mastin, Patrick Jaillet, Sang Chin arxiv:1401.7043v1 [cs.dm] 27 Jan 2014 January 29, 2014 Abstract The minmax regret problem
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationarxiv: v1 [math.oc] 18 Jul 2011
arxiv:1107.3493v1 [math.oc] 18 Jul 011 Tchebycheff systems and extremal problems for generalized moments: a brief survey Iosif Pinelis Department of Mathematical Sciences Michigan Technological University
More informationRobust Growth-Optimal Portfolios
Robust Growth-Optimal Portfolios! Daniel Kuhn! Chair of Risk Analytics and Optimization École Polytechnique Fédérale de Lausanne rao.epfl.ch 4 Technology Stocks I 4 technology companies: Intel, Cisco,
More informationOptimal Uncertainty Quantification of Quantiles
Optimal Uncertainty Quantification of Quantiles Application to Industrial Risk Management. Merlin Keller 1, Jérôme Stenger 12 1 EDF R&D, France 2 Institut Mathématique de Toulouse, France ETICS 2018, Roscoff
More information1 PROBLEM DEFINITION. i=1 z i = 1 }.
Algorithms for Approximations of Nash Equilibrium (003; Lipton, Markakis, Mehta, 006; Kontogiannis, Panagopoulou, Spirakis, and 006; Daskalakis, Mehta, Papadimitriou) Spyros C. Kontogiannis University
More informationData-Driven Distributionally Robust Chance-Constrained Optimization with Wasserstein Metric
Data-Driven Distributionally Robust Chance-Constrained Optimization with asserstein Metric Ran Ji Department of System Engineering and Operations Research, George Mason University, rji2@gmu.edu; Miguel
More informationDistirbutional robustness, regularizing variance, and adversaries
Distirbutional robustness, regularizing variance, and adversaries John Duchi Based on joint work with Hongseok Namkoong and Aman Sinha Stanford University November 2017 Motivation We do not want machine-learned
More informationReformulation of chance constrained problems using penalty functions
Reformulation of chance constrained problems using penalty functions Martin Branda Charles University in Prague Faculty of Mathematics and Physics EURO XXIV July 11-14, 2010, Lisbon Martin Branda (MFF
More informationIntroduction to Convex Optimization
Introduction to Convex Optimization Daniel P. Palomar Hong Kong University of Science and Technology (HKUST) ELEC5470 - Convex Optimization Fall 2018-19, HKUST, Hong Kong Outline of Lecture Optimization
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
More informationSecond-order cone programming formulation for two player zero-sum game with chance constraints
Second-order cone programming formulation for two player zero-sum game with chance constraints Vikas Vikram Singh a, Abdel Lisser a a Laboratoire de Recherche en Informatique Université Paris Sud, 91405
More informationChance constrained optimization - applications, properties and numerical issues
Chance constrained optimization - applications, properties and numerical issues Dr. Abebe Geletu Ilmenau University of Technology Department of Simulation and Optimal Processes (SOP) May 31, 2012 This
More informationRisk-Averse Dynamic Optimization. Andrzej Ruszczyński. Research supported by the NSF award CMMI
Research supported by the NSF award CMMI-0965689 Outline 1 Risk-Averse Preferences 2 Coherent Risk Measures 3 Dynamic Risk Measurement 4 Markov Risk Measures 5 Risk-Averse Control Problems 6 Solution Methods
More informationOptimal Uncertainty Quantification
SIAM REVIEW Vol. 55, No. 2, pp. 271 345 c 2013 Society for Industrial and Applied Mathematics Optimal Uncertainty Quantification H. Owhadi C. Scovel T. J. Sullivan M. McKerns M. Ortiz Abstract. We propose
More informationEstimation and Optimization: Gaps and Bridges. MURI Meeting June 20, Laurent El Ghaoui. UC Berkeley EECS
MURI Meeting June 20, 2001 Estimation and Optimization: Gaps and Bridges Laurent El Ghaoui EECS UC Berkeley 1 goals currently, estimation (of model parameters) and optimization (of decision variables)
More informationModel-Based Rigorous Uncertainty Quantification in Complex Systems
Model-Based Rigorous Uncertainty Quantification in Complex Systems California Institute of Technology Congress on Numerical Methods in Engineering (CMNE 2011) University of Coimbra Coimbra, Portugal, June
More informationFast Direct Methods for Gaussian Processes
Fast Direct Methods for Gaussian Processes Mike O Neil Departments of Mathematics New York University oneil@cims.nyu.edu December 12, 2015 1 Collaborators This is joint work with: Siva Ambikasaran Dan
More informationHyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems. Joanna Zwierzchowska
Decision Making in Manufacturing and Services Vol. 9 05 No. pp. 89 00 Hyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems Joanna Zwierzchowska Abstract.
More informationRobust portfolio selection under norm uncertainty
Wang and Cheng Journal of Inequalities and Applications (2016) 2016:164 DOI 10.1186/s13660-016-1102-4 R E S E A R C H Open Access Robust portfolio selection under norm uncertainty Lei Wang 1 and Xi Cheng
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationWavelet-Based Numerical Homogenization for Scaled Solutions of Linear Matrix Equations
International Journal of Discrete Mathematics 2017; 2(1: 10-16 http://www.sciencepublishinggroup.com/j/dmath doi: 10.11648/j.dmath.20170201.13 Wavelet-Based Numerical Homogenization for Scaled Solutions
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
More informationA note on scenario reduction for two-stage stochastic programs
A note on scenario reduction for two-stage stochastic programs Holger Heitsch a and Werner Römisch a a Humboldt-University Berlin, Institute of Mathematics, 199 Berlin, Germany We extend earlier work on
More informationapproximation algorithms I
SUM-OF-SQUARES method and approximation algorithms I David Steurer Cornell Cargese Workshop, 201 meta-task encoded as low-degree polynomial in R x example: f(x) = i,j n w ij x i x j 2 given: functions
More informationTime inconsistency of optimal policies of distributionally robust inventory models
Time inconsistency of optimal policies of distributionally robust inventory models Alexander Shapiro Linwei Xin Abstract In this paper, we investigate optimal policies of distributionally robust (risk
More informationRandom sets. Distributions, capacities and their applications. Ilya Molchanov. University of Bern, Switzerland
Random sets Distributions, capacities and their applications Ilya Molchanov University of Bern, Switzerland Molchanov Random sets - Lecture 1. Winter School Sandbjerg, Jan 2007 1 E = R d ) Definitions
More informationOptimal Transport Methods in Operations Research and Statistics
Optimal Transport Methods in Operations Research and Statistics Jose Blanchet (based on work with F. He, Y. Kang, K. Murthy, F. Zhang). Stanford University (Management Science and Engineering), and Columbia
More informationA projection algorithm for strictly monotone linear complementarity problems.
A projection algorithm for strictly monotone linear complementarity problems. Erik Zawadzki Department of Computer Science epz@cs.cmu.edu Geoffrey J. Gordon Machine Learning Department ggordon@cs.cmu.edu
More informationCanonical Problem Forms. Ryan Tibshirani Convex Optimization
Canonical Problem Forms Ryan Tibshirani Convex Optimization 10-725 Last time: optimization basics Optimization terology (e.g., criterion, constraints, feasible points, solutions) Properties and first-order
More informationNumerical Solutions of ODEs by Gaussian (Kalman) Filtering
Numerical Solutions of ODEs by Gaussian (Kalman) Filtering Hans Kersting joint work with Michael Schober, Philipp Hennig, Tim Sullivan and Han C. Lie SIAM CSE, Atlanta March 1, 2017 Emmy Noether Group
More informationNested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model
Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Xiaowei Chen International Business School, Nankai University, Tianjin 371, China School of Finance, Nankai
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationSignal Recovery from Permuted Observations
EE381V Course Project Signal Recovery from Permuted Observations 1 Problem Shanshan Wu (sw33323) May 8th, 2015 We start with the following problem: let s R n be an unknown n-dimensional real-valued signal,
More informationCORC Technical Report TR Ambiguous chance constrained problems and robust optimization
CORC Technical Report TR-2004-0 Ambiguous chance constrained problems and robust optimization E. Erdoğan G. Iyengar September 7, 2004 Abstract In this paper we study ambiguous chance constrained problems
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationDecomposition Algorithms for Two-Stage Distributionally Robust Mixed Binary Programs
Decomposition Algorithms for Two-Stage Distributionally Robust Mixed Binary Programs Manish Bansal Grado Department of Industrial and Systems Engineering, Virginia Tech Email: bansal@vt.edu Kuo-Ling Huang
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationNear-Potential Games: Geometry and Dynamics
Near-Potential Games: Geometry and Dynamics Ozan Candogan, Asuman Ozdaglar and Pablo A. Parrilo January 29, 2012 Abstract Potential games are a special class of games for which many adaptive user dynamics
More informationExperiences with Model Reduction and Interpolation
Experiences with Model Reduction and Interpolation Paul Constantine Stanford University, Sandia National Laboratories Qiqi Wang (MIT) David Gleich (Purdue) Emory University March 7, 2012 Joe Ruthruff (SNL)
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationSupermodular ordering of Poisson arrays
Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore
More informationAN INEQUALITY FOR TAIL PROBABILITIES OF MARTINGALES WITH BOUNDED DIFFERENCES
Lithuanian Mathematical Journal, Vol. 4, No. 3, 00 AN INEQUALITY FOR TAIL PROBABILITIES OF MARTINGALES WITH BOUNDED DIFFERENCES V. Bentkus Vilnius Institute of Mathematics and Informatics, Akademijos 4,
More informationStochastic optimal control with rough paths
Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction
More informationOn the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. xx, No. x, Xxxxxxx 00x, pp. xxx xxx ISSN 0364-765X EISSN 156-5471 0x xx0x 0xxx informs DOI 10.187/moor.xxxx.xxxx c 00x INFORMS On the Power of Robust Solutions in
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationStrong Duality in Robust Semi-Definite Linear Programming under Data Uncertainty
Strong Duality in Robust Semi-Definite Linear Programming under Data Uncertainty V. Jeyakumar and G. Y. Li March 1, 2012 Abstract This paper develops the deterministic approach to duality for semi-definite
More informationLecture 1. Stochastic Optimization: Introduction. January 8, 2018
Lecture 1 Stochastic Optimization: Introduction January 8, 2018 Optimization Concerned with mininmization/maximization of mathematical functions Often subject to constraints Euler (1707-1783): Nothing
More informationA direct formulation for sparse PCA using semidefinite programming
A direct formulation for sparse PCA using semidefinite programming A. d Aspremont, L. El Ghaoui, M. Jordan, G. Lanckriet ORFE, Princeton University & EECS, U.C. Berkeley Available online at www.princeton.edu/~aspremon
More informationUsing what we know: Inference with Physical Constraints
Using what we know: Inference with Physical Constraints Chad M. Schafer & Philip B. Stark Department of Statistics University of California, Berkeley {cschafer,stark}@stat.berkeley.edu PhyStat 2003 8 September
More informationMaster of Science in Statistics A Proposal
1 Master of Science in Statistics A Proposal Rationale of the Program In order to cope up with the emerging complexity on the solutions of realistic problems involving several phenomena of nature it is
More information