852 QIN Kaihuai, CHANG Zhengyi et al Vol7 are generated by the recursive subdivision Each subdivision step will bring new vertices, new edges, as well
|
|
- Juliana Patterson
- 6 years ago
- Views:
Transcription
1 Vol7 No6 J Comput Sci & Technol Nov 22 Physics-ased Loop Surface Modeling QIN Kaihuai ( ΠΞ), CHANG Zhengyi ( ffω), WANG Huawei (ΦΛΨ) and LI Denggao (± ) Department of Computer Science and Technology, Tsinghua University, eijing 84, PR China Abstract qkh-dcstsinghuaeducn Received July 5, 2; revised September 2, 2 Strongly inspired by the research on physics-based dynamic models for surfaces, we propose a new method for precisely evaluating the dynamic parameters (mass, damping and stiffness matrices, and dynamic forces) for Loop surfaces without recursive subdivision regardless of regular or irregular faces It is shown that the thin-plate-energy of Loop surfaces can be evaluated precisely and efficiently, even though there are extraordinary points in the initial meshes, unlike the previous dynamic Loop surface scheme Hence, the new method presented for Loop surfaces is much more efficient than the previous schemes Keywords subdivision, Loop surface, physics-based modeling Introduction Free-form curves and surfaces are very important and common in CAGD and CAD/CAM The Non-Uniform Rational -Spline (NURS) is a very popular tool, but there are too many parameters in NURS curves and surfaces Generally, a surface of arbitrary topology cannot be represented by a single NURS surface It is difficult for users to know how to set the knot spacings, how to manipulate the control vertices and how to modify the weights In addition, NURS surfaces are defined on a rectangular parametric domain This has hampered their widespread uses in modeling shapes On the other hand, a subdivision surface can represent any complex surface of arbitrary topology Physics-based subdivision surfaces allow users to physically modify shapes of surfaces at desired locations via forces applied Furthermore, physics-based modeling gives the users an intuitive and natural feeling for geometric modeling with clay dough Thus, in recent years, attention has been paid to physics-based dynamic subdivision surfaces [ 3] In 997, Qin, Mandal and Vemuri introduced the physical quantity" into dynamic Catmull-Clark surfaces and successfully applied it in visualization of medical data [2] ut in case of extraordinary points, the faces that contain an extraordinary point were subdivided recursively to get the mass, damping and stiffness matrices Hence not only does the number of equations of the numerical integrals increase quickly, but also the result is approximate This is because their approach used for computing the thin-plate-energy of Catmull-Clark surfaces diverges sometimes [2;4] In 2, Mandal, Qin and Vemuri [] gave an approximate solution to the physical matrices of Loop surfaces They failed to get an exactly analytic solution of the dynamic matrices like mass, damping and stiffness ones, and only got the approximations of such physical matrices by a spring-mass system Qin, Wang and Li, et al presented efficient methods for precisely computing the dynamic parameters of physics-based catmull-clark surfaces [3] In this paper, a new method for physics-based Loop surfaces is presented, the thin-plate-energy of Loop surfaces is evaluated precisely, and the physical matrices are evaluated exactly and analytically 2 Loop Subdivision Surfaces In 987, Loop generalized the box splines to irregular meshes [5] The surfaces defined by the control meshes that have extraordinary points can be processed efficiently by Loop's scheme Patches Supported by the National Natural Science Foundation of China (Grant Nos and ) and the Research Fund for the Doctoral Program of Higher Education from China Ministry of Education (Grant No2348)
2 852 QIN Kaihuai, CHANG Zhengyi et al Vol7 are generated by the recursive subdivision Each subdivision step will bring new vertices, new edges, as well as new triangles Loop's subdivision rules are listed as follows [5;6] : ffl New vertex points Suppose that for each non-boundary vertex V whose valence is N A new vertex is introduced corresponding to the old one, and the positions of the new vertex is given by (ff(n)v + V + V V N )=(ff(n)+n); ff(n) = (N( fi(n)))=fi(n); fi(n) = 5=8 (3 + 2 cos(2ß=n)) 2 =64: V is the position of vertex V, and V ; V 2 ;:::;V N are the vertex positions of N vertices connected to the vertex V ffl New edge points For each non-boundary edge E, a new edge point is introduced by (3(V + V 2 )+V F + V F2 )=8, V and V 2 are the end points of the edge E, V F and V F2 are the third points of the two triangles sharing the edge E, respectively ffl New edges New edges are formed by connecting each new vertex point to the new edge points corresponding to the old edges incident on the old vertices, and connecting each new edge point to the new edge points of other edges in two faces which share the original edge ffl New triangles New triangles are defined as faces enclosed by the new edges 3 Physics-ased Deformable Loop Surface Model In order to make the derivation below as clear and compact as possible we adopt the following notational conventions: ffl All the matrices are denoted by bold capital letters, eg, M; D; K ffl The element of the i-th row, j-th column of matrix M is denoted by (M) ij ffl The vectors are denoted by bold lowercase letters, eg, c 3 Parametric Expression of Loop Surfaces Each patch of a Loop subdivision surface [5;6] can be defined by 2 vertices Lai [7] got the basis functions of the 2 vertices by transforming the box splines into triangular ézier patches All the basis functions are shown in Appendix A In Fig, the shaded Loop patch has the formula [6] as follows: s(v; w) = C T s b(v; w); (v; w) 2 Ω; Ω = f(v; w)jv 2 [; ] and w 2 [; v]g, C T s = (c ;:::;c 2 ) is a 3 2 matrix whose elements are the coordinates of the 2 vertices, b(v; w) is a vector composed of the 2 basic functions (see Appendix A) For an irregular patch (with a vertex whose valence is not equal to 6) as shown in Fig2 Let K = N + 6 and M = K + 6, N is the valence of the extraordinary point, K is the number of vertices that define the shaded face The coordinates of the surface at any parameters (v; w) can be computed by eigenanalysis as follows [6] : s(v; w) = C T s V T Λ n (P k AV ) T b(t k;n (v; w)) () Λ is a diagonal matrix composed of the eigenvalues, V is the eigenvector matrix corresponding to the eigenvalues (see Appendix C), A is the extended subdivision matrix, P k (k = ; 2; 3) are 2 M picking" matrices, and the transformations t k;n (k = ; 2; 3) map Ω n k to the unit triangles (see Fig3): t ;n (v; w) = (2 n v ; 2 n w); t ;n 2 Ω n ; t 2;n (v; w) = ( 2 n v; 2 n w); t 2;n 2 Ω n 2 ; t 3;n (v; w) = (2 n v; 2 n w ); t 3;n 2 Ω n 3 : The definitions of A, A and P k can be found in Appendix
3 No6 Physics-ased Loop Surface Modeling 853 Fig A Loop patch Fig2 An irregular triangular patch defined by K = N +6= 3 control vertices with an extraordinary vertex of valence 7 () can be rewritten as follows: s T (v; w) = J Λ C s ; J Λ = b T (t k;n (v; w))x k Λ n V ; X k = P k AV : If the mesh has m vertices and d faces, then C T s = J i C s = J = J Λ R i C = JC ; Fig3 The parametric domain of Loop subdivision J Λ R i : (2) = (c ;:::;c m ) is a matrix whose elements are the coordinates of the m vertices R i is a k m picking" matrices (each row has only one " element, the others are ) 32 Dynamics Equation of Loop Subdivision Surfaces Consider the vertex positions C as time-varying variables, then the dynamics equation can be expressed as follows: M C + D _C + KC = F p (3) M; D and K are the mass, damping and stiffness matrices of the physical model, respectively, and F p is the generalized force vector Let μ be the density of mass, fl be the coefficient of damping, then the mass and damping matrices are expressed respectively as follows: 8 >< >: M = D = μj T Jdw flj T Jdw The stiffness matrix can be obtained by computing the thin-plate-energy" of each face: K = (ff J T v J v + ff 22 J T wj w + fi J T vvj vv + fi 2 J T vwj vw + fi 22 J T wwj ww )dw ff ii and fi ij are the tension and the rigidity functions Generally, μ; fl; ff ii ;fi ij are the functions of v and w Let f(v; w; t) be the force imposed at the point (v; w), then the force vector F p can be expressed as follows: F p = J T f T (v; w; t)dw: In general, the force f(v; w; t) is equal to a constant on a small face subdivided (4)
4 854 QIN Kaihuai, CHANG Zhengyi et al Vol7 33 Evaluation of Loop Surfaces Suppose the triangular faces subdivided are small, such that μ; fl; ff ii and fi ij can be regarded as constants Combining (2) and (4) yields M = μj T Jdw = μ(j Λ R i ) T (J Λ R i )dw = μ i is the mass density of the i-th face Similarly, we can obtain the damping matrix D: D = fl i R T i J ΛT J Λ dw Ri fl i is the damping coefficient of the i-th face y evaluating μ i R T i J ΛT J Λ dw Ri J ΛT J Λ dw, one can get the matrices M and D Note that, the integral can be calculated on the sub-domains Ω i, Ω i 2 and Ω i 3 (see Fig3) Hence we have the following expressions: Ω Ω 2 J ΛT J Λ dw = J ΛT J Λ dw = Ω V T Λ X T b(2v ; 2w)b T (2v ; 2w)X Λ V dw Z = V T Λ X T = V T Λ X T 4 Ω 2 = V T Λ X T = V T Λ X T =2 Z b(2v ; 2w)b T (2v ; 2w)dw X Λ V X Λ V ; V T Λ X T b(4v ; 4w)b T (4v ; 4w)X Λ V dw Z =2 =4 4 2 Z Thus, from the above we can deduce that Ω n J ΛT J Λ dw = 4 n V T Λ n X T Similarly, we can obtain the following formulae: Ω n 2 Ω n 3 J ΛT J Λ dw = 4 V T n Λ n X T 2 J ΛT J Λ dw = 4 n V T Λ n X T 3 Z =2 v Z Z Z Hence, we can obtain the following unified expression: Ω n k J ΛT J Λ dw = 4 V T n Λ n X T k Let Z k = X T k Z Z b(4v ; 4w)b T (4v ; 4w)dw X Λ V X Λ V : X Λ n V : X 2 Λ n V ; X 3 Λ n V : Xk Λ n V ; k = ; 2; 3: Xk and Q P k = n= 4 n Λn Z k Λ n, the elements of the i-th row, j-th column of the matrices Q k and Z k are denoted by (Q k ) ij and (Z k ) ij, respectively, then we have J ΛT J Λ dw = 3X X k= n= Ω n k J ΛT J Λ dw = 3X k= V T Q k V
5 No6 Physics-ased Loop Surface Modeling 855 It is clear that (Q k ) ij = 4 X n= i j n (Zk ) ij (Z k ) ij : 4 4 i j J ΛT J Λ dw can be evaluated precisely and efficiently Since the resulting matrices are only dependent on the valence N of the extraordinary point, they can be evaluated during preprocessing Similarly, we can precisely and efficiently evaluate K and F p by computing in aance J ΛT v J Λ v dw, J ΛT w J Λ w dw, J ΛT vv J Λ vv dw, J ΛT vwj Λ vw dw, J ΛT wwj Λ ww dw and J ΛT dw, respectively Thus, we can precisely and quickly evaluate the dynamic parameters M, D, K, F p for Loop surfaces, and this is a novel contribution of this paper 34 Solving the Differential Equations We can use the difference method to solve (3) If t is small enough, then C (t + t) = C (t + t) 2C (t)+c (t t) ; _C t 2 (t + t) = C (t + t) C (t t) 2 t Then, the differential equation can be written as (2M + D t + 2 t 2 K)C (t + t) = 2 t 2 F p (t + t)+(d t 2M)C (t t)+4mc (t): The above linear system of equations can be solved quickly by the LU decomposition of the coefficient matrix Since F p is timely-varying, F p should be re-evaluated at each iteration step 35 Application Examples There are two examples dynamically generated by the presented method as follws: Fig4 A vase deformed (a) Initial shape (b) Deformed shape after seconds (c) Deformed shape after 5 seconds Fig5 Dynamically modeling a kettle (a) Initial shape (b) Deformed shape after 5 seconds (c) Deformed shape after seconds
6 856 QIN Kaihuai, CHANG Zhengyi et al Vol7 Example Fig4 shows a vase The initial control vertices, the Loop surface and 6 initial forces on the boundary of the surface are shown in Fig4(a); the deformed shapes are shown in Fig4(b) and (c), respectively Example 2 Fig5 shows a kettle" generated dynamically Fig5(a) shows the initial shape, control vertices and force on the extraordinary point whose valence is 5; Fig5(b) and (c) show the deformed surfaces after 5 and seconds, respectively 4 Computational Complexity First, let us recall the previous physics-based Loop surfaces [] proposed by Mandal, Qin and Vemuri They used 4 j triangular faces obtained after j subdivision steps to approximate the smooth patch in the limit surface corresponding to a face in the initial mesh In their method, the mass matrix M is calculated using the following formula (the same way for calculating D, K and f p ): M abc = kx μ(v j i )fj abc (vj i )gt f j abc (vj i )g; k is the number of vertices of the 4 j faces, v j i is the i-th vertex of the 4 j faces in the j-th level mesh, j is the collection of abc basic functions at the vertices v, which is of dimension 3r and is abc the concatenation of the (x; y; z) positions for all r vertices that are within the 2-neighourhood of the vertices a, b, c in the initial mesh as show in Fig6 Note that the cost for calculating M is dependent on the subdivision level j, all the matrices (M; D; K and f p ) must be calculated Fig6 Points a; b; c and their 2-neighbourhoods again in the next subdivision step, and the results are approximative At level j of the subdivision, the number of vertices of the 4 j faces k = 3j+ + 2 Suppose that 2 there are m triangular faces in the initial mesh, then the computational complexity of M is O(m3 j ), increasing exponentially However, the complexity of the method proposed in the paper for calculating M is O(m), depending only on the number of faces in the initial mesh In addition, all the matrices (M; D; K and f p ) can be calculated exactly and saved in aance So the computational complexity of the new method is a constant, much better than the previous method 5 Conclusions and Future Work In this paper, based on the dynamic model, we can obtain the precise solution of physics-based Loop surfaces We have proposed a unified method to evaluate the physical parameters of the dynamic model regardless of a regular or an irregular face The matrices M; D; K and F p can be obtained in much less time than the previous schemes, because they can be calculated using the analytic formulae in our method without recursive subdivision and numerical integrals [] So our method is much more efficient The thin-plate-energy of Loop surfaces can be evaluated precisely and analytically, even though there are extraordinary points in the control meshes The method proposed can be used easily for modeling any complex shapes of surfaces of arbitrary topology, and the shapes can be modified conveniently It is easy to be implemented It can also be used in medical visualization, which is our future work References [] Mandal C, Qin H, Vemuri C A novel FEM-based dynamic framework for subdivision surfaces Computer-Aided Design, 2, 32: [2] Qin H, Mandal C, Vemuri C Dynamic Catmull-Clark subdivision surfaces IEEE Transactions on Visualization and Computer Graphics, 998, 4(3): [3] Qin K, Wang H, Li D, Kikinis R, Halle M Physics-based subdivision surface modeling for medical imaging and simulation In Proc of MIAR'2, IEEE Computer Society Press, Hong Kong, 2, pp7 24
7 No6 Physics-ased Loop Surface Modeling 857 [4] Halstead M, Kass M, DeRose T Efficient, fair interpolation using Catmull-Clark surfaces Computer Graphics (Proceedings of SIGGRAPH'93), 993, pp35 44 [5] Loop C Smooth subdivision surfaces based on triangles [Thesis] Department of Mathematics, University of Utah, 987 [6] Stam J Evaluation of Loop subdivision surface In Proceedings CD of SIGGRAPH'98, 998 [7] Lai M J Fortrain subroutines for -nets of box splines on three- and four-directional meshes Numerical Algorithms, 992, 2: QIN Kaihuai is a professor in Department of Computer Science and Technology, Tsinghua University His research interests include graphics, computer data visualization, image processing, computer aided geometric design, physics-based geometric modeling, virtual reality and CAD/CAM CHANG Zhengyi is a software engineer He got his bachelor's degree in Hunan University in 988 He is now a graduate student in the Department of Computer Science and Technology, Tsinghua University His research interests include computer graphics and computer aided geometric modeling WANG Huawei is a PhD candidate in the Department of Computer Science and Technology, Tsinghua University His research interests include computer graphics, computer aided geometric design and physicsbased geometric modeling LI Denggao is a PhD candidate in the Department of Computer Science and Technology, Tsinghua University His research interests include computer graphics, computer aided geometric design and physicsbased computer aid geometric modeling Appendix A The vector b(v; w) contains 2 bi-quadratic ezier basis functions as follows [7] : b T (v; w) = 2 (u4 +2u 3 v; u 4 +2u 3 w; u 4 +2u 3 w +6u 3 v +6u 2 vw +2u 2 v 2 +6uv 2 w +6uv 3 +2v 3 w + v 4 ; 6u 4 u = v w +24u 3 w +24u 2 w 2 +8uw 3 + w 4 +24u 3 v +6u 2 vw +36uvw 2 +6vw 3 +24u 2 v 2 +36uv 2 w +2v 2 w 2 +8uv 3 +6v 3 w + v 4 ;u 4 +6u 3 w +2u 2 w 2 +6uw 3 + w 4 +2u 3 v +6u 2 vw +6uvw 2 +2vw 3 ; 2uv 3 + v 4 ;u 4 +6u 3 w +2u 2 w 2 +6uw 3 + w 4 +8u 3 v +36u 2 vw +36uvw 2 +8vw 3 +24u 2 v 2 +6uv 2 w +24v 2 w 2 +24uv 3 +24v 3 w +6v 4 ;u 4 +8u 3 w +24u 2 w 2 +24uw 3 +6w 4 +6u 3 v +36u 2 vw +6uvw 2 +24vw 3 +2u 2 v 2 +36uv 2 w +24v 2 w 2 +6uv 3 +8v 3 w + v 4 ; 2uw 3 + w 4 ; 2v 3 w + v 4 ; 2uw 3 + w 4 +6uvw 2 +6vw 3 +6uv 2 w +2v 2 w 2 +2uv 3 +6v 3 w + v 4 ;w 4 +2vw 3 ); Appendix The K K subdivision matrix A and the M K extended subdivision matrix A respectively have the following block forms: ψ S S A = ; A = S S 2 the matrix S is an (N +) (N +)matrix as follows [6] : S =! S S 2 S 2 S 22 a N b N b N b N b N b N b N b N b N b N c c d d c d c d c d c d c d c d c d c d c d d c C A ;
8 858 QIN Kaihuai, CHANG Zhengyi et al Vol7 a N = ff N ; b N = ff N =N; c = 3=8; d = =8; ff N = 5 (3 + 2 cos(2ß=n))2 ; S 2 = 6 2 C A ; S = C A ; S 2 = C A ; S22 = C A : 3 3 P k (k = ; 2; 3) is a 2 M picking" matrix, there is only one " in each line of P k, and other elements equal P k (k = ; 2; 3) are defined by introducing the following 3 vectors: q = (3; ;N +4; 2;N +;N +9;N +3;N +2;N +5;N +8;N +7;N + ); q 2 = (N +;N +7;N +5;N +2;N +3;N +6;N +; 2;N +4;N;; 3); q 3 = (;N;2;N +;N +6;N +3;N +2;N +5;N +2;N +7;N +;N + ): If the i-th element of q k equals j, then (P k ) ij = and other elements of the i-th row of P k all equal Appendix C The eigenstructure of matrix A is denoted by (V ; Λ), Λ is the diagonal matrix containing the eigenvalues of A, and V is an invertible matrix whose columns are the corresponding eigenvectors, that is AV = V Λ, [6] Λ = ±; ;U ; W can be evaluated as follows: ± ; V = U : U W ± = ρ diag(;μ2;μ 3;μ 3;:::;μ (N )=2;μ (N )=2); diag(;μ 2;μ 3;μ 3;:::;μ N=2 ;μ N=2 ; =8); if N is odd; if N is even μ = ;μ 2 = 5 ffn ; μ3 = f();:::;μn+ = f(n ); 8 and f(k) = ßk cos( 8 N ); = diag 8 ; 8 ; 8 ; 6 ; ; 6 W = C A ; U = C(k) = cos(2ßk=n); S(k) = sin(2ßk=n); 8 ffn 3 C() S() C(2) C(2) S(2) C(4) C(N ) S(N ) C(2(N )) C A ; U is computed by solving the linear system of equation S U + S 2U = U ±:
Evaluation of Loop Subdivision Surfaces
Evaluation of Loop Subdivision Surfaces Jos Stam Alias wavefront, Inc. 8 Third Ave, 8th Floor, Seattle, WA 980, U.S.A. jstam@aw.sgi.com Abstract This paper describes a technique to evaluate Loop subdivision
More informationG 2 Interpolation for Polar Surfaces
1 G 2 Interpolation for Polar Surfaces Jianzhong Wang 1, Fuhua Cheng 2,3 1 University of Kentucky, jwangf@uky.edu 2 University of Kentucky, cheng@cs.uky.edu 3 National Tsinhua University ABSTRACT In this
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 informationSubdivision based Interpolation with Shape Control
Subdivision based Interpolation with Shape Control Fengtao Fan University of Kentucky Deparment of Computer Science Lexington, KY 40506, USA ffan2@uky.edu Fuhua (Frank) Cheng University of Kentucky Deparment
More informationNon-Uniform Recursive Doo-Sabin Surfaces (NURDSes)
Non-Uniform Recursive Doo-Sabin Surfaces Zhangjin Huang 1 Guoping Wang 2 1 University of Science and Technology of China 2 Peking University, China SIAM Conference on Geometric and Physical Modeling Doo-Sabin
More informationModified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation
Modified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation Karl-Heinz Brakhage RWTH Aachen, 55 Aachen, Deutschland, Email: brakhage@igpm.rwth-aachen.de Abstract In this paper
More informationUsing Semi-Regular 4 8 Meshes for Subdivision Surfaces
Using Semi-Regular 8 Meshes for Subdivision Surfaces Luiz Velho IMPA Instituto de Matemática Pura e Aplicada Abstract. Semi-regular 8 meshes are refinable triangulated quadrangulations. They provide a
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationSubdivision Curves and Surfaces
Subdivision Surfaces or How to Generate a Smooth Mesh?? Subdivision Curves and Surfaces Subdivision given polyline(2d)/mesh(3d) recursively modify & add vertices to achieve smooth curve/surface Each iteration
More informationExact Evaluation Of Catmull-Clark Subdivision Surfaces At Arbitrary Parameter Values
Exact Evaluation Of Catmull-Clark Subdivision Surfaces At Arbitrary Parameter Values Jos Stam Alias wavefront Inc Abstract In this paper we disprove the belief widespread within the computer graphics community
More informationNon-Uniform Recursive Doo-Sabin Surfaces
Non-Uniform Recursive Doo-Sabin Surfaces Zhangjin Huang a,b,c,, Guoping Wang d,e a School of Computer Science and Technology, University of Science and Technology of China, PR China b Key Laboratory of
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 informationPolar Embedded Catmull-Clark Subdivision Surface
Polar Embedded Catmull-Clark Subdivision Surface Anonymous submission Abstract In this paper, a new subdivision scheme with Polar embedded Catmull-Clark mesh structure is presented. In this new subdivision
More informationu 0+u 2 new boundary vertex
Combined Subdivision Schemes for the design of surfaces satisfying boundary conditions Adi Levin School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel. Email:fadilev@math.tau.ac.ilg
More informationHoneycomb Subdivision
Honeycomb Subdivision Ergun Akleman and Vinod Srinivasan Visualization Sciences Program, Texas A&M University Abstract In this paper, we introduce a new subdivision scheme which we call honeycomb subdivision.
More informationSubdivision Depth Computation for Extra-Ordinary Catmull-Clark Subdivision Surface Patches
Subdivision Depth Computation for Extra-Ordinary Catmull-Clark Subdivision Surface Patches Fuhua Frank Cheng,GangChen, and Jun-Hai Yong University of Kentucky, Lexington, KY, USA Tsinghua University, Beijing,
More informationLocal Modification of Subdivision Surfaces Based on Curved Mesh
Local Modification of Subdivision Surfaces Based on Curved Mesh Yoshimasa Tokuyama Tokyo Polytechnic University tokuyama@image.t-kougei.ac.jp Kouichi Konno Iwate University konno@cis.iwate-u.ac.jp Junji
More informationSubdivision. Outline. Key Questions. Subdivision Surfaces. Advanced Computer Graphics (Spring 2013) Video: Geri s Game (outside link)
Advanced Computer Graphics (Spring 03) CS 83, Lecture 7: Subdivision Ravi Ramamoorthi http://inst.eecs.berkeley.edu/~cs83/sp3 Slides courtesy of Szymon Rusinkiewicz, James O Brien with material from Denis
More informationFEM-Based Dynamic Subdivision Splines
FEM-Based Dynamic Subdivision Splines Hong Qin Department of Computer Science State University of New York at Stony Brook Stony Brook, NY 11794 4400 Tel: (631) 632 8450; Fax: (631) 632 8334 Email: qin@cs.sunysb.edu
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 informationUNIVERSITY OF CALGARY. Subdivision Surfaces. Advanced Geometric Modeling Faramarz Samavati
Subdivision Surfaces Surfaces Having arbitrary Topologies Tensor Product Surfaces Non Tensor Surfaces We can t find u-curves and v-curves in general surfaces General Subdivision Coarse mesh Subdivision
More informationSurface Quality Assessment of Subdivision Surfaces on Programmable Graphics Hardware
Sur Quality Assessment of Subdivision Surs on Programmable Graphics Hardware Yusuke Yasui Takashi Kanai Keio University SFC Faculty of Environmental Information 53 Endo, Fujisawa, Kanagawa, 5-850, JAPAN.
More informationUntil now we have worked with flat entities such as lines and flat polygons. Fit well with graphics hardware Mathematically simple
Curves and surfaces Escaping Flatland Until now we have worked with flat entities such as lines and flat polygons Fit well with graphics hardware Mathematically simple But the world is not composed of
More informationA Continuous 3-D Medial Shape Model with Branching
A Continuous 3-D Medial Shape Model with Branching Timothy B. Terriberry Guido Gerig Outline Introduction The Generic 3-D Medial Axis Review of Subdivision Surfaces Boundary Reconstruction Edge Curves
More informationSubdivision Curves and Surfaces: An Introduction
Subdivision Curves and Surfaces: An Introduction Corner Cutting De Casteljau s and de Boor s algorithms all use corner-cutting procedures. Corner cutting can be local or non-local. A cut is local if it
More informationEstimating normal vectors and curvatures by centroid weights
Computer Aided Geometric Design 21 (2004) 447 458 www.elsevier.com/locate/cagd Estimating normal vectors and curvatures by centroid weights Sheng-Gwo Chen, Jyh-Yang Wu Department of Mathematics, National
More informationGrid Generation and Grid Conversion by Subdivision Schemes
Grid Generation and Grid Conversion by Subdivision Schemes Karl Heinz Brakhage Institute for Geometry and Applied Mathematics RWTH Aachen University D-55 Aachen brakhage@igpm.rwth-aachen.de Abstract In
More informationPhysically-Based Modeling and Animation. University of Missouri at Columbia
Overview of Geometric Modeling Overview 3D Shape Primitives: Points Vertices. Curves Lines, polylines, curves. Surfaces Triangle meshes, splines, subdivision surfaces, implicit surfaces, particles. Solids
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationBS-Patch: Constrained Bezier Parametric Patch
BS-Patch: Constrained Bezier Parametric Patch VACLAV SKALA, VIT ONDRACKA Department of Computer Science and Engineering University of West Bohemia Univerzitni 8, CZ 06 14 Plzen CZECH REPUBLIC skala@kiv.zcu.cz
More informationPhysics-Based Shape Modeling And Shape Recovery Using Multiresolution Subdivision Surfaces
Physics-Based Shape Modeling And Shape Recovery Using Multiresolution Subdivision Surfaces Chhandomay Mandal Λ Hong Qin y Baba C. Vemuri Λ Λ Department of Computer and Information Science and Engineering
More informationTriangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property
Triangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property Charles Loop cloop@microsoft.com February 1, 2001 Technical Report MSR-TR-2001-24 The masks for Loop s triangle subdivision
More informationApproximating Catmull-Clark Subdivision Surfaces with Bicubic Patches
Approximating Catmull-Clark Subdivision Surfaces with Bicubic Patches Charles Loop Microsoft Research Scott Schaefer Texas A&M University April 24, 2007 Technical Report MSR-TR-2007-44 Microsoft Research
More informationInterpolatory 3-Subdivision
EUROGRAPHICS 2000 / M. Gross and F.R.A. Hopgood (Guest Editors) Volume 19 (2000), Number 3 Interpolatory 3-Subdivision U. Labsik G. Greiner Computer Graphics Group University of Erlangen-Nuremberg Am Weichselgarten
More informationTernary Butterfly Subdivision
Ternary Butterfly Subdivision Ruotian Ling a,b Xiaonan Luo b Zhongxian Chen b,c a Department of Computer Science, The University of Hong Kong b Computer Application Institute, Sun Yat-sen University c
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationSmooth Patching of Refined Triangulations
Smooth Patching of Refined Triangulations Jörg Peters July, 200 Abstract This paper presents a simple algorithm for associating a smooth, low degree polynomial surface with triangulations whose extraordinary
More informationCurves and Surfaces 2
Curves and Surfaces 2 Computer Graphics Lecture 17 Taku Komura Today More about Bezier and Bsplines de Casteljau s algorithm BSpline : General form de Boor s algorithm Knot insertion NURBS Subdivision
More informationAn Efficient Data Structure for Representing Trilateral/Quadrilateral Subdivision Surfaces
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 3, No 3 Sofia 203 Print ISSN: 3-9702; Online ISSN: 34-408 DOI: 0.2478/cait-203-0023 An Efficient Data Structure for Representing
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 information1. Introduction. 2. Parametrization of General CCSSs. 3. One-Piece through Interpolation. 4. One-Piece through Boolean Operations
Subdivision Surface based One-Piece Representation Shuhua Lai Department of Computer Science, University of Kentucky Outline. Introduction. Parametrization of General CCSSs 3. One-Piece through Interpolation
More informationRecursive Subdivision Surfaces for Geometric Modeling
Recursive Subdivision Surfaces for Geometric Modeling Weiyin Ma City University of Hong Kong, Dept. of Manufacturing Engineering & Engineering Management Ahmad Nasri American University of Beirut, Dept.
More informationAdvanced Graphics. Subdivision Surfaces. Alex Benton, University of Cambridge Supported in part by Google UK, Ltd
Advanced Graphics Subdivision Surfaces Alex Benton, University of Cambridge A.Benton@damtp.cam.ac.uk Supported in part by Google UK, Ltd NURBS patches aren t the greatest NURBS patches are nxm, forming
More informationSurfaces for CAGD. FSP Tutorial. FSP-Seminar, Graz, November
Surfaces for CAGD FSP Tutorial FSP-Seminar, Graz, November 2005 1 Tensor Product Surfaces Given: two curve schemes (Bézier curves or B splines): I: x(u) = m i=0 F i(u)b i, u [a, b], II: x(v) = n j=0 G
More informationApproximate Geodesics on Smooth Surfaces of Arbitrary Topology
Approximate Geodesics on Smooth Surfaces of Arbitrary Topology Paper ID: 418 Category: Technical Paper The 6th International Symposium on Visual Computing (ISCV10) Las Vegas, Nevada, November 29 - December
More informationRemoving Polar Rendering Artifacts in Subdivision Surfaces
This is an electronic version of an article published in Journal of Graphics, GPU, and Game Tools, Volume 14, Issue 2 pp. 61-76, DOI: 10.1080/2151237X.2009.10129278. The Journal of Graphics, GPU, and Game
More informationLecture 2.2 Cubic Splines
Lecture. Cubic Splines Cubic Spline The equation for a single parametric cubic spline segment is given by 4 i t Bit t t t i (..) where t and t are the parameter values at the beginning and end of the segment.
More informationG 2 Bezier Crust on Quad Subdivision Surfaces
Pacific Graphics (2013) B. Levy, X. Tong, and K. Yin (Editors) Short Papers G 2 Bezier Crust on Quad Subdivision Surfaces paper 1348 Figure 1: Two examples of Bezier crust applied on Catmull-Clark subdivision
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 informationCurve Corner Cutting
Subdivision ision Techniqueses Spring 2010 1 Curve Corner Cutting Take two points on different edges of a polygon and join them with a line segment. Then, use this line segment to replace all vertices
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationNear-Optimum Adaptive Tessellation of General Catmull-Clark Subdivision Surfaces
Near-Optimum Adaptive Tessellation of General Catmull-Clark Subdivision Surfaces Shuhua Lai and Fuhua (Frank) Cheng (University of Kentucky) Graphics & Geometric Modeling Lab, Department of Computer Science,
More informationCS354 Computer Graphics Surface Representation IV. Qixing Huang March 7th 2018
CS354 Computer Graphics Surface Representation IV Qixing Huang March 7th 2018 Today s Topic Subdivision surfaces Implicit surface representation Subdivision Surfaces Building complex models We can extend
More informationSURFACE FAIRING FOR SHIP HULL DESIGN
SURFACE FAIRING FOR SHIP HULL DESIGN Xoán A. Leiceaga Eva Soto GED, Universidad de Vigo, Vigo, España leiceaga@uvigo.es Oscar E. Ruiz Carlos A. Vanegas Laboratorio CAD/CAM/CAE, Universidad EAFIT, Medellín,
More informationTiled Textures What if Miro Had Painted a Sphere
Tiled Textures What if Miro Had Painted a Sphere ERGUN AKLEMAN, AVNEET KAUR and LORI GREEN Visualization Sciences Program, Department of Architecture Texas A&M University December 26, 2005 Abstract We
More informationStudy on Improving the Quality of Reconstructed NURBS Surfaces
Study on Improving the Quality of Reconstructed NURBS Surfaces Shufeng jiang, Shigang Wang, Yong Yan School of Mechatronic Engineering, Qiqihar University, Qiqihar 161006, China Abstract In aspect of surface
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationHigh-Quality Shape Fitting for Aesthetic Modeling Based on Class A Condition
632 High-Quality Shape Fitting for Aesthetic Modeling Based on Class A Condition Tetsuo Oya 1, Hiroki Amemiya 2, Hideki Aoyama 3 and Masatake Higashi 4 1 Keio University, oya@sd.keio.ac.jp 2 Keio University,
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationApproximating Catmull-Clark Subdivision Surfaces with Bicubic Patches
Approximating Catmull-Clark Subdivision Surfaces with Bicubic Patches CHARLES LOOP Microsoft Research and SCOTT SCHAEFER Texas A&M University We present a simple and computationally efficient algorithm
More informationAn interpolating 4-point C 2 ternary stationary subdivision scheme
Computer Aided Geometric Design 9 (2002) 8 www.elsevier.com/locate/comaid An interpolating 4-point C 2 ternary stationary subdivision scheme M.F Hassan a,, I.P. Ivrissimitzis a, N.A. Dodgson a,m.a.sabin
More informationSubdivision surfaces for CAD: integration through parameterization and local correction
Workshop: New trends in subdivision and related applications September 4 7, 212 Department of Mathematics and Applications, University of Milano-Bicocca, Italy Subdivision surfaces for CAD: integration
More informationSurface Modeling. Polygon Tables. Types: Generating models: Polygon Surfaces. Polygon surfaces Curved surfaces Volumes. Interactive Procedural
Surface Modeling Types: Polygon surfaces Curved surfaces Volumes Generating models: Interactive Procedural Polygon Tables We specify a polygon surface with a set of vertex coordinates and associated attribute
More informationSmooth Surface Reconstruction using Doo-Sabin Subdivision Surfaces
Smooth Surface Reconstruction using Doo-Sabin Subdivision Surfaces Fuhua (Frank) Cheng, Fengtao Fan, Conglin Huang, Jiaxi Wang Department of Computer Science, University of Kentucky, Lexington, KY 40506,
More informationGeometric modeling 1
Geometric Modeling 1 Look around the room. To make a 3D model of a room requires modeling every single object you can see. Leaving out smaller objects (clutter) makes the room seem sterile and unrealistic
More informationParametrization of General Catmull-Clark Subdivision Surfaces and its Applications
Parametrization of General Catmull-Clark Subdivision Surfaces and its pplications Shuhua Lai and Fuhua (Frank) Cheng Graphics & Geometric Modeling Lab, Department of Computer Science University of Kentucky,
More informationSubdivision Surfaces. Homework 1: Last Time? Today. Bilinear Patch. Tensor Product. Spline Surfaces / Patches
Homework 1: Questions/Comments? Subdivision Surfaces Last Time? Curves & Surfaces Continuity Definitions Spline Surfaces / Patches Tensor Product Bilinear Patches Bezier Patches Trimming Curves C0, G1,
More informationCS354 Computer Graphics Surface Representation III. Qixing Huang March 5th 2018
CS354 Computer Graphics Surface Representation III Qixing Huang March 5th 2018 Today s Topic Bspline curve operations (Brief) Knot Insertion/Deletion Subdivision (Focus) Subdivision curves Subdivision
More informationExample: Loop Scheme. Example: Loop Scheme. What makes a good scheme? recursive application leads to a smooth surface.
Example: Loop Scheme What makes a good scheme? recursive application leads to a smooth surface 200, Denis Zorin Example: Loop Scheme Refinement rule 200, Denis Zorin Example: Loop Scheme Two geometric
More informationFall CSCI 420: Computer Graphics. 4.2 Splines. Hao Li.
Fall 2014 CSCI 420: Computer Graphics 4.2 Splines Hao Li http://cs420.hao-li.com 1 Roller coaster Next programming assignment involves creating a 3D roller coaster animation We must model the 3D curve
More informationTechnical Report. Removing polar rendering artifacts in subdivision surfaces. Ursula H. Augsdörfer, Neil A. Dodgson, Malcolm A. Sabin.
Technical Report UCAM-CL-TR-689 ISSN 1476-2986 Number 689 Computer Laboratory Removing polar rendering artifacts in subdivision surfaces Ursula H. Augsdörfer, Neil A. Dodgson, Malcolm A. Sabin June 2007
More informationHierarchical Grid Conversion
Hierarchical Grid Conversion Ali Mahdavi-Amiri, Erika Harrison, Faramarz Samavati Abstract Hierarchical grids appear in various applications in computer graphics such as subdivision and multiresolution
More informationINF3320 Computer Graphics and Discrete Geometry
INF3320 Computer Graphics and Discrete Geometry More smooth Curves and Surfaces Christopher Dyken, Michael Floater and Martin Reimers 10.11.2010 Page 1 More smooth Curves and Surfaces Akenine-Möller, Haines
More informationKnot Insertion and Reparametrization of Interval B-spline Curves
International Journal of Video&Image Processing and Network Security IJVIPNS-IJENS Vol:14 No:05 1 Knot Insertion and Reparametrization of Interval B-spline Curves O. Ismail, Senior Member, IEEE Abstract
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationOn Smooth Bicubic Surfaces from Quad Meshes
On Smooth Bicubic Surfaces from Quad Meshes Jianhua Fan and Jörg Peters Dept CISE, University of Florida Abstract. Determining the least m such that one m m bi-cubic macropatch per quadrilateral offers
More informationFINAL REPORT. Tessellation, Fairing, Shape Design, and Trimming Techniques for Subdivision Surface based Modeling
FINAL REPORT Tessellation, Fairing, Shape Design, and Trimming Techniques for Subdivision Surface based Modeling (DMS-0422126) PI: Fuhua (Frank) Cheng Department of Computer Science College of Engineering
More informationJoe Warren, Scott Schaefer Rice University
Joe Warren, Scott Schaefer Rice University Polygons are a ubiquitous modeling primitive in computer graphics. Their popularity is such that special purpose graphics hardware designed to render polygons
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationMatrix-valued 4-point Spline and 3-point Non-spline Interpolatory Curve Subdivision Schemes
Matrix-valued 4-point Spline and -point Non-spline Interpolatory Curve Subdivision Schemes Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis
More informationAn Algorithm of 3D Mesh Reconstructing Based on the Rendering Pipeline
3rd International Conference on Mechatronics and Information Technology (ICMIT 2016) An Algorithm of 3D Mesh Reconstructing Based on the Rendering Pipeline Zhengjie Deng1, a, Shuqian He1,b, Chun Shi1,c,
More informationNote on Industrial Applications of Hu s Surface Extension Algorithm
Note on Industrial Applications of Hu s Surface Extension Algorithm Yu Zang, Yong-Jin Liu, and Yu-Kun Lai Tsinghua National Laboratory for Information Science and Technology, Department of Computer Science
More informationNatural Quartic Spline
Natural Quartic Spline Rafael E Banchs INTRODUCTION This report describes the natural quartic spline algorithm developed for the enhanced solution of the Time Harmonic Field Electric Logging problem As
More informationPairs of Bi-Cubic Surface Constructions Supporting Polar Connectivity
Pairs of Bi-Cubic Surface Constructions Supporting Polar Connectivity Ashish Myles a, Kestutis Karčiauskas b Jörg Peters a a Department of CISE, University of Florida b Department of Mathematics and Informatics,
More informationAdvanced Geometric Modeling CPSC789
Advanced Geometric Modeling CPSC789 Fall 2004 General information about the course CPSC 789 Advanced Geometric Modeling Fall 2004 Lecture Time and Place ENF 334 TR 9:30 10:45 Instructor : Office: MS 618
More informationSubdivision overview
Subdivision overview CS4620 Lecture 16 2018 Steve Marschner 1 Introduction: corner cutting Piecewise linear curve too jagged for you? Lop off the corners! results in a curve with twice as many corners
More informationTo appear in Computer-Aided Design Revised June 18, J-splines
To appear in Computer-Aided Design Revised June 18, 2008 J-splines Jarek Rossignac School of Interactive Computing, College of Computing, Georgia Institute of Technology, Atlanta, GA http://www.gvu.gatech.edu/~jarek
More informationSubdivision Surfaces. Homework 1: Questions on Homework? Last Time? Today. Tensor Product. What s an illegal edge collapse?
Homework 1: Questions/Comments? Subdivision Surfaces Questions on Homework? Last Time? What s an illegal edge collapse? Curves & Surfaces Continuity Definitions 2 3 C0, G1, C1, C 1 a b 4 Interpolation
More informationModeling. Simulating the Everyday World
Modeling Simulating the Everyday World Three broad areas: Modeling (Geometric) = Shape Animation = Motion/Behavior Rendering = Appearance Page 1 Geometric Modeling 1. How to represent 3d shapes Polygonal
More informationFast Rendering of Subdivision Surfaces
Fast Rendering of Subdivision Surfaces Kari Pulli (Univ. of Washington, Seattle, WA) Mark Segal (SGI) Abstract Subdivision surfaces provide a curved surface representation that is useful in a number of
More informationSubdivision Surfaces. Homework 1: Questions/Comments?
Subdivision Surfaces Homework 1: Questions/Comments? 1 Questions on Homework? What s an illegal edge collapse? 1 2 3 a b 4 7 To be legal, the ring of vertex neighbors must be unique (have no duplicates)!
More information3D Modeling techniques
3D Modeling techniques 0. Reconstruction From real data (not covered) 1. Procedural modeling Automatic modeling of a self-similar objects or scenes 2. Interactive modeling Provide tools to computer artists
More informationProgressive Surface Modeling Based On 3D Motion Sketch
Progressive Surface Modeling Based On 3D Motion Sketch SHENGFENG IN, and DAVID K WRIGHT School of Engineering and Design Brunel University Uxbridge, Middlesex UB8 3PH UK Abstract: - This paper presents
More informationMathematical Tools in Computer Graphics with C# Implementations Table of Contents
Mathematical Tools in Computer Graphics with C# Implementations by Hardy Alexandre, Willi-Hans Steeb, World Scientific Publishing Company, Incorporated, 2008 Table of Contents List of Figures Notation
More informationA Heuristic Offsetting Scheme for Catmull-Clark Subdivision Surfaces
1 A Heuristic Offsetting Scheme for Catmull-Clark Subdivision Surfaces Jianzhong Wang 1, Fuhua Cheng 1,2 1 University of Kentucky, jwangf@uky.edu, cheng@cs.uky.edu 2 ational Tsing Hua University, cheng@cs.uky.edu
More informationQuasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters
Quasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters Reenu Sharma Assistant Professor, Department of Mathematics, Mata Gujri Mahila Mahavidyalaya, Jabalpur, Madhya Pradesh, India
More informationNeed for Parametric Equations
Curves and Surfaces Curves and Surfaces Need for Parametric Equations Affine Combinations Bernstein Polynomials Bezier Curves and Surfaces Continuity when joining curves B Spline Curves and Surfaces Need
More informationReparametrization of Interval Curves on Rectangular Domain
International Journal of Video&Image Processing and Network Security IJVIPNS-IJENS Vol:15 No:05 1 Reparametrization of Interval Curves on Rectangular Domain O. Ismail, Senior Member, IEEE Abstract The
More informationRefinable bivariate quartic and quintic C 2 -splines for quadrilateral subdivisions
Refinable bivariate quartic and quintic C 2 -splines for quadrilateral subdivisions Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis St. Louis,
More information