p 1 q 1 p 3 q 3 q 2 p 2
|
|
- Evangeline Cooper
- 5 years ago
- Views:
Transcription
1 Automatic Grid Generation with Almost regular Delaunay Tetrahedra Alexander Fuchs Department of Mathematics A University of Stuttgart Pfaenwaldring 57, D Stuttgart, Germany fuchs@mathematik.uni-stuttgart.de Abstract. An algorithm for constructing almost regular triangulations (ARTs) for three-dimensional polygonal domains is described. Ideally such triangulations consist entirely of congruent tetrahedra with nearly equal edges. The new feature of this method is that the position of the vertices is adjusted, before any connecting edges are assigned. This leads to grids with very few irregular vertices, i.e. most interior vertices are shared by 24 tetrahedra as for regular tessellations of R 3 with congruent tetrahedra. keywords. delaunay tetrahedra, tessellation of R 3, almost regularity Introduction Grid generation algorithms are essential tools in many areas of modern applied mathematics. Typical examples are nite element methods (FEM, [2]) and geometric modelling (CAGD, [4]). In both applications the quality of the grids strongly inuences the accuracy of numerical computations. For example, irregularly shaped tetrahedra result in large condition numbers of the stiness matrices for Galerkin methods. Similarly, spline approximation on irregular partitions requires non standard smoothness conditions [8, 5, 4, 5]. Ideally, a triangulation should consist entirely of equilateral tetrahedra. Unfortunately, this is not possible. However, as will be shown in section 2, there exist congruent tetrahedra having nearly equal edges. These tessellations have a regular lattice structure, i.e. every vertex has 4 neighbours and is shared by 24 tetrahedra. While such regularity can seldom be achieved it is natural to insist on \almost regularity". An almost regular triangulation (ART) should have only few irregular vertices, i.e. most of the interior vertices are shared by exactly 24 tetrahedra. Such combinatorial optimization also leads to good geometric properties. Clearly for an ART smoothing of the vertices usually results in very small deviations of the angles from the optimal value. Figure shows an example. The triangulation consists of 860 vertices and has just 83 irregular vertices. The percentage of angles less than 2:0 is lower than 0:2% and the maximum angle is 68:2. Ohand the construction of ARTs seems rather dicult since the combinatorial requirements may result in global constraints on the tetrahedra. We therefore employ an indirect approach, and extend our 2D-algorithm (see [0]) to 3D. We rst construct a set of inner vertices according to a density function %. This function is dened over the entire domain and describes the edge-length of the triangulation. In the second step we optimize the position of the vertices without assigning any connecting edges. Then, we form a Delaunay triangulation. More precisely, our algorithm for constructing an ART of a polygonal domain consists of three steps: (ART) construct an initial conguration and scale the density function, (ART2) adjust the position of the vertices by minimizing a penalty functional, (ART3) form a Delaunay triangulation. As usual, the density function controls the size of the tetrahedra which may vary throughout the domain. We do not discuss its construction here, which may depend on the geometry as well as on the local error of nite element calculations. The functional in the second step is minimal for regular triangulations. The optimization will therefore yield an almost regular conguration of the vertices. After these preprocessing steps the Delaunay triangulation gives excellent results, although the geometric advantages of two-dimensional Delaunay ([9],[7]) do not hold for three dimensions. supported by SFB 404, University of Stuttgart
2 Figure : Example of an ART In sections 3{5 the steps of the ART-Algorithm and the underlying theory will be described in detail. In section 6 we discuss a number of examples illustrating the performance of our method. We conclude with some remarks on possible extensions. 2 Regular Tessellations with Tetrahedra A regular triangulation is a partition of R 3 into congruent tetrahedra T i, so that their intersection is either empty or a vertex, edge or face. Unlike in two dimensions, there exists no canonical regular triangulation since one cannot partition R 3 with equilateral tetrahedra. Of course for applications in nite element methods and geometric modelling, the tetrahedra should be nearly equilateral. Another natural requirement is invariance under subdivision. If a tetrahedron T is split into 8 tetrahedra by halving the edges as is indicated in Figure 2, the small tetrahedra should by congruent to T scaled by the factor. Subdivision invariance is particularly important for grid renement. 2 Moreover, it is easy to show that a subdivision invariant tetrahedron generates a regular triangulation of R 3. One should note, however, that neither an equilateral tetrahedron nor the one produced by the tessellation of the standard cube (corners (0,0,0),(,0,0), (,,0) and (,,)) are subdivision invariant. p 4 q 6 q 4 q 5 p 3 T : p ; q ; q 3; q 4 T 2 : p 2; q ; q 2; q 5 T 3 : p 3; q 2; q 3; q 6 T 4 : p 4; q 4; q 5; q 6 T 5 : q ; q 2; q 4; q 5 T 6 : q ; q 2; q 4; q 3 T 7 : q 2; q 4; q 6; q 5 T 8 : q 2; q 4; q 6; q 3 q 3 q 2 p q p 2 Figure 2: Subdivision of a tetrahedron T.
3 With the following theorem we characterize all subdivision invariant tetrahedra. Theorem A tetrahedron in standard coordinates with vertices p = 0; p 2 = v = 0 0 A ; p3 = v 2 = x y 0 0 A and p4 = v 3 x2 y 2 z having x ; x 2; y 2 0 and y ; z > 0 (cf. Figure 2), is subdivision invariant, i one of the following 4 conditions is fullled: A ; Condition : 0 x < 2 x 2 = 2x y 2 =? x 2 y 2 ( + x)(? 2x)2 2 =? x Condition 2: y 2 > 0 z 2 ( + x)(? 2x) =? x x = 3 y = 2y 2 x 2 = 2 3 z = p 3y 2 Condition 3: 0 x < 2 x 2 =? x y 2 =? x 2 y 2 2 = x2 ( + x )? x Condition 4: 0 < x < 3 z 2 ( + x)(? 2x) =? x x 2 = 2 (x + ) y2 = x (3? x ) y2 2 = x(3? x) 4 z2 = (3? x) 4 Proof. It is easy to show that tetrahedra T, T 2, T 3 and T 4 are congruent to T. In the interior of T there are 2 2 dierent tetrahedra which have to be mapped to the tetrahedron T by suitable Euclidean movements. This results in the conditions to 4. The explicit computation of the congruence mappings is straightforward but rather tedious. Since the arguments are not crucial for the subsequent development of the grid generation algorithm we refer to [9] v v 2 v 3 Figure 3: Tetrahedra forming a parallelepiped. Figure 4: Tetrahedron T. Any of the tetrahedra dened by Theorem gives rise to a regular triangulation. Analyzing the structure of such triangulations in more detail, we rst notice that 6 tetrahedra form a parallelepiped as shown in Figure 3. The labels at the vertices in Figure 3 refer to the number of tetrahedra sharing this vertex. Moreover, we see that 24 tetrahedra meet in a regular vertex of a partition with parallelepipeds and each regular vertex has 4 neighbours. As a prototype for our grid generation algorithm we choose a subdivision invariant tetrahedron which diers from an equilateral tetrahedron as little as possible. Figure 5 shows the standard deviation of the dihedral angles of the tetrahedra satisfying the conditions of Theorem as a function of the free parameters. We select a tetrahedron T satisfying condition with x = =3 (see Figure 5).
4 Condition x x Condition Condition x Figure 5: Standard deviation of dihedral angles T has 3 congruent faces with angles 54:74 ; 54:74 ; 70:52 and the standard deviation to 60 is 7:44. The lengths of the edges of the tetrahedron T are ; ; ; ; 2= p 3 :5; 2= p 3 with a standard deviation of (2? p 3)=3 0:089 from the equilateral case, and the dihedral angles are 60 ; 60 ; 60 ; 60 ; 90 ; 90, having a standard deviation of 4:5 from 70:53, the dihedral angle for the equilateral tetrahedron. 3 Construction of an initial Conguration We can obtain an initial conguration with the following simple algorithm. We choose an approximate centre C 0 of the domain and construct the polyhedron formed by the 24 tetrahedra sharing C 0 (see Figure 6). Now we scale the polyhedron, until it covers the entire domain. Then we divide the tetrahedra which, according to the density function, are too large and repeat this, until all tetrahedra are small enough. For a constant density function, we get a fully regular triangulation. Finally we project points outside the domain, which have connection to an inner point, on the boundary of. All points within and on the boundary of the domain are then used as an initial set of vertices. At this point the exact position of the points on the boundary is not crucial, since no connecting edges are assigned in the subsequent optimization procedure. The number of vertices obtained in this way will, in general, not exactly match the density function. Since the density function plays an important role in the following optimization process (cf. Section 4), we have to rescale the density
5 Figure 6: Polyhedron, formed by 24 tetrahedra (solid and upper half transparent) function %! s % (s ) with an appropriate scaling factor s, in order to match the number of the constructed vertices best. Therefore, we have to analyse the relation between the density function and the number of vertices. The polygonal domain to be triangulated can be multiply connected and is specied by the edges and vertices making up the boundary polygons. For example the domain in Figure 7 has 2 boundary polyhedrons with a total of 24 edges and 6 vertices. Figure 7: Polygonal Domain To estimate the number of interior vertices X j, j = ; : : : ; n i in terms of the density function %, we assume that the triangulation is regular, i.e. that 24 tetrahedra share an interior and 2 tetrahedra share a at boundary vertex. The number of tetrahedra sharing a vertex on a boundary edge varies as can be seen in Figure 3. But in general all \types" of edge vertices occur with the same frequency, so we assume that 6 tetrahedra share an edge vertex. Denoting by n i, n b and n e the number of interior, boundary and edge vertices of the triangulation, it follows that the total number of tetrahedra is n t = (24n i + 2n b + 6n e)=4: () From Figure 3 it can be seen that 6 triangles share a boundary vertex and an edge vertex. It follows that the number of triangles on the boundary of the domain is n f = (6n b + 6n e)=3: (2) Hence we obtain for the number of vertices n = n i + n b + n e the formula n = 6 nt + 4 n f + ne: (3) 4
6 Now we relate n t and n f to the density function %, dropping our assumption that the triangulation is regular. The length of the edges should be approximately equal to the value of % at their midpoints. If we denote by V 0 the volume of the tetrahedron T (cf. Section 2), a regular triangulation with constant density function yields n t V 0 % 3 = vol(); and we can compute n t in terms of the volume of the domain. For non constant % we use the approximation Z % 3 X vol() % 3 X V 0 = V 0n t; (4) where the sums are taken over all tetrahedra of the triangulation and % denotes the average value of % on. In order to relate n f to the density function, we note that there may exist boundary triangles with dierent areas, depending on the tetrahedron T chosen (for our choice all boundary triangles have the same area). Since in general all triangles should occur with roughly the same frequency we have % 2 F0n f ; (5) where F 0 is the average area of the boundary triangles. Similarly we estimate the number of edge vertices by denoting by E 0 the average edge length of the tetrahedron T. E0ne; (6) % Combining equations (3), (4), (5) and (6) we obtain for the number of vertices Z Z Z n 6V 0 % + 3 4F 0 % + 2 4E 0 % : (7) We can now use equation (7) to scale the density function, %! s %, so that equation (7) is maintained for the number of vertices constructed with our algorithm. This estimate was derived for regular triangulations only. However, it remains fairly accurate for an ART, in particular when the number of irregular vertices is small. In any case, a slight deviation from the optimal number of points has no major impact on the quality of the triangulation as will become apparent from the next section. Merely the length of the edges will not correspond exactly to the density and the triangulation might be slightly distorted at 4 Adjusting the Vertices Adjusting the vertices to improve the geometric characteristics of a triangulation is a standard procedure [, 6]. Unlike to smoothing in previous work (cf. e.g. [7, 8, 6]), the new feature of our approach is to we optimize the position of the vertices, before assigning any connections to them. There associated Delaunay triangulation should be almost regular. This can be achieved by minimizing a functional f which penalizes congurations producing irregular vertices. The choice of our functional is based on the following conjecture. A vector H is a permutation of D (H; D 2 R n ) if the equations (h + + h n) = = (d + + d n) = hold. (h h 2 n) =2 = (d d 2 n) =2. (8) (h n + + h n n) =n = (d n + + d n n) =n A proof for this assumption has yet not been found, but we can prove the following stronger theorem.
7 Theorem 2 A vector H is a permutation of a vector D i the equations hold. (h + + h n) = = (d + + d n) = (h h 2 n) =2 = (d d 2 n) =2. (9) (h 2n + + h 2n n ) =(2n) = (d 2n + + d 2n n ) =(2n) To proof this, we need the following lemma. Lemma Equations (9) imply for every polynom p(x) of degree 2n p(h ) + p(h 2) + + p(h n) = p(d ) + p(d 2) + + p(d n) : (0) Proof. Let p(x) = a 2nx 2n + a 2n?x 2n? + + a 0 be the monomial representation of the polynom, then the following holds p(h ) + + p(h n) and thus (0) holds. = a 2n(h 2n + + h 2n n ) + a n?(h 2n? + + h 2n? n ) + + n a 0 = a 2n (h 2n + + h 2n n ) =(2n) 2n + an? (h 2n? + + h 2n? n ) =(2n?) 2n? + + n a0 (because of (9)) = a 2n (d 2n + + d 2n n ) =(2n) 2n + an? (d 2n? + + d 2n? n ) =(2n?) 2n? + + n a0 = a 2n(d 2n + + d 2n n ) + a n?(d 2n? + + d 2n? n ) + + n a 0 = p(d ) + + p(d n) Having this lemma, the proof of Theorem 2 is quite simple. Proof of Theorem 2. \)" : If H is a permutation of D, then the equations (9) are trivially fullled. \(" : (complete induction) For n = we have h = d, and thus H a permutation of D. Now use Lemma with the polynom p(x) := (x? d ) 2 (x? d 2) 2 (x? d n) 2 : Then for a solution of equations (9) we have p(h ) + + p(h n) = p(d ) + + p(d n) = 0 : Because p(x) 0 it has to be p(h i) = 0 and thus 8h i9d j : h i = d j. Without loss of generality let now be h n = d n, then we can reduce the dimension of (9) by one. Now we can dene a functional, which penalizes the error of equations (8). We denote by d d 2 the distances of the vertices of a regular triangulation to a xed vertex. For an almost regular triangulation, the corresponding distances kx i? X k k; k = ; 2; : : : ; n should be close to a permutation of D. This means that 0 =j 0 B kx i? X k C B C k j? n k= k6=i n X k= k6=i d j ka=j is small for each i and j. This motivates the following denition of our penalty functional f:
8 f(x ; ; X n) := nx i= j= 00 mx n k= k6=i '(h i;k C ) j A =j? 0 n k= k6=i '(d k C ) j A =j2 C A ; h i;k := kx i? X k k (%(X i) + %(X k ))=2 : We have replaced kx i? X k k by the normalized distance h i;k which is a good approximation to the corresponding component of D for neighbouring points since (%(X i) + %(X j))=2 describes the average side length; h i;k is approximately equal to the number of tetrahedra between the i-th and k-th vertex. For a constant density %, which results in an uniformal spacing of the vertices, this replacement is just a normalization of the distances. If we have a non constant density function, which results in a non-uniformal spacing of the vertices, then h i;k normalizes the distances of neighbouring vertices in the lattice. Of course, small distances are most relevant. Therefore, we introduce a cut-o function ' which vanishes for arguments larger than d l. d l+ d l d l d l+ d l d l+ Figure 8: Obvious but poor choice of '. Figure 9: Good choice of '. The obvious choice of ' is shown in Figure 8. The drawback of this choice is, that we cannot distinguish a distance in the interval [d l d l+ ] and [d l 0]. Therefore we prefer another (monotone) choice of ' as shown in Figure 9. We notice that we only have to compute l + distances in vector D, because '(x) = 0 for x d l+. In practice it suces to take one layer of neighbouring vertices, so that l = 4. For our choice of tetrahedron T we have D 2 = [d ; ; d l+ ] = [ 4=3 4=3 4=3 4=3 4=3 4=3 8=3]. Because of the cut-o function '() most of the terms in equations (8) are zero. Therefore it makes sense not to take the sum over j = ; ; n and to introduce a cut-o parameter m, having l m n. In practice we get good results for m = 2 l. One notices that, because of the local character of the 'cut-o' function ', the sum over k will contain only few nonzero terms. Therefore, with an appropriate data structure, the function f can be computed in O(n) operations. One also has to take into account, that the vector D diers for vertices inside and at the boundary of, thus the vertices X i have to be split into 2 sets. Of course there is great deal of heuristic in our choice of the penalty functional f. Strictly speaking, our analysis is valid only for a regular grid which is compatible with the boundary. However, because of the local character of ', one can expect that f will be nearly minimal for ARTs if the density function is smooth. In fact, in all numerical tests (cf. Section 6) we have obtained very good results. The numerical implementation of the minimization procedure is straightforward from a mathematical point of view. We employ a conjugate gradient method, and project the \search-direction" for vertices on the boundary. The actual C++-Code ( lines) is rather complex. More generally than described in this paper, the algorithm can handle free an xed vertices on the boundary as well as internal constraints. A number of additional strategies are employed to improve the eciency of the program. To speed up the convergence of the iteration we can compute
9 rst a coarse triangulation and obtain improved start congurations via subdivision. Also, we can parallelize most of the computations. Figure 0 shows a few typical steps of the smoothing process. Start conguration Point conguration after one step of CG 5 steps of CG 0 steps of CG Figure 0: Smoothing Process 5 Triangulation A Delaunay tessellation can be obtained as dual graph of the Voronoi-diagram, which consists of the sets of closest points to a given set of vertices [9]. To determine the Delaunay tetrahedra we use the following criterion. Denition (Delaunay Tetrahedron) The vertices P,Q,R,S form a Delaunay tetrahedron, if there is no further vertex in its circumscribed ball. Methods for a sequential construction of the Delaunay tetrahedra sequentially are well known (cf. e.g. [9]). Since our vertices are constructed with respect to a density function, we can implement an ecient parallel algorithm. This algorithm will be described in detail in []. In order to keep the resulting triangulation conform to the boundary, we use a constrained Delaunay triangulation (cf. e.g. [3]). Of course, the placement of the edges is not essential on smooth parts of the boundary, where any polygon is just an approximation. For the example considered in the previous section, Figures and 2 show the resulting Delaunay triangulation. For complex domains or domains with a lot of interior faces we can construct a density function, which \recognizes" the local complexity of the domain, and adapts the density. This has already been done for two dimensions (cf. Figure 3) and will be described for three dimensions in a forthcoming paper.
10 Figure : View from above Figure 2: Cut through the domain 2d-net zoomed net Figure 3: Two dimensional net with interior edges and automatically adapted density 6 Examples In the following we discuss various examples, illustrating the performance of our algorithm. For evaluating the quality of the results we use the following criteria which compare the ARTs with regular triangulations. Irregular vertices By n irr we denote the number of interior vertices which do not have 4 neighbours. For an ART, n irr should be small compared to n i, the total number of interior vertices. Dihedral angles To emphasize the distribution of the dihedral angles, we give the minimum ( min) and maximum ( max) angle, as well as the percentage of angles lower than 6 ; 2 and bigger than 68 ; 74.
11 Condition of FEM matrices We denote by C the condition number of the stiness matrix for the Galerkin approximation of the Poisson problem using piecewise linear basis functions. This condition number is compared with the corresponding condition numbers C o and C i for regular triangulations of a subset of a bounding domain. More precisely, we choose for C o the tetrahedra which have an intersection with the domain (cf. Figure 4) and for C i those which lie fully inside the domain (cf. Figure 5). The edge length is chosen as the average edge length of the triangulation. Figure 4: Outbounding Net Figure 5: Inbounding Net Figure 6 shows the triangulation of a cube with a hole. Since for this example the percentage of boundary vertices is rather high, there is not much exibility to position the interior vertices. Hence, there are many irregular vertices. Nevertheless, the minimum and maximum angle is very good ( min = 7:6, max = 68:2 ). Also the condition number for the FE-matrix is comparable to the one for regular meshes. Vertices Angle % Dihedral angles Condition Total Inner Boundary Irreg. min max < 6 < 2 > 68 > 74 C C o C i :6 68: Figure 6: Cube with a Hole
12 Figure 7 shows the triangulation of a ball. Surprisingly the triangulation consists entirely of regular vertices. Vertices Angle % Dihedral angles Condition Total Inner Boundary Irreg. min max < 6 < 2 > 68 > 74 C C o C i :4 48: Figure 7: Triangulation of a ball with cross section In Figure 9 we show a constraint triangulation. Since for this case most of the interior vertices have to be positioned on the inner plane, the number of irregular vertices is rather high. Nevertheless, the minimum ( min = 6:9 ) and maximum ( max = 69:3 ) dihedral angle is acceptable and Figure 8 shows the distribution of the inner vertices and angles for this triangulation n inner vertices angles Figure 8: Distribution of inner vertices and angles.
13 Vertices Angle % Dihedral angles Condition Total Inner Boundary Irreg. min max < 6 < 2 > 68 > 74 C C o C i :9 69: Figure 9: Cube with an inner plane
14 The last example (Figure 20) shows a graded triangulation. Vertices Angle % Dihedral angles Condition Total Inner Boundary Irreg. min max < 6 < 2 > 68 > 74 C C o C i :9 74: Figure 20: Graded Triangulation 7 Conclusion We have presented a method for constructing almost regular triangulations (ARTs). The key idea is to adjust the vertices of Delaunay triangulations using a functional which penalizes irregular vertices. Our actual implementation of the triangulation algorithm is slightly more exible than described in this paper. For example, one may choose features like interior faces and one has a number of options to control the optimization process. For such details we refer to [2] and the online description of our program which is available via anonymous ftp from ftp.mathematik.uni-stuttgart.de or via WWW from A corresponding algorithm for 2-dimensional problems is described in [0]. An advantage of our method compared to many existing algorithms is, that the smoothing and optimization of the point conguration is done before any connections are assigned. This provides greater exibility, in particular when propagating irregularities of the boundary towards the interior of the domain. We intend to generalize our algorithm to NURBS solids and trimmed surfaces. Acknowledgements This work was performed at the University of Stuttgart. This work was supported by the SFB 404, University of Stuttgart. References [] Scott A. Canann, Micheal B. Stephenson, and Ted Blacker. Optismoothing: An optimization-driven approach to mesh smoothing. Finite Elements in Analysis and Design, 3:85{90, 993.
15 [2] Susanne C. Brenner and Scott L. Ridgway. The Mathematical Theory of Finite Element Methods. Springer- Verlag, 994. [3] H. Borouchaki, F. Hecht, E. Saltel and P.L. George. Reasonably Ecient Delaunay Based Mesh Generator in 3 Dimensions. In Proceedings, 4th International Meshing Roundtable, pages 3{4. Sandia National Laboratories, 995. [4] G. Farin. Geometric Modelling: Algorithms and new Trends. SIAM, 987. [5] G. Farin. Curves and Surfaces for Computer Aided Geometric Design. Academic Press Inc. (London) Ltd., edition, 988. [6] Lori A. Freitag, Mark Jones, and Paul Plassmann. An Ecient Parallel Algorithm for Mesh Smoothing. In Proceedings, 4th International Meshing Roundtable, pages 47{58. Sandia National Laboratories, 995. [7] Lori A. Freitag and Carl Ollivier-Gooch. A Comparison of Tetrahedral Mesh Improvement Techniques. In Proceedings, 5th International Meshing Roundtable, pages 87{00. Sandia National Laboratories, 996. [8] Lori A. Freitag. Tetrahedral Mesh Improvement Using Swapping and Smoothing. To appear in International Journal for Numerical Methods in Engineering. [9] A. Fuchs. Optimierung von Delaunay-Triangulierungen. Master's thesis, Universitat Stuttgart, 996. [0] A. Fuchs. Almost regular Delaunay-Triangulations. Int. J. for Numerical Methods in Engineering, 40:4595{460, 997. [] A. Fuchs. Parallel Delaunay triangulation. in preparation, University of Stuttgart, 998. [2] NAGD Research Group. Software. Report 96-5, University of Stuttgart, 996. [3] J. Horner. Condition Numbers of FE-Matrices. private communication. [4] K. Hollig and H. Mogerle. G-splines. Computer Aided Geometric Design, 7:97{207, 990. [5] J. Hoschek and D. Lasser. Grundlagen der Geometrischen Datenverarbeitung. B.G. Teubner, Stuttgart, 989. [6] N. A. Golias and T. D. Tsiboukis. An Approach to Rening Three-Dimensional Tetrahedral Meshes Based on Delaunay Transformations. Int. J. for Numerical Methods in Engineering, 37:793{82, 994. [7] D. T. Lee and B. J. Schachter. Two Algorithms for Constructing a Delaunay Triangulation. Int. J. Comput. Inf. Sci. 9, 9:29{242, 980. [8] D. Liu and J. Hoschek. GC Continuity Conditions between Adjacent Rectangular and Triangular Bezier Surface Patches. CAD, 2(4):94{200, May 989. [9] Franco P. Preparata and Michael Shamos. Computational Geometry: An Introduction. Texts and Monographs in Computer Science. Springer-Verlag, 985.
N. Hitschfeld. Blanco Encalada 2120, Santiago, CHILE.
Generalization of modied octrees for geometric modeling N. Hitschfeld Dpto. Ciencias de la Computacion, Univ. de Chile Blanco Encalada 2120, Santiago, CHILE E-mail: nancy@dcc.uchile.cl Abstract. This paper
More informationWe consider the problem of rening quadrilateral and hexahedral element meshes. For
Rening quadrilateral and hexahedral element meshes R. Schneiders RWTH Aachen Lehrstuhl fur Angewandte Mathematik, insb. Informatik Ahornstr. 55, 5056 Aachen, F.R. Germany (robert@feanor.informatik.rwth-aachen.de)
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 information2) For any triangle edge not on the boundary, there is exactly one neighboring
Triangulating Trimmed NURBS Surfaces Chang Shu and Pierre Boulanger Abstract. This paper describes techniques for the piecewise linear approximation of trimmed NURBS surfaces. The problem, called surface
More informationTiling Three-Dimensional Space with Simplices. Shankar Krishnan AT&T Labs - Research
Tiling Three-Dimensional Space with Simplices Shankar Krishnan AT&T Labs - Research What is a Tiling? Partition of an infinite space into pieces having a finite number of distinct shapes usually Euclidean
More information13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY
13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY Lecture 23 Dr. W. Cho Prof. N. M. Patrikalakis Copyright c 2003 Massachusetts Institute of Technology Contents 23 F.E. and B.E. Meshing Algorithms 2
More informationOn a nested refinement of anisotropic tetrahedral grids under Hessian metrics
On a nested refinement of anisotropic tetrahedral grids under Hessian metrics Shangyou Zhang Abstract Anisotropic grids, having drastically different grid sizes in different directions, are efficient and
More informationTilings of the Euclidean plane
Tilings of the Euclidean plane Yan Der, Robin, Cécile January 9, 2017 Abstract This document gives a quick overview of a eld of mathematics which lies in the intersection of geometry and algebra : tilings.
More information8.B. The result of Regiomontanus on tetrahedra
8.B. The result of Regiomontanus on tetrahedra We have already mentioned that Plato s theory that the five regular polyhedra represent the fundamental elements of nature, and in supplement (3.D) to the
More informationNormals of subdivision surfaces and their control polyhedra
Computer Aided Geometric Design 24 (27 112 116 www.elsevier.com/locate/cagd Normals of subdivision surfaces and their control polyhedra I. Ginkel a,j.peters b,,g.umlauf a a University of Kaiserslautern,
More informationWe have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance.
Solid geometry We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance. First, note that everything we have proven for the
More informationAnswer Key: Three-Dimensional Cross Sections
Geometry A Unit Answer Key: Three-Dimensional Cross Sections Name Date Objectives In this lesson, you will: visualize three-dimensional objects from different perspectives be able to create a projection
More informationPebble Sets in Convex Polygons
2 1 Pebble Sets in Convex Polygons Kevin Iga, Randall Maddox June 15, 2005 Abstract Lukács and András posed the problem of showing the existence of a set of n 2 points in the interior of a convex n-gon
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationChapter 8. Voronoi Diagrams. 8.1 Post Oce Problem
Chapter 8 Voronoi Diagrams 8.1 Post Oce Problem Suppose there are n post oces p 1,... p n in a city. Someone who is located at a position q within the city would like to know which post oce is closest
More informationRay shooting from convex ranges
Discrete Applied Mathematics 108 (2001) 259 267 Ray shooting from convex ranges Evangelos Kranakis a, Danny Krizanc b, Anil Maheshwari a;, Jorg-Rudiger Sack a, Jorge Urrutia c a School of Computer Science,
More informationLinear Complexity Hexahedral Mesh Generation
Linear Complexity Hexahedral Mesh Generation David Eppstein Department of Information and Computer Science University of California, Irvine, CA 92717 http://www.ics.uci.edu/ eppstein/ Tech. Report 95-51
More informationCollars and Intestines: Practical Conforming Delaunay Refinement
Collars and Intestines: Practical Conforming Delaunay Refinement Alexander Rand and Noel Walkington Carnegie Mellon University Summary. While several existing Delaunay refinement algorithms allow acute
More informationThree applications of Euler s formula. Chapter 10
Three applications of Euler s formula Chapter 10 A graph is planar if it can be drawn in the plane R without crossing edges (or, equivalently, on the -dimensional sphere S ). We talk of a plane graph if
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More information[8] that this cannot happen on the projective plane (cf. also [2]) and the results of Robertson, Seymour, and Thomas [5] on linkless embeddings of gra
Apex graphs with embeddings of face-width three Bojan Mohar Department of Mathematics University of Ljubljana Jadranska 19, 61111 Ljubljana Slovenia bojan.mohar@uni-lj.si Abstract Aa apex graph is a graph
More informationHyperbolic structures and triangulations
CHAPTER Hyperbolic structures and triangulations In chapter 3, we learned that hyperbolic structures lead to developing maps and holonomy, and that the developing map is a covering map if and only if the
More informationCourse Number: Course Title: Geometry
Course Number: 1206310 Course Title: Geometry RELATED GLOSSARY TERM DEFINITIONS (89) Altitude The perpendicular distance from the top of a geometric figure to its opposite side. Angle Two rays or two line
More informationEXTREME POINTS AND AFFINE EQUIVALENCE
EXTREME POINTS AND AFFINE EQUIVALENCE The purpose of this note is to use the notions of extreme points and affine transformations which are studied in the file affine-convex.pdf to prove that certain standard
More informationimplicit surfaces, approximate implicitization, B-splines, A- patches, surface fitting
24. KONFERENCE O GEOMETRII A POČÍTAČOVÉ GRAFICE ZBYNĚK ŠÍR FITTING OF PIECEWISE POLYNOMIAL IMPLICIT SURFACES Abstrakt In our contribution we discuss the possibility of an efficient fitting of piecewise
More informationImproved Results on Geometric Hitting Set Problems
Improved Results on Geometric Hitting Set Problems Nabil H. Mustafa nabil@lums.edu.pk Saurabh Ray saurabh@cs.uni-sb.de Abstract We consider the problem of computing minimum geometric hitting sets in which,
More informationCAD & Computational Geometry Course plan
Course plan Introduction Segment-Segment intersections Polygon Triangulation Intro to Voronoï Diagrams & Geometric Search Sweeping algorithm for Voronoï Diagrams 1 Voronoi Diagrams Voronoi Diagrams or
More informationVoronoi Diagram. Xiao-Ming Fu
Voronoi Diagram Xiao-Ming Fu Outlines Introduction Post Office Problem Voronoi Diagram Duality: Delaunay triangulation Centroidal Voronoi tessellations (CVT) Definition Applications Algorithms Outlines
More informationSimple Mesh Examples to Illustrate Specific Finite Element Mesh Requirements
Simple Mesh Examples to Illustrate Specific Finite Element Mesh Requirements Peter Fleischmann, Robert Kosik, Bernhard Haindl, and Siegfried Selberherr Institute for Microelectronics, TU Vienna, Gu6hausstrafie
More informationCOMPUTING CONSTRAINED DELAUNAY
COMPUTING CONSTRAINED DELAUNAY TRIANGULATIONS IN THE PLANE By Samuel Peterson, University of Minnesota Undergraduate The Goal The Problem The Algorithms The Implementation Applications Acknowledgments
More informationTopological Issues in Hexahedral Meshing
Topological Issues in Hexahedral Meshing David Eppstein Univ. of California, Irvine Dept. of Information and Computer Science Outline I. What is meshing? Problem statement Types of mesh Quality issues
More informationEnumeration of Full Graphs: Onset of the Asymptotic Region. Department of Mathematics. Massachusetts Institute of Technology. Cambridge, MA 02139
Enumeration of Full Graphs: Onset of the Asymptotic Region L. J. Cowen D. J. Kleitman y F. Lasaga D. E. Sussman Department of Mathematics Massachusetts Institute of Technology Cambridge, MA 02139 Abstract
More informationspline structure and become polynomials on cells without collinear edges. Results of this kind follow from the intrinsic supersmoothness of bivariate
Supersmoothness of bivariate splines and geometry of the underlying partition. T. Sorokina ) Abstract. We show that many spaces of bivariate splines possess additional smoothness (supersmoothness) that
More information1 Automatic Mesh Generation
1 AUTOMATIC MESH GENERATION 1 1 Automatic Mesh Generation 1.1 Mesh Definition Mesh M is a discrete representation of geometric model in terms of its geometry G, topology T, and associated attributes A.
More informationarxiv: v1 [math.co] 17 Jan 2014
Regular matchstick graphs Sascha Kurz Fakultät für Mathematik, Physik und Informatik, Universität Bayreuth, Germany Rom Pinchasi Mathematics Dept., Technion Israel Institute of Technology, Haifa 2000,
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 informationA technique for adding range restrictions to. August 30, Abstract. In a generalized searching problem, a set S of n colored geometric objects
A technique for adding range restrictions to generalized searching problems Prosenjit Gupta Ravi Janardan y Michiel Smid z August 30, 1996 Abstract In a generalized searching problem, a set S of n colored
More informationInterlinked Isohedral Tilings of 3D Space
Interlinked Isohedral Tilings of 3D Space Roman Fuchs Carlo H. Séquin Electrical Engineering and Computer Sciences University of California at Berkeley Technical Report No. UCB/EECS-2008-83 http://www.eecs.berkeley.edu/pubs/techrpts/2008/eecs-2008-83.html
More informationOrthogonal art galleries with holes: a coloring proof of Aggarwal s Theorem
Orthogonal art galleries with holes: a coloring proof of Aggarwal s Theorem Pawe l Żyliński Institute of Mathematics University of Gdańsk, 8095 Gdańsk, Poland pz@math.univ.gda.pl Submitted: Sep 9, 005;
More informationVariational Design and Fairing of Spline Surfaces
EUROGRAPHICS 94 / M. Daehlen and L. Kjelldahl (Guest Editors), Blackwell Publishers Eurographics Association, 1994 Volume 13, (1994), number 3 Variational Design and Fairing of Spline Surfaces Günther
More informationMath 443/543 Graph Theory Notes 11: Graph minors and Kuratowski s Theorem
Math 443/543 Graph Theory Notes 11: Graph minors and Kuratowski s Theorem David Glickenstein November 26, 2008 1 Graph minors Let s revisit some de nitions. Let G = (V; E) be a graph. De nition 1 Removing
More informationAdaptive Tessellation for Trimmed NURBS Surface
Adaptive Tessellation for Trimmed NURBS Surface Ma YingLiang and Terry Hewitt 2 Manchester Visualization Centre, University of Manchester, Manchester, M3 9PL, U.K. may@cs.man.ac.uk 2 W.T.Hewitt@man.ac.uk
More informationNormals of subdivision surfaces and their control polyhedra
Normals of subdivision surfaces and their control polyhedra I. Ginkel, a, J. Peters b, and G. Umlauf a, a University of Kaiserslautern, Germany b University of Florida, Gainesville, FL, USA Abstract For
More informationCONSTRUCTIONS OF QUADRILATERAL MESHES: A COMPARATIVE STUDY
South Bohemia Mathematical Letters Volume 24, (2016), No. 1, 43-48. CONSTRUCTIONS OF QUADRILATERAL MESHES: A COMPARATIVE STUDY PETRA SURYNKOVÁ abstrakt. Polygonal meshes represent important geometric structures
More informationcoding of various parts showing different features, the possibility of rotation or of hiding covering parts of the object's surface to gain an insight
Three-Dimensional Object Reconstruction from Layered Spatial Data Michael Dangl and Robert Sablatnig Vienna University of Technology, Institute of Computer Aided Automation, Pattern Recognition and Image
More informationA new 8-node quadrilateral spline finite element
Journal of Computational and Applied Mathematics 195 (2006) 54 65 www.elsevier.com/locate/cam A new 8-node quadrilateral spline finite element Chong-Jun Li, Ren-Hong Wang Institute of Mathematical Sciences,
More informationMa/CS 6b Class 26: Art Galleries and Politicians
Ma/CS 6b Class 26: Art Galleries and Politicians By Adam Sheffer The Art Gallery Problem Problem. We wish to place security cameras at a gallery, such that they cover it completely. Every camera can cover
More informationQuadratic and cubic b-splines by generalizing higher-order voronoi diagrams
Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24
More informationThe Global Standard for Mobility (GSM) (see, e.g., [6], [4], [5]) yields a
Preprint 0 (2000)?{? 1 Approximation of a direction of N d in bounded coordinates Jean-Christophe Novelli a Gilles Schaeer b Florent Hivert a a Universite Paris 7 { LIAFA 2, place Jussieu - 75251 Paris
More informationComputational Geometry
Motivation Motivation Polygons and visibility Visibility in polygons Triangulation Proof of the Art gallery theorem Two points in a simple polygon can see each other if their connecting line segment is
More information2003/2010 ACOS MATHEMATICS CONTENT CORRELATION GEOMETRY 2003 ACOS 2010 ACOS
CURRENT ALABAMA CONTENT PLACEMENT G.1 Determine the equation of a line parallel or perpendicular to a second line through a given point. G.2 Justify theorems related to pairs of angles, including angles
More informationPlanar Graphs. 1 Graphs and maps. 1.1 Planarity and duality
Planar Graphs In the first half of this book, we consider mostly planar graphs and their geometric representations, mostly in the plane. We start with a survey of basic results on planar graphs. This chapter
More information3 Identify shapes as two-dimensional (lying in a plane, flat ) or three-dimensional ( solid ).
Geometry Kindergarten Identify and describe shapes (squares, circles, triangles, rectangles, hexagons, cubes, cones, cylinders, and spheres). 1 Describe objects in the environment using names of shapes,
More informationLecture 25: Bezier Subdivision. And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10
Lecture 25: Bezier Subdivision And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10 1. Divide and Conquer If we are going to build useful
More informationREINER HORST. License or copyright restrictions may apply to redistribution; see
MATHEMATICS OF COMPUTATION Volume 66, Number 218, April 1997, Pages 691 698 S 0025-5718(97)00809-0 ON GENERALIZED BISECTION OF n SIMPLICES REINER HORST Abstract. A generalized procedure of bisection of
More informationEgemen Tanin, Tahsin M. Kurc, Cevdet Aykanat, Bulent Ozguc. Abstract. Direct Volume Rendering (DVR) is a powerful technique for
Comparison of Two Image-Space Subdivision Algorithms for Direct Volume Rendering on Distributed-Memory Multicomputers Egemen Tanin, Tahsin M. Kurc, Cevdet Aykanat, Bulent Ozguc Dept. of Computer Eng. and
More informationInside hole. Inside wire. Input: and. d y. Inside vertex. Outside loop
Energy-Minimizing Approach to Meshing Curved Wire-Frame Models Atsushi Yamada Kenji Shimada Takayuki Itoh Tokyo Research Laboratory, IBM Japan, Ltd., 1623-14, Shimotsuruma, Yamato-shi, Kanagawa-ken 242,
More informationSimplicial Hyperbolic Surfaces
Simplicial Hyperbolic Surfaces Talk by Ken Bromberg August 21, 2007 1-Lipschitz Surfaces- In this lecture we will discuss geometrically meaningful ways of mapping a surface S into a hyperbolic manifold
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 informationCLASSIFICATION OF SURFACES
CLASSIFICATION OF SURFACES JUSTIN HUANG Abstract. We will classify compact, connected surfaces into three classes: the sphere, the connected sum of tori, and the connected sum of projective planes. Contents
More informationCSG obj. oper3. obj1 obj2 obj3. obj5. obj4
Solid Modeling Solid: Boundary + Interior Volume occupied by geometry Solid representation schemes Constructive Solid Geometry (CSG) Boundary representations (B-reps) Space-partition representations Operations
More informationcontrol polytope. These points are manipulated by a descent method to compute a candidate global minimizer. The second method is described in Section
Some Heuristics and Test Problems for Nonconvex Quadratic Programming over a Simplex Ivo Nowak September 3, 1998 Keywords:global optimization, nonconvex quadratic programming, heuristics, Bezier methods,
More informationSurface Reconstruction. Gianpaolo Palma
Surface Reconstruction Gianpaolo Palma Surface reconstruction Input Point cloud With or without normals Examples: multi-view stereo, union of range scan vertices Range scans Each scan is a triangular 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 information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More informationMiddle School Math Course 3
Middle School Math Course 3 Correlation of the ALEKS course Middle School Math Course 3 to the Texas Essential Knowledge and Skills (TEKS) for Mathematics Grade 8 (2012) (1) Mathematical process standards.
More informationACTUALLY DOING IT : an Introduction to Polyhedral Computation
ACTUALLY DOING IT : an Introduction to Polyhedral Computation Jesús A. De Loera Department of Mathematics Univ. of California, Davis http://www.math.ucdavis.edu/ deloera/ 1 What is a Convex Polytope? 2
More informationHowever, m pq is just an approximation of M pq. As it was pointed out by Lin [2], more precise approximation can be obtained by exact integration of t
FAST CALCULATION OF GEOMETRIC MOMENTS OF BINARY IMAGES Jan Flusser Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodarenskou vez 4, 82 08 Prague 8, Czech
More informationNumerical representation of the quality measures of triangles and. triangular meshes
Numerical representation of the quality measures of triangles and triangular meshes J. Sarrate, J. Palau and A. Huerta Departament de Matemàtica Aplicada III, E.T.S. de Ingenieros de Caminos, Canales y
More informationGeometry Curriculum Map
Quadrilaterals 7.1 Interior Angle Sum Theorem 7.2 Exterior Angle Sum Theorem 7.3 Using Interior and Exterior Angles to Solve Problems Define the Angle Sum Theorem. Illustrate interior angles with the Angle
More informationOn the Relationships between Zero Forcing Numbers and Certain Graph Coverings
On the Relationships between Zero Forcing Numbers and Certain Graph Coverings Fatemeh Alinaghipour Taklimi, Shaun Fallat 1,, Karen Meagher 2 Department of Mathematics and Statistics, University of Regina,
More informationDefinitions. Topology/Geometry of Geodesics. Joseph D. Clinton. SNEC June Magnus J. Wenninger
Topology/Geometry of Geodesics Joseph D. Clinton SNEC-04 28-29 June 2003 Magnus J. Wenninger Introduction Definitions Topology Goldberg s polyhedra Classes of Geodesic polyhedra Triangular tessellations
More informationLecture 2 - Introduction to Polytopes
Lecture 2 - Introduction to Polytopes Optimization and Approximation - ENS M1 Nicolas Bousquet 1 Reminder of Linear Algebra definitions Let x 1,..., x m be points in R n and λ 1,..., λ m be real numbers.
More informationγ 2 γ 3 γ 1 R 2 (b) a bounded Yin set (a) an unbounded Yin set
γ 1 γ 3 γ γ 3 γ γ 1 R (a) an unbounded Yin set (b) a bounded Yin set Fig..1: Jordan curve representation of a connected Yin set M R. A shaded region represents M and the dashed curves its boundary M that
More informationStudents are not expected to work formally with properties of dilations until high school.
Domain: Geometry (G) Cluster: Understand congruence and similarity using physical models, transparencies, or geometry software. Standard: 8.G.1. Verify experimentally the properties of rotations, reflections,
More informationMath 462: Review questions
Math 462: Review questions Paul Hacking 4/22/10 (1) What is the angle between two interior diagonals of a cube joining opposite vertices? [Hint: It is probably quickest to use a description of the cube
More informationLocalization in Graphs. Richardson, TX Azriel Rosenfeld. Center for Automation Research. College Park, MD
CAR-TR-728 CS-TR-3326 UMIACS-TR-94-92 Samir Khuller Department of Computer Science Institute for Advanced Computer Studies University of Maryland College Park, MD 20742-3255 Localization in Graphs Azriel
More informationMatching and Planarity
Matching and Planarity Po-Shen Loh June 010 1 Warm-up 1. (Bondy 1.5.9.) There are n points in the plane such that every pair of points has distance 1. Show that there are at most n (unordered) pairs of
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 informationApproximation of 3D-Parametric Functions by Bicubic B-spline Functions
International Journal of Mathematical Modelling & Computations Vol. 02, No. 03, 2012, 211-220 Approximation of 3D-Parametric Functions by Bicubic B-spline Functions M. Amirfakhrian a, a Department of Mathematics,
More informationCombinatorial Tilings of the Sphere by Pentagons
Combinatorial Tilings of the Sphere by Pentagons Min Yan Department of Mathematics Hong Kong University of Science and Technology Kowloon, Hong Kong mamyan@ust.hk Submitted: Nov 16, 2012; Accepted: Mar
More informationNon-extendible finite polycycles
Izvestiya: Mathematics 70:3 1 18 Izvestiya RAN : Ser. Mat. 70:3 3 22 c 2006 RAS(DoM) and LMS DOI 10.1070/IM2006v170n01ABEH002301 Non-extendible finite polycycles M. Deza, S. V. Shpectorov, M. I. Shtogrin
More informationAcute Triangulations of Polygons
Europ. J. Combinatorics (2002) 23, 45 55 doi:10.1006/eujc.2001.0531 Available online at http://www.idealibrary.com on Acute Triangulations of Polygons H. MAEHARA We prove that every n-gon can be triangulated
More informationGeometry Geometry Grade Grade Grade
Grade Grade Grade 6.G.1 Find the area of right triangles, special quadrilaterals, and polygons by composing into rectangles or decomposing into triangles and other shapes; apply these techniques in the
More informationUNSTRUCTURED GRIDS ON NURBS SURFACES. The surface grid can be generated either in a parameter. surfaces. Generating grids in a parameter space is
UNSTRUCTURED GRIDS ON NURBS SURFACES Jamshid Samareh-Abolhassani 1 Abstract A simple and ecient computational method is presented for unstructured surface grid generation. This method is built upon an
More informationLecture notes: Object modeling
Lecture notes: Object modeling One of the classic problems in computer vision is to construct a model of an object from an image of the object. An object model has the following general principles: Compact
More informationThe Art Gallery Problem
The Art Gallery Problem Connor Bohlken Nebraska Wesleyan University cbohlken@nebrwesleyan.edu May 11, 2016 Abstract Given a polygonal art gallery of n vertices, how many security guards are necessary to
More informationRegularity Analysis of Non Uniform Data
Regularity Analysis of Non Uniform Data Christine Potier and Christine Vercken Abstract. A particular class of wavelet, derivatives of B-splines, leads to fast and ecient algorithms for contours detection
More informationRigidity of ball-polyhedra via truncated Voronoi and Delaunay complexes
!000111! NNNiiinnnttthhh IIInnnttteeerrrnnnaaatttiiiooonnnaaalll SSSyyymmmpppooosssiiiuuummm ooonnn VVVooorrrooonnnoooiii DDDiiiaaagggrrraaammmsss iiinnn SSSccciiieeennnccceee aaannnddd EEEnnngggiiinnneeeeeerrriiinnnggg
More informationPoint Cloud Simplification and Processing for Path-Planning. Master s Thesis in Engineering Mathematics and Computational Science DAVID ERIKSSON
Point Cloud Simplification and Processing for Path-Planning Master s Thesis in Engineering Mathematics and Computational Science DAVID ERIKSSON Department of Mathematical Sciences Chalmers University of
More informationThe Full Survey on The Euclidean Steiner Tree Problem
The Full Survey on The Euclidean Steiner Tree Problem Shikun Liu Abstract The Steiner Tree Problem is a famous and long-studied problem in combinatorial optimization. However, the best heuristics algorithm
More informationMISCELLANEOUS SHAPES
MISCELLANEOUS SHAPES 4.1. INTRODUCTION Five generic shapes of polygons have been usefully distinguished in the literature: convex, orthogonal, star, spiral, and monotone. 1 Convex polygons obviously do
More informationA Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
Scientific Papers, University of Latvia, 2010. Vol. 756 Computer Science and Information Technologies 207 220 P. A Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
More informationObtaining the H and T Honeycomb from a Cross-Section of the 16-cell Honeycomb
Bridges 2017 Conference Proceedings Obtaining the H and T Honeycomb from a Cross-Section of the 16-cell Honeycomb Hideki Tsuiki Graduate School of Human and Environmental Studies, Kyoto University Yoshida-Nihonmatsu,
More informationLATEST TRENDS on APPLIED MATHEMATICS, SIMULATION, MODELLING
3D surface reconstruction of objects by using stereoscopic viewing Baki Koyuncu, Kurtuluş Küllü bkoyuncu@ankara.edu.tr kkullu@eng.ankara.edu.tr Computer Engineering Department, Ankara University, Ankara,
More informationMonotone Paths in Geometric Triangulations
Monotone Paths in Geometric Triangulations Adrian Dumitrescu Ritankar Mandal Csaba D. Tóth November 19, 2017 Abstract (I) We prove that the (maximum) number of monotone paths in a geometric triangulation
More informationAdvanced Modeling 2. Katja Bühler, Andrej Varchola, Eduard Gröller. March 24, x(t) z(t)
Advanced Modeling 2 Katja Bühler, Andrej Varchola, Eduard Gröller March 24, 2014 1 Parametric Representations A parametric curve in E 3 is given by x(t) c : c(t) = y(t) ; t I = [a, b] R z(t) where x(t),
More informationConstrained Delaunay Triangulations (CDT)
Constrained Delaunay Triangulations (CDT) Recall: Definition 1 (Constrained triangulation). A triangulation with prespecified edges or breaklines between nodes. Motivation: Geological faults in oil recovery
More informationGeneralizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that g
Generalizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that generalizes the four-directional box spline of class
More informationGeometry I Can Statements I can describe the undefined terms: point, line, and distance along a line in a plane I can describe the undefined terms:
Geometry I Can Statements I can describe the undefined terms: point, line, and distance along a line in a plane I can describe the undefined terms: point, line, and distance along a line in a plane I can
More information