In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model f...In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.展开更多
In this paper, we use a univariate multiquadric quasi-interpolation scheme to solve the one-dimensional nonlinear sine-Gordon equation that is related to many physical phenomena. We obtain a numerical scheme by using ...In this paper, we use a univariate multiquadric quasi-interpolation scheme to solve the one-dimensional nonlinear sine-Gordon equation that is related to many physical phenomena. We obtain a numerical scheme by using the derivative of the quasi-interpolation to approximate the spatial derivative and a difference scheme to approximate the temporal derivative. The advantage of the obtained scheme is that the algorithm is very simple so that it is very easy to implement. The results of numerical experiments are presented and compared with analytical solutions to confirm the good accuracy of the presented scheme.展开更多
In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This st...In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.展开更多
Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the movi...Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.展开更多
Node interpolation cell method（NICM）is a micromechanics method employing the virtual displacement principle and the representative volume element（RVE）scheme to obtain the relationship between the global and the lo...Node interpolation cell method（NICM）is a micromechanics method employing the virtual displacement principle and the representative volume element（RVE）scheme to obtain the relationship between the global and the local strain.Mechanical properties of 2-D textile fabric reinforced ceramic matrix composites are predicted by NICM.Microstructures of 2-D woven and braided fabric reinforced composite are modeled by two kinds of RVE scheme.NICM is used to predict the macroscopic mechanical properties.The fill and warp yarns are simulated with cubic B-spline and their undulating forms are approximated by sinusoid.The effect of porosity on the fiber and matrix are considered as a reduction of elastic module.The connection of microstructure parameters and fiber volume fraction is modeled to investigate the reflection on the mechanical properties.The results predicted by NICM are compared with that by the finite element method（FEM）.The comparison shows that NICM is a valid and feasible method for predicting the mechanics properties of 2-D woven and braided fabric reinforced ceramic matrix composites.展开更多
A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines...A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.展开更多
In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been use...In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been used for this purpose.Traditional particle interpolation (TPI) is simple and easy to do, but its low accuracy has become an obstacle to its wider application.This can be seen in the cases of particle disorder arrangements and derivative calculations.There are many different methods to improve accuracy, with the moving least square (MLS) method one of the most important meshless interpolation methods.Unfortunately, it requires complex matrix computing and so is quite time-consuming.The authors developed a simpler scheme, called higher-order particle interpolation (HPI).This scheme can get more accurate derivatives than the MLS method, and its function value and derivatives can be obtained simultaneously.Although this scheme was developed for the SPH method, it has been found useful for other meshless methods.展开更多
The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the ...The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.展开更多
In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
An interpolating reproducing kernel particle method for two-dimensional （2D） scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in th...An interpolating reproducing kernel particle method for two-dimensional （2D） scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.展开更多
Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity proble...Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares （IMLS） method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin （IEFG） method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.展开更多
A B-spline active contour model based on finite element method is presented, into which the advantages of a B-spline active contour attributing to its fewer parameters and its smoothness is built accompanied with redu...A B-spline active contour model based on finite element method is presented, into which the advantages of a B-spline active contour attributing to its fewer parameters and its smoothness is built accompanied with reduced computational complexity and better numerical stability resulted from the finite element method. In this model, a cubic B-spline segment is taken as an element, and the finite element method is adopted to solve the energy minimization problem of the B-spline active contour, thus to implement image segmentation. Experiment results verify that this method is efficient for B-spline active contour, which attains stable, accurate and faster convergence.展开更多
In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation （RLW） is integrated fully by using an exponentia...In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation （RLW） is integrated fully by using an exponential B-spline Galerkin method in space together with Crank-Nicolson method in time. Three numerical examples related to propagation of sin- gle solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.展开更多
Making an exact computation of added resistance in sea waves is of high interest due to the economic effects relating to ship design and operation. In this paper, a B-spline based method is developed for computation o...Making an exact computation of added resistance in sea waves is of high interest due to the economic effects relating to ship design and operation. In this paper, a B-spline based method is developed for computation of added resistance. Based on the potential flow assumption, the velocity potential is computed using Green＇s formula. The Kochin function is applied to compute added resistance using Maruo＇s far-field method, the body surface is described by a B-spline curve and potentials and normal derivation of potentials are also described by B-spline basis functions and B-spline derivations. A collocation approach is applied for numerical computation, and integral equations are then evaluated by applying Gauss–Legendre quadrature. Computations are performed for a spheroid and different hull forms; results are validated by a comparison with experimental results. All results obtained with the present method show good agreement with experimental results.展开更多
The work studies model reduction method for nonlinear systems based on proper orthogonal decomposition (POD)and discrete empirical interpolation method (DEIM). Instead of using the classical DEIM to directly approxima...The work studies model reduction method for nonlinear systems based on proper orthogonal decomposition (POD)and discrete empirical interpolation method (DEIM). Instead of using the classical DEIM to directly approximate thenonlinear term of a system, our approach extracts the main part of the nonlinear term with a linear approximation beforeapproximating the residual with the DEIM. We construct the linear term by Taylor series expansion and dynamic modedecomposition (DMD), respectively, so as to obtain a more accurate reconstruction of the nonlinear term. In addition, anovel error prediction model is devised for the POD-DEIM reduced systems by employing neural networks with the aid oferror data. The error model is cheaply computable and can be adopted as a remedy model to enhance the reduction accuracy.Finally, numerical experiments are performed on two nonlinear problems to show the performance of the proposed method.展开更多
A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalize...A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalized inductiveelectromotive force was investigated. According to the turning point, the transient process is divided into the earlyphase, the turning point, and the late phase. Afterwards, apparent resistivity is obtained through inverse spline interpo-lation in the early and the late phases, respectively. Finally, the resistivities of the early-time and the late-time wereconnected together by the turning point. The result shows that the inverse spline method is feasible and the method alsolays a foundation for initial model construction in the TEM automatic inversion.展开更多
基金This work was supported by the National Natural Science Foundation of China(10071037)
文摘In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.
基金supported by the State Key Development Program for Basic Research of China (Grant No 2006CB303102)Science and Technology Commission of Shanghai Municipality,China (Grant No 09DZ2272900)
文摘In this paper, we use a univariate multiquadric quasi-interpolation scheme to solve the one-dimensional nonlinear sine-Gordon equation that is related to many physical phenomena. We obtain a numerical scheme by using the derivative of the quasi-interpolation to approximate the spatial derivative and a difference scheme to approximate the temporal derivative. The advantage of the obtained scheme is that the algorithm is very simple so that it is very easy to implement. The results of numerical experiments are presented and compared with analytical solutions to confirm the good accuracy of the presented scheme.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China(Grant No.10902076)the Natural Science Foundation of Shanxi Province of China(Grant No.2007011009)+1 种基金the Scientific Research and Development Program of the Shanxi Higher Education Institutions(Grant No.20091131)the Doctoral Startup Foundation of Taiyuan University of Science and Technology(Grant No.200708)
文摘In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Natural Science Foundation of Ningbo City,China(GrantNo.2013A610103)+2 种基金the Natural Science Foundation of Zhejiang Province,China(Grant No.Y6090131)the Disciplinary Project of Ningbo City,China(GrantNo.SZXL1067)the K.C.Wong Magna Fund in Ningbo University,China
文摘Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.
基金Supported by the Aviation Science Foundationof China(2009ZB5052)the Specialized Research Foundation for the Doctor Program of Higher Education(20070287039)~~
文摘Node interpolation cell method（NICM）is a micromechanics method employing the virtual displacement principle and the representative volume element（RVE）scheme to obtain the relationship between the global and the local strain.Mechanical properties of 2-D textile fabric reinforced ceramic matrix composites are predicted by NICM.Microstructures of 2-D woven and braided fabric reinforced composite are modeled by two kinds of RVE scheme.NICM is used to predict the macroscopic mechanical properties.The fill and warp yarns are simulated with cubic B-spline and their undulating forms are approximated by sinusoid.The effect of porosity on the fiber and matrix are considered as a reduction of elastic module.The connection of microstructure parameters and fiber volume fraction is modeled to investigate the reflection on the mechanical properties.The results predicted by NICM are compared with that by the finite element method（FEM）.The comparison shows that NICM is a valid and feasible method for predicting the mechanics properties of 2-D woven and braided fabric reinforced ceramic matrix composites.
基金supported by the National Natural Science Foundation of China(11871312,12131014)the Natural Science Foundation of Shandong Province,China(ZR2023MA086)。
文摘A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.
基金Supported by the National Natural Science Foundation of China under Grant No.10572041,50779008Doctoral Fund of Ministry of Education of China under Grant No.20060217009
文摘In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been used for this purpose.Traditional particle interpolation (TPI) is simple and easy to do, but its low accuracy has become an obstacle to its wider application.This can be seen in the cases of particle disorder arrangements and derivative calculations.There are many different methods to improve accuracy, with the moving least square (MLS) method one of the most important meshless interpolation methods.Unfortunately, it requires complex matrix computing and so is quite time-consuming.The authors developed a simpler scheme, called higher-order particle interpolation (HPI).This scheme can get more accurate derivatives than the MLS method, and its function value and derivatives can be obtained simultaneously.Although this scheme was developed for the SPH method, it has been found useful for other meshless methods.
基金supported by the Key Program of the National Natural Science Foundation of China (Grand No. 51138001)the China-German Cooperation Project (Grand No. GZ566)+1 种基金the Innovative Research Groups Funded by the National Natural Science Foundation of China (Grand No. 51121005)the Special Funds for the Basic Scientific Research Expenses for the Central University (Grant No. DUT13LK16)
文摘The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金supported by the National Natural Science Foundation of China(Grant No.11171208)the Natural Science Foundation of Shanxi Province,China(Grant No.2013011022-6)
文摘An interpolating reproducing kernel particle method for two-dimensional （2D） scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)
文摘Based on the improved interpolating moving least-squares （ⅡMLS） method and the Galerkin weak form, an improved interpolating element-free Galerkin （ⅡEFG） method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares （IMLS） method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin （IEFG） method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.
基金the National Natural Science Foundation of China (No.59975057).
文摘A B-spline active contour model based on finite element method is presented, into which the advantages of a B-spline active contour attributing to its fewer parameters and its smoothness is built accompanied with reduced computational complexity and better numerical stability resulted from the finite element method. In this model, a cubic B-spline segment is taken as an element, and the finite element method is adopted to solve the energy minimization problem of the B-spline active contour, thus to implement image segmentation. Experiment results verify that this method is efficient for B-spline active contour, which attains stable, accurate and faster convergence.
基金supported by the Scientific and Technological Research Council of Turkey(Grant No.113F394)
文摘In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation （RLW） is integrated fully by using an exponential B-spline Galerkin method in space together with Crank-Nicolson method in time. Three numerical examples related to propagation of sin- gle solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.
文摘Making an exact computation of added resistance in sea waves is of high interest due to the economic effects relating to ship design and operation. In this paper, a B-spline based method is developed for computation of added resistance. Based on the potential flow assumption, the velocity potential is computed using Green＇s formula. The Kochin function is applied to compute added resistance using Maruo＇s far-field method, the body surface is described by a B-spline curve and potentials and normal derivation of potentials are also described by B-spline basis functions and B-spline derivations. A collocation approach is applied for numerical computation, and integral equations are then evaluated by applying Gauss–Legendre quadrature. Computations are performed for a spheroid and different hull forms; results are validated by a comparison with experimental results. All results obtained with the present method show good agreement with experimental results.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11871400 and 11971386)the Natural Science Foundation of Shaanxi Province,China(Grant No.2017JM1019).
文摘The work studies model reduction method for nonlinear systems based on proper orthogonal decomposition (POD)and discrete empirical interpolation method (DEIM). Instead of using the classical DEIM to directly approximate thenonlinear term of a system, our approach extracts the main part of the nonlinear term with a linear approximation beforeapproximating the residual with the DEIM. We construct the linear term by Taylor series expansion and dynamic modedecomposition (DMD), respectively, so as to obtain a more accurate reconstruction of the nonlinear term. In addition, anovel error prediction model is devised for the POD-DEIM reduced systems by employing neural networks with the aid oferror data. The error model is cheaply computable and can be adopted as a remedy model to enhance the reduction accuracy.Finally, numerical experiments are performed on two nonlinear problems to show the performance of the proposed method.
基金Project 40344022 supported by National Natural Science Foundation of China
文摘A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalized inductiveelectromotive force was investigated. According to the turning point, the transient process is divided into the earlyphase, the turning point, and the late phase. Afterwards, apparent resistivity is obtained through inverse spline interpo-lation in the early and the late phases, respectively. Finally, the resistivities of the early-time and the late-time wereconnected together by the turning point. The result shows that the inverse spline method is feasible and the method alsolays a foundation for initial model construction in the TEM automatic inversion.
