Convex Optimization for Simplicial Mesh Improvement
|
|
- Matilda Harper
- 6 years ago
- Views:
Transcription
1 Convex Optimization for Simplicial Mesh Improvement Bryan Feldman December 9, 25 Abstract In this project I explore several strategies utilizing convex optimization for the purpose of improving the quality of a simplicial mesh. The intended application for the meshes is for performing finite element/volume simulations and therefore mesh quality and in particular the minimum element quality is of fundamental importance. In the first part I look at methods for creating valid meshes from ones that contain inverted elements. In the second part I look at improving the minimum element quality of the mesh. 2 Introduction Creating good quality meshes is important because poorly shaped elements can severely hinder the performance of finite element or finite volume simulations that use these meshes to discretize the simulation domain. A large number of researchers have investigated this problem using various strategies and within the literature several definitions for quality have arisen. The work of Schewcuk, [4] has shed some light on the quality of these various quality metrics and in what situations (if any) the metrics are appropriate. Quality can be improved in two ways, by optimizing the topology of elements or by optimizing the positions of the nodes. Delaunay triangulation (2D) or tetrahedralization (3D) is known to minimize several mesh quality metrics, and is widely regarded as the defacto standard method for generation of the underlying mesh connectivity. In this work I focus on the problem of optimizing the location of nodes in a mesh, referred to as mesh smoothing. Within the field of mesh smoothing there are generally two approaches, heuristic-based and optimization-based. The heuristic based methods are often relatively quick and simple. The most common method, Laplacian Smoothing, simply adjusts a node position to the average of its neighbors. However, heuristic methods offer no guarantees to their performance. In fact, these methods can decrease mesh quality even to the point of creating invalid meshes. Alternatively, optimization-based approaches, while slower, generally create nicer meshes as they are guided by the objective one wishes to maximize. Since the number of nodes in a
2 mesh is typically quite large, almost exclusively these algorithms operate by optimizing the mesh quality by moving a single node at a time, and repeating for all (or some subset of) nodes, perhaps several times. 3 Mesh Untangling In order to have any hope of a simulation producing reasonable results the mesh that the simulation is performed on need be valid, i.e. not contain any inverted or tangled components. One can envision obtaining an invalid mesh in the situation where mesh nodes are moved over a time step of the simulation without changing mesh topology such that some elements of the mesh become inverted. A quick indication for an inverted element is if the area (2D) or volume (3D) of the element is negative. In both 2D and 3D the area/volume can be computed by the determinant of a matrix formed by 2/3 edges of the element and so it is linear with respect to a single node position. As show in [2] a mesh can be untangled by optimizing the position nodes in the mesh belonging to a tangled elements. Since the optimization is performed for a single node at a time, several sweeps through the mesh may be required to untangle all nodes but this cost is substantially less than forming and optimizing the problem as a bundle. The area of the triangles surrounding a particular node, i at position, x i are the components of a vector, a which is found by a = A i x i + c Here obviously A i contains the linear multiples in the area calculations and c has all the required constant components. In [2], x i is optimized for an objective function that maximizes the minimum area max min j N a j This can be achieved by forming a matrix A c that is the concatenation of A i with a row of all negative ones, and solving for a vector, π which is the concatenation of x i and the minimum area (volume) local element. We can define slack variables s, as the difference between each local elements area and the minimum area (the last component of π) which from definitions yields: A c π + c = s With this we can maximize the minimum area element by the following LP maximize b T π subject to A c π + c = s s where b is a vector of all zeros except for the last term which is one. Since the last term of π is the minimum element area, maximizing this quantity achieves the desired optimization function. The constraints ensure that the difference between all element areas and the minimum is. 2
3 Additionally, I tested a method for untangling meshes that minimizes the two norm between the vector of local areas and some desired areas, â. Desired areas can be found by setting them to the average of existing local triangles, resulting in the following optimization, minimize A i x i + c â 2 where â = N ( T a ). This minimization can easily be solved as it is a standard Least-Squared problem. 3. Results The above described algorithms were implemented in MATLAB using the Convex Optimization package cvx As shown in Figure the methods used to untangle the meshes all successfully created valid meshes. They differ in the run time and ending quality of the meshes. There is apparently an inverse relationship between complexity of the algorithm and quality/speed of obtaining the results for these 2-D examples. Laplacian Smoothing, the dumbest and by far the fastest was able eventually untangle all test meshes I gave it. While it sometimes inverted a triangle subsequent applications of smoothing to either adjacent nodes or the same node in the next pass, seemed to undo the problem and resulted in a rather nice looking mesh. The algorithm which found the node position that made triangles have near equal area performed slightly worse (and slower) while the hard-core, maximize the minimum area algorithm created a rather poor quality mesh, although to be fair it did do the required job of untangling the mesh. I should point out that all tests are performed on 2-D meshes which the literature reports are substantially easier to handle than their notoriously fickle 3-D counterparts. Therefore, it stands to reason that in 3-D the more sophisticated algorithms may be required. This endeavor remains as an area for future work. 4 Node Position Optimization While the work in the previous section was able to successfully untangle initial meshes that included many severely inverted triangles, the quality of the triangles are often far from ideal. This shouldn t be surprising as area is not a good indicator of mesh quality, so optimizations involving only the area have no preferred tendency to be high quality (just not inverted). Thus, in this next part of the project I investigate a method for optimizing a sensible mesh quality metric. As reported in [4] there are a plethora of mesh quality metrics that have been used. One particularly nice metric suggested in the above mentioned work is q i = 4 3 A i l 2 rms boyd/cvx/ (2D) q i = V i A rms l rms (3D) 3
4 .5 tangled nodepos Freitag nodepos Laplacian smoothing nodepos avg Area nodepos Figure : Comparison of methods to untangle mesh. Top Left) original tangled mesh. Top Right)Untangling by maximizing minimum area. Lower Left) Laplacian Smoothing to untangle. Lower Right) Minimizing distance from average area. where q i is the quality of element i, l is length of an edge, A is area, V volume and the subscript rms indicates the root mean squared of all edges, or areas belonging to that element. Also note that these metrics include normalization factors such that the ideal (equilateral) elements has a quality of one. This metric is convex, smooth, and differentiable and actually has a nice intuitive meaning as being the condition number of the matric which transforms the element to the ideal element. While overall (or average) mesh quality is important, the accuracy of a simulation can be ruined by a single bad element as it has the capability to ill-condition matrices. Therefore, it is important to create meshes that have elements of at least some minimal quality. Using ideas presented in Optimization Strategies in Unstructured Mesh Generation [5] I implemented a system to maximize an objective function, Ψ(x), the minimum quality element in the mesh, Ψ(x) = min i=...n (q i) 4
5 with my implementation using a different quality metric than the one in [5]. In this form the problem is particularly nasty for a number of reasons. First, the problem size can be quite large, meshes can contain hundreds of thousands or even millions of nodes, rendering any optimization intractable. The second nastiness of this problem is the fact that the objective is non-differentiable. In order to cope with these problems [5] suggest several relaxations to the problem. In order to reduce the dimension of the problem, first the current worst quality element is found and only the positions of the nodes in some local neighborhood are used in the optimization, x. In order to deal with the non-differentiability of the objective they suggest a modified Steepest Decent Algorithm that is taken from [3]. The approach relaxes an indicator variable, I, to be a real valued vector, µ. The steps of the optimization are:. Compute the step direction δ = m µ i q i (x ) where µ is computed as the solution of Quadratic Programing sub-problem i µ = argmin [ m i µ k [Ψ(x ) q i (x )]] 2 + L 2 m i µ k q i (x ) 2 s.t. µ T µ = with L set to an average local edge length size. 2. Compute the step size, λ using Armijo s rule λ = max λ = {β k (k =,,...) Ψ(x + λδ) Ψ(x ) λα δ 2 } 3. Take step (adjust positions of local neighborhood of nodes) x = x + λδ Because solving the QP sub-problem is expensive another relaxation is used to reduce the dimension of the problem. It is assumed that µ is zero for all local elements with q(x ) Ψ(x ) ɛ with a typical value for ɛ =.. This reduces the size of the QP substantially, potentially even removing the need for it if no other elements have low quality. This results in a significant speedup of the optimization. 4. Results The algorithm successfully improved many of the poor quality elements. As shown in Figure 2 the resulting mesh looks substantially better than the initial one, although in fact some of the worst elements in the initial mesh are difficult to see since they are so thin. Also in Figure 3 we see that the minimum element quality is increased from.78 to.443. Simulation on 5
6 the initial mesh would probably create problems while the improved mesh should do fine. Comparing the optimized mesh with the one obtained using Laplacian Smoothing, again we see that that Laplacian Smoothing results in a superior mesh. This brings to mind a quip from John Strain of our math department who said, You should resist the urge to be clever. In the case of 2-D meshes it seems his message rings true. In the defense of cleverness I must again point out that the examples take place in 2-D and that complexity in 3-D is so much more that there is probably a place for some cleverness when dealing with 3-D meshes. In particular, a good system for mesh improvement would probably include both the simple and robust method. One can envision using Laplacian Smoothing to quickly improve large portions of the mesh, and then relying on the more robust approaches of optimization-based methods to deal with problems remaining/created by the heuristic-based methods. before after Figure 2: Mesh before and after optimization by method described in Section 4 References [] L. Freitag and P. Knupp. Tetrahedral mesh improvement of the element condition number. Inter. Journal for Numerical Methods in Engineering, 53:377 39, 22. [2] L. Freitag and P. Plassmann. Local optimization-based simplicial mesh untangling and improvement. 2nd Symposium on Trends in Unstructured Mesh Generation,
7 min mean Figure 3: Improvement in minimum and mean quality [3] Elijah Polak. On the mathematical foundations of nondifferentiable optimization in engineering design. Technical Report UCB/ERL M85/7, EECS Department, University of California, Berkeley, 985. [4] J.R. Shewchuk. What is a good linear element? interpolation, conditioning and quality measures. th International Meshing Roundtable, 22. [5] P.D. Zavattieri, E.A. Dari, and G.C. Buscaglia. Optimization strategies in unstructured mesh generation. Int. Journal of Numerical Methods, 39:255 27,
Lecture 2 Unstructured Mesh Generation
Lecture 2 Unstructured Mesh Generation MIT 16.930 Advanced Topics in Numerical Methods for Partial Differential Equations Per-Olof Persson (persson@mit.edu) February 13, 2006 1 Mesh Generation Given a
More informationSolving for dynamic user equilibrium
Solving for dynamic user equilibrium CE 392D Overall DTA problem 1. Calculate route travel times 2. Find shortest paths 3. Adjust route choices toward equilibrium We can envision each of these steps as
More informationCSE 417 Dynamic Programming (pt 4) Sub-problems on Trees
CSE 417 Dynamic Programming (pt 4) Sub-problems on Trees Reminders > HW4 is due today > HW5 will be posted shortly Dynamic Programming Review > Apply the steps... 1. Describe solution in terms of solution
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationSegmentation and Grouping
Segmentation and Grouping How and what do we see? Fundamental Problems ' Focus of attention, or grouping ' What subsets of pixels do we consider as possible objects? ' All connected subsets? ' Representation
More informationInterval velocity estimation through convex optimization
Stanford Exploration Project, Report 125, January 16, 2007, pages 1?? Interval velocity estimation through convex optimization Ben Witten and Michael Grant ABSTRACT Convex optimization is an optimization
More informationWeighted and Continuous Clustering
John (ARC/ICAM) Virginia Tech... Math/CS 4414: http://people.sc.fsu.edu/ jburkardt/presentations/ clustering weighted.pdf... ARC: Advanced Research Computing ICAM: Interdisciplinary Center for Applied
More informationSimplicial Global Optimization
Simplicial Global Optimization Julius Žilinskas Vilnius University, Lithuania September, 7 http://web.vu.lt/mii/j.zilinskas Global optimization Find f = min x A f (x) and x A, f (x ) = f, where A R n.
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationUnderstanding Gridfit
Understanding Gridfit John R. D Errico Email: woodchips@rochester.rr.com December 28, 2006 1 Introduction GRIDFIT is a surface modeling tool, fitting a surface of the form z(x, y) to scattered (or regular)
More information4 Integer Linear Programming (ILP)
TDA6/DIT37 DISCRETE OPTIMIZATION 17 PERIOD 3 WEEK III 4 Integer Linear Programg (ILP) 14 An integer linear program, ILP for short, has the same form as a linear program (LP). The only difference is that
More informationGeneration of Triangle Meshes from Time-of-Flight Data for Surface Registration
Generation of Triangle Meshes from Time-of-Flight Data for Surface Registration Thomas Kilgus, Thiago R. dos Santos, Alexander Seitel, Kwong Yung, Alfred M. Franz, Anja Groch, Ivo Wolf, Hans-Peter Meinzer,
More informationA Patch-Based Mesh Optimization Algorithm for Partitioned Meshes
A Patch-Based Mesh Optimization Algorithm for Partitioned Meshes Nicholas Voshell 1, Suzanne Shontz 1, Lori Diachin 2, Patrick Knupp 3, and Todd Munson 4 1 Department of Computer Science and Engineering
More informationChapter 2 Basic Structure of High-Dimensional Spaces
Chapter 2 Basic Structure of High-Dimensional Spaces Data is naturally represented geometrically by associating each record with a point in the space spanned by the attributes. This idea, although simple,
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationParallel Unstructured Mesh Generation by an Advancing Front Method
MASCOT04-IMACS/ISGG Workshop University of Florence, Italy Parallel Unstructured Mesh Generation by an Advancing Front Method Yasushi Ito, Alan M. Shih, Anil K. Erukala, and Bharat K. Soni Dept. of Mechanical
More informationAs a consequence of the operation, there are new incidences between edges and triangles that did not exist in K; see Figure II.9.
II.4 Surface Simplification 37 II.4 Surface Simplification In applications it is often necessary to simplify the data or its representation. One reason is measurement noise, which we would like to eliminate,
More informationExploiting a database to predict the in-flight stability of the F-16
Exploiting a database to predict the in-flight stability of the F-16 David Amsallem and Julien Cortial December 12, 2008 1 Introduction Among the critical phenomena that have to be taken into account when
More informationPARALLEL OPTIMIZATION OF TETRAHEDRAL MESHES
6 th European Conference on Computational Mechanics (ECCM 6) 7 th European Conference on Computational Fluid Dynamics (ECFD 7) 11... 15 June 2018, Glasgow, UK PARALLEL OPTIMIZATION OF TETRAHEDRAL MESHES
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationSome Open Problems in Graph Theory and Computational Geometry
Some Open Problems in Graph Theory and Computational Geometry David Eppstein Univ. of California, Irvine Dept. of Information and Computer Science ICS 269, January 25, 2002 Two Models of Algorithms Research
More informationA Brief Overview of Optimization Problems. Steven G. Johnson MIT course , Fall 2008
A Brief Overview of Optimization Problems Steven G. Johnson MIT course 18.335, Fall 2008 Why optimization? In some sense, all engineering design is optimization: choosing design parameters to improve some
More informationGEOMETRY MODELING & GRID GENERATION
GEOMETRY MODELING & GRID GENERATION Dr.D.Prakash Senior Assistant Professor School of Mechanical Engineering SASTRA University, Thanjavur OBJECTIVE The objectives of this discussion are to relate experiences
More informationRecent Developments in Model-based Derivative-free Optimization
Recent Developments in Model-based Derivative-free Optimization Seppo Pulkkinen April 23, 2010 Introduction Problem definition The problem we are considering is a nonlinear optimization problem with constraints:
More informationA New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes
A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes Sanjay Kumar Khattri Department of Mathematics, University of Bergen, Norway sanjay@mi.uib.no http://www.mi.uib.no/ sanjay Abstract. Mesh
More informationTopological Classification of Data Sets without an Explicit Metric
Topological Classification of Data Sets without an Explicit Metric Tim Harrington, Andrew Tausz and Guillaume Troianowski December 10, 2008 A contemporary problem in data analysis is understanding the
More informationUntangling and Smoothing of Quadrilateral and Hexahedral Meshes
Paper 36 Untangling and Smoothing of Quadrilateral and Hexahedral Meshes Civil-Comp Press, 2012 Proceedings of the Eighth International Conference on Engineering Computational Technology, B.H.V. Topping,
More informationSome Applications of Graph Bandwidth to Constraint Satisfaction Problems
Some Applications of Graph Bandwidth to Constraint Satisfaction Problems Ramin Zabih Computer Science Department Stanford University Stanford, California 94305 Abstract Bandwidth is a fundamental concept
More informationLecture 11: Clustering and the Spectral Partitioning Algorithm A note on randomized algorithm, Unbiased estimates
CSE 51: Design and Analysis of Algorithms I Spring 016 Lecture 11: Clustering and the Spectral Partitioning Algorithm Lecturer: Shayan Oveis Gharan May nd Scribe: Yueqi Sheng Disclaimer: These notes have
More informationOutline of the presentation
Surface Reconstruction Petra Surynková Charles University in Prague Faculty of Mathematics and Physics petra.surynkova@mff.cuni.cz Outline of the presentation My work up to now Surfaces of Building Practice
More informationEdge Detection (with a sidelight introduction to linear, associative operators). Images
Images (we will, eventually, come back to imaging geometry. But, now that we know how images come from the world, we will examine operations on images). Edge Detection (with a sidelight introduction to
More informationA Brief Overview of Optimization Problems. Steven G. Johnson MIT course , Fall 2008
A Brief Overview of Optimization Problems Steven G. Johnson MIT course 18.335, Fall 2008 Why optimization? In some sense, all engineering design is optimization: choosing design parameters to improve some
More informationNote Set 4: Finite Mixture Models and the EM Algorithm
Note Set 4: Finite Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine Finite Mixture Models A finite mixture model with K components, for
More informationHow and what do we see? Segmentation and Grouping. Fundamental Problems. Polyhedral objects. Reducing the combinatorics of pose estimation
Segmentation and Grouping Fundamental Problems ' Focus of attention, or grouping ' What subsets of piels do we consider as possible objects? ' All connected subsets? ' Representation ' How do we model
More informationAutomatic hybrid mesh generation for the boundary face method
Boundary Elements and Other Mesh Reduction Methods XXXVI 139 Automatic hybrid mesh generation for the boundary face method Cheng Huang & Jianming Zhang State Key Laboratory of Advanced Design and Manufacturing
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
More informationLarge Mesh Deformation Using the Volumetric Graph Laplacian
Large Mesh Deformation Using the Volumetric Graph Laplacian Kun Zhou1 Jin Huang2 John Snyder3 Xinguo Liu1 Hujun Bao2 Baining Guo1 Heung-Yeung Shum1 1 Microsoft Research Asia 2 Zhejiang University 3 Microsoft
More informationA Course in Machine Learning
A Course in Machine Learning Hal Daumé III 13 UNSUPERVISED LEARNING If you have access to labeled training data, you know what to do. This is the supervised setting, in which you have a teacher telling
More informationEdge and corner detection
Edge and corner detection Prof. Stricker Doz. G. Bleser Computer Vision: Object and People Tracking Goals Where is the information in an image? How is an object characterized? How can I find measurements
More informationLecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1
CME 305: Discrete Mathematics and Algorithms Instructor: Professor Aaron Sidford (sidford@stanfordedu) February 6, 2018 Lecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1 In the
More informationMathematical Morphology and Distance Transforms. Robin Strand
Mathematical Morphology and Distance Transforms Robin Strand robin.strand@it.uu.se Morphology Form and structure Mathematical framework used for: Pre-processing Noise filtering, shape simplification,...
More informationLecture 27: Fast Laplacian Solvers
Lecture 27: Fast Laplacian Solvers Scribed by Eric Lee, Eston Schweickart, Chengrun Yang November 21, 2017 1 How Fast Laplacian Solvers Work We want to solve Lx = b with L being a Laplacian matrix. Recall
More informationRobust Realignment of fmri Time Series Data
Robust Realignment of fmri Time Series Data Ben Dodson bjdodson@stanford.edu Olafur Gudmundsson olafurg@stanford.edu December 12, 2008 Abstract FMRI data has become an increasingly popular source for exploring
More informationChapter 1. Numeric Artifacts. 1.1 Introduction
Chapter 1 Numeric Artifacts 1.1 Introduction Virtually all solutions to problems in electromagnetics require the use of a computer. Even when an analytic or closed form solution is available which is nominally
More informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More information05 - Surfaces. Acknowledgements: Olga Sorkine-Hornung. CSCI-GA Geometric Modeling - Daniele Panozzo
05 - Surfaces Acknowledgements: Olga Sorkine-Hornung Reminder Curves Turning Number Theorem Continuous world Discrete world k: Curvature is scale dependent is scale-independent Discrete Curvature Integrated
More informationA new method of quality improvement for quadrilateral mesh based on small polygon reconnection
Acta Mech. Sin. (2012) 28(1):140 145 DOI 10.1007/s10409-012-0022-x RESEARCH PAPER A new method of quality improvement for quadrilateral mesh based on small polygon reconnection Jian-Fei Liu Shu-Li Sun
More informationK-Means Clustering Using Localized Histogram Analysis
K-Means Clustering Using Localized Histogram Analysis Michael Bryson University of South Carolina, Department of Computer Science Columbia, SC brysonm@cse.sc.edu Abstract. The first step required for many
More informationTELCOM2125: Network Science and Analysis
School of Information Sciences University of Pittsburgh TELCOM2125: Network Science and Analysis Konstantinos Pelechrinis Spring 2015 2 Part 4: Dividing Networks into Clusters The problem l Graph partitioning
More informationSurfaces, meshes, and topology
Surfaces from Point Samples Surfaces, meshes, and topology A surface is a 2-manifold embedded in 3- dimensional Euclidean space Such surfaces are often approximated by triangle meshes 2 1 Triangle mesh
More informationLAPLACIAN MESH SMOOTHING FOR TETRAHEDRA BASED VOLUME VISUALIZATION 1. INTRODUCTION
JOURNAL OF MEDICAL INFORMATICS & TECHNOLOGIES Vol.4/2002, ISSN 642-6037 Rafał STĘGIERSKI *, Paweł MIKOŁAJCZAK * volume data,triangle mesh generation, mesh smoothing, marching tetrahedra LAPLACIAN MESH
More informationSmoothing and local refinement techniques for improving tetrahedral mesh quality
Smoothing and local refinement techniques for improving tetrahedral mesh quality J.M. Escobar, R. Montenegro, G. Montero, E. Rodríguez, J.M. González-Yuste University Institute of Intelligent Systems and
More informationTowards a unified measurement of quality for mixed elements
Towards a unified measurement of qualit for mixed elements technical report n 015/01, di utfsm Claudio Lobos Departmento de Informática Universidad Técnica Federico Santa María Santiago, Chile clobos@inf.utfsm.cl
More informationInteractive segmentation, Combinatorial optimization. Filip Malmberg
Interactive segmentation, Combinatorial optimization Filip Malmberg But first... Implementing graph-based algorithms Even if we have formulated an algorithm on a general graphs, we do not neccesarily have
More informationA Data Dependent Triangulation for Vector Fields
A Data Dependent Triangulation for Vector Fields Gerik Scheuermann Hans Hagen Institut for Computer Graphics and CAGD Department of Computer Science University of Kaiserslautern, Postfach 3049, D-67653
More informationSurface Mesh Generation
Surface Mesh Generation J.-F. Remacle Université catholique de Louvain September 22, 2011 0 3D Model For the description of the mesh generation process, let us consider the CAD model of a propeller presented
More informationShape Modeling and Geometry Processing
252-0538-00L, Spring 2018 Shape Modeling and Geometry Processing Discrete Differential Geometry Differential Geometry Motivation Formalize geometric properties of shapes Roi Poranne # 2 Differential Geometry
More informationAnnouncements. Edges. Last Lecture. Gradients: Numerical Derivatives f(x) Edge Detection, Lines. Intro Computer Vision. CSE 152 Lecture 10
Announcements Assignment 2 due Tuesday, May 4. Edge Detection, Lines Midterm: Thursday, May 6. Introduction to Computer Vision CSE 152 Lecture 10 Edges Last Lecture 1. Object boundaries 2. Surface normal
More informationDigital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering
Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing
More informationCS 523: Computer Graphics, Spring Shape Modeling. Differential Geometry of Surfaces
CS 523: Computer Graphics, Spring 2011 Shape Modeling Differential Geometry of Surfaces Andrew Nealen, Rutgers, 2011 2/22/2011 Differential Geometry of Surfaces Continuous and Discrete Motivation Smoothness
More informationAdvanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material
More informationMöbius Transformations in Scientific Computing. David Eppstein
Möbius Transformations in Scientific Computing David Eppstein Univ. of California, Irvine School of Information and Computer Science (including joint work with Marshall Bern from WADS 01 and SODA 03) Outline
More informationLecture 2 September 3
EE 381V: Large Scale Optimization Fall 2012 Lecture 2 September 3 Lecturer: Caramanis & Sanghavi Scribe: Hongbo Si, Qiaoyang Ye 2.1 Overview of the last Lecture The focus of the last lecture was to give
More informationScientific Computing: Interpolation
Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25
More informationReporting Mesh Statistics
Chapter 15. Reporting Mesh Statistics The quality of a mesh is determined more effectively by looking at various statistics, such as maximum skewness, rather than just performing a visual inspection. Unlike
More informationApproximate Linear Programming for Average-Cost Dynamic Programming
Approximate Linear Programming for Average-Cost Dynamic Programming Daniela Pucci de Farias IBM Almaden Research Center 65 Harry Road, San Jose, CA 51 pucci@mitedu Benjamin Van Roy Department of Management
More informationMAT128A: Numerical Analysis Lecture Two: Finite Precision Arithmetic
MAT128A: Numerical Analysis Lecture Two: Finite Precision Arithmetic September 28, 2018 Lecture 1 September 28, 2018 1 / 25 Floating point arithmetic Computers use finite strings of binary digits to represent
More informationStructural and Syntactic Pattern Recognition
Structural and Syntactic Pattern Recognition Selim Aksoy Department of Computer Engineering Bilkent University saksoy@cs.bilkent.edu.tr CS 551, Fall 2017 CS 551, Fall 2017 c 2017, Selim Aksoy (Bilkent
More informationSpectral Surface Reconstruction from Noisy Point Clouds
Spectral Surface Reconstruction from Noisy Point Clouds 1. Briefly summarize the paper s contributions. Does it address a new problem? Does it present a new approach? Does it show new types of results?
More informationComparing different interpolation methods on two-dimensional test functions
Comparing different interpolation methods on two-dimensional test functions Thomas Mühlenstädt, Sonja Kuhnt May 28, 2009 Keywords: Interpolation, computer experiment, Kriging, Kernel interpolation, Thin
More informationSOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS
SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented
More informationABOUT THE GENERATION OF UNSTRUCTURED MESH FAMILIES FOR GRID CONVERGENCE ASSESSMENT BY MIXED MESHES
VI International Conference on Adaptive Modeling and Simulation ADMOS 2013 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ABOUT THE GENERATION OF UNSTRUCTURED MESH FAMILIES FOR GRID CONVERGENCE
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationHIGH-ORDER MESH GENERATION ON CAD GEOMETRIES
VI International Conference on Adaptive Modeling and Simulation ADMOS 2013 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) HIGH-ORDER MESH GENERATION ON CAD GEOMETRIES A. GARGALLO-PEIRÓ,
More informationClassical Gradient Methods
Classical Gradient Methods Note simultaneous course at AMSI (math) summer school: Nonlin. Optimization Methods (see http://wwwmaths.anu.edu.au/events/amsiss05/) Recommended textbook (Springer Verlag, 1999):
More informationFinal Project, Digital Geometry Processing
Final Project, Digital Geometry Processing Shayan Hoshyari Student #: 81382153 December 2016 Introduction In this project an adaptive surface remesher has been developed based on the paper [1]. An algorithm
More information16.410/413 Principles of Autonomy and Decision Making
16.410/413 Principles of Autonomy and Decision Making Lecture 17: The Simplex Method Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology November 10, 2010 Frazzoli (MIT)
More informationComputational Methods. Constrained Optimization
Computational Methods Constrained Optimization Manfred Huber 2010 1 Constrained Optimization Unconstrained Optimization finds a minimum of a function under the assumption that the parameters can take on
More informationSamuel Coolidge, Dan Simon, Dennis Shasha, Technical Report NYU/CIMS/TR
Detecting Missing and Spurious Edges in Large, Dense Networks Using Parallel Computing Samuel Coolidge, sam.r.coolidge@gmail.com Dan Simon, des480@nyu.edu Dennis Shasha, shasha@cims.nyu.edu Technical Report
More informationA Constrained Delaunay Triangle Mesh Method for Three-Dimensional Unstructured Boundary Point Cloud
International Journal of Computer Systems (ISSN: 2394-1065), Volume 03 Issue 02, February, 2016 Available at http://www.ijcsonline.com/ A Constrained Delaunay Triangle Mesh Method for Three-Dimensional
More informationSimplicial Complexes: Second Lecture
Simplicial Complexes: Second Lecture 4 Nov, 2010 1 Overview Today we have two main goals: Prove that every continuous map between triangulable spaces can be approximated by a simplicial map. To do this,
More informationLevel Set Extraction from Gridded 2D and 3D Data
Level Set Extraction from Gridded 2D and 3D Data David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationShape fitting and non convex data analysis
Shape fitting and non convex data analysis Petra Surynková, Zbyněk Šír Faculty of Mathematics and Physics, Charles University in Prague Sokolovská 83, 186 7 Praha 8, Czech Republic email: petra.surynkova@mff.cuni.cz,
More informationConformal Flattening ITK Filter
=1 Conformal Flattening ITK Filter Release 0.00 Yi Gao 1, John Melonakos 1, and Allen Tannenbaum 1 July 10, 2006 1 Georgia Institute of Technology, Atlanta, GA Abstract This paper describes the Insight
More informationNative mesh ordering with Scotch 4.0
Native mesh ordering with Scotch 4.0 François Pellegrini INRIA Futurs Project ScAlApplix pelegrin@labri.fr Abstract. Sparse matrix reordering is a key issue for the the efficient factorization of sparse
More informationAM205: lecture 2. 1 These have been shifted to MD 323 for the rest of the semester.
AM205: lecture 2 Luna and Gary will hold a Python tutorial on Wednesday in 60 Oxford Street, Room 330 Assignment 1 will be posted this week Chris will hold office hours on Thursday (1:30pm 3:30pm, Pierce
More informationQ. Wang National Key Laboratory of Antenna and Microwave Technology Xidian University No. 2 South Taiba Road, Xi an, Shaanxi , P. R.
Progress In Electromagnetics Research Letters, Vol. 9, 29 38, 2009 AN IMPROVED ALGORITHM FOR MATRIX BANDWIDTH AND PROFILE REDUCTION IN FINITE ELEMENT ANALYSIS Q. Wang National Key Laboratory of Antenna
More informationLP-Modelling. dr.ir. C.A.J. Hurkens Technische Universiteit Eindhoven. January 30, 2008
LP-Modelling dr.ir. C.A.J. Hurkens Technische Universiteit Eindhoven January 30, 2008 1 Linear and Integer Programming After a brief check with the backgrounds of the participants it seems that the following
More informationI. Meshing and Accuracy Settings
Guidelines to Set CST Solver Accuracy and Mesh Parameter Settings to Improve Simulation Results with the Time Domain Solver and Hexahedral Meshing System illustrated with a finite length horizontal dipole
More informationSummary of Raptor Codes
Summary of Raptor Codes Tracey Ho October 29, 2003 1 Introduction This summary gives an overview of Raptor Codes, the latest class of codes proposed for reliable multicast in the Digital Fountain model.
More informationUsing ɛ-dominance for Hidden and Degenerated Pareto-Fronts
IEEE Symposium Series on Computational Intelligence Using ɛ-dominance for Hidden and Degenerated Pareto-Fronts Heiner Zille Institute of Knowledge and Language Engineering University of Magdeburg, Germany
More informationA surface mesh smoothing and untangling method independent of the CAD parameterization
A surface mesh smoothing and untangling method independent of the CAD parameterization The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationUsing a Divide and Conquer Method for Routing in a PC Vehicle Routing Application. Abstract
Using a Divide and Conquer Method for Routing in a PC Vehicle Routing Application Brenda Cheang Department of Management Information Systems University College Dublin Belfield, Dublin 4, Ireland. Sherlyn
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationReal-Time Shape Editing using Radial Basis Functions
Real-Time Shape Editing using Radial Basis Functions, Leif Kobbelt RWTH Aachen Boundary Constraint Modeling Prescribe irregular constraints Vertex positions Constrained energy minimization Optimal fairness
More informationNetworks: A Deterministic Approach with Constant Overhead
Trace-Routing in 3D Wireless Sensor Networks: A Deterministic Approach with Constant Overhead Su Xia, Hongyi Wu, and Miao Jin! School of Computing and Informatics! University of Louisiana at Lafayette
More informationPOLYHEDRAL GEOMETRY. Convex functions and sets. Mathematical Programming Niels Lauritzen Recall that a subset C R n is convex if
POLYHEDRAL GEOMETRY Mathematical Programming Niels Lauritzen 7.9.2007 Convex functions and sets Recall that a subset C R n is convex if {λx + (1 λ)y 0 λ 1} C for every x, y C and 0 λ 1. A function f :
More informationProject Updates Short lecture Volumetric Modeling +2 papers
Volumetric Modeling Schedule (tentative) Feb 20 Feb 27 Mar 5 Introduction Lecture: Geometry, Camera Model, Calibration Lecture: Features, Tracking/Matching Mar 12 Mar 19 Mar 26 Apr 2 Apr 9 Apr 16 Apr 23
More informationA Global Laplacian Smoothing Approach with Feature Preservation
A Global Laplacian Smoothing Approach with Feature Preservation hongping Ji Ligang Liu Guojin Wang Department of Mathematics State Key Lab of CAD&CG hejiang University Hangzhou, 310027 P.R. China jzpboy@yahoo.com.cn,
More informationSubdivision Surfaces. Course Syllabus. Course Syllabus. Modeling. Equivalence of Representations. 3D Object Representations
Subdivision Surfaces Adam Finkelstein Princeton University COS 426, Spring 2003 Course Syllabus I. Image processing II. Rendering III. Modeling IV. Animation Image Processing (Rusty Coleman, CS426, Fall99)
More information