a hamiltonian preserving discontinuous galerkin method for the generalized korteweg de vries equation the invariant preserving property is one of the guiding principles for numerical algorithms in solving wave equations in order to minimize phase and amplitude errors after long time simulation in this paper we design analyze and numerically validate a hamiltonian preserving discontinuous galerkin method for solving the korteweg de vries kdv equation for the generalized kdv equation the semi discrete formulation is shown to preserve both the first and the third conserved integrals and approximately preserve the second conserved integral for the linearized kdv equation all the first three conserved integrals are preserved and optimal error estimates are obtained for polynomials of even degree the preservation properties are also maintained by the fully discrete dg scheme our numerical experiments demonstrate both high accuracy of convergence and preservation of all three conserved integrals for the generalized kdv equation we also show that the shape of the solution after long time simulation is well preserved due to the hamiltonian preserving property

on the numerical experiments of the cauchy problem for semi linear klein gordon equations in the de sitter spacetime abstract the computational analysis of the cauchy problem for semi linear klein gordon equations in the de sitter spacetime is considered several simulations are performed to show the time global behaviors of the solutions of the equations in the spacetime based on the structure preserving scheme it is remarked that the sufficiently large hubble constant yields the strong diffusion effect which gives the long and stable simulations for the defocusing semi linear terms the reliability of the simulations is confirmed by the preservation of the numerically modified hamiltonian of the equations

a dynamically consistent method to solve nonlinear multidimensional advection reaction equations with fractional diffusion abstract this work is motivated by an extension of both the burgers fisher and the burgers huxley equations in multiple dimensions considering riesz fractional diffusion initial boundary conditions which are positive and bounded are imposed on a closed and bounded set and a finite difference method is proposed to approximate the solutions of the fractional model the methodology is a linear and implicit technique which is based on fractional centered differences we show in this manuscript that the method can be expressed in vector form using a minkowski matrix under suitable conditions the main properties of minkowski matrices are used then to establish the existence and the uniqueness of the solutions of the finite difference method as well as the capability of the technique to preserve the positivity and the boundedness additionally we show that the method is a consistent technique which is stable and convergent with first order of convergence in time and second order in space some illustrative simulations show that the scheme is capable of preserving the positivity and the boundedness of the numerical approximations

a moderate deviation principle for 2 d stochastic navier stokes equations in this paper we prove a central limit theorem and establish a moderate deviation principle for twodimensional stochastic navier stokes equations with multiplicative noise the weak convergence method plays an important role

conservative nonlinear difference scheme for the cahn hilliard equation ii abstract numerical solutions for the cahn hilliard equation is considered using the crank nicolson type finite difference method existence of the solution for the difference scheme has been shown by brouwer fixed point theorem stability convergence and error analysis of the scheme are shown we also show that the scheme preserves the discrete mass even though the linearized scheme in 1 is conditionally stable and does not preserve the mass

bracket formulations and energy and helicity preserving numerical methods for the three dimensional vorticity equation abstract the vorticity equation for three dimensional viscous incompressible fluid flows is formulated within different bracket formalisms using the poisson or nambu bracket together with a dissipative bracket the budgets of kinetic energy helicity and enstrophy derived from the bracket formulations are properly inherited by the finite difference equations obtained by invoking the discrete variational derivative method combined with the mimetic finite difference method in particular energy and helicity are conserved precisely in inviscid flow computations the energy and enstrophy dissipate properly owing to viscosity in viscous flow computations and the enstrophy is appropriately produced by the vortex stretching effect in both inviscid and viscous flow computations the relationships between the stream function velocity and vorticity as well as the solenoidal conditions on the velocity and vorticity fields are also inherited numerical experiments on a periodic array of rolls that permits analytical solutions have been done to examine the properties and usefulness of the proposed method

a roadmap for discretely energy stable schemes for dissipative systems based on a generalized auxiliary variable with guaranteed positivity abstract we present a framework for devising discretely energy stable schemes for general dissipative systems based on a generalized auxiliary variable the auxiliary variable a scalar number can be defined in terms of the energy functional by a general class of functions not limited to the square root function adopted in previous approaches the current method has another remarkable property the computed values for the generalized auxiliary variable are guaranteed to be positive on the discrete level regardless of the time step sizes or the external forces this property of guaranteed positivity is not available in previous approaches a unified procedure for treating the dissipative governing equations and the generalized auxiliary variable on the discrete level has been presented the discrete energy stability of the proposed numerical scheme and the positivity of the computed auxiliary variable have been proved for general dissipative systems the current method termed gpav generalized positive auxiliary variable requires only the solution of linear algebraic equations within a time step with appropriate choice of the operator in the algorithm the resultant linear algebraic systems upon discretization involve only constant and time independent coefficient matrices which only need to be computed once and can be pre computed several specific dissipative systems are studied in relative detail using the gpav framework ample numerical experiments are presented to demonstrate the performance of the method and the robustness of the scheme at large time step sizes

numerically induced phase shift in the kdv soliton abstract when using a finite difference scheme to study the motion of a kdv solition a shift in the position of the solition from the exact solution is detected in this paper we retain the lowest order terms in the truncation error and treat them analytically as a perturbation of the kdv equation it is found that perturbation theory can be used to determine the numerically induced shift

a robust numerical integrator for the short pulse equation near criticality abstract the short pulse equation was introduced by schafer wayne 2004 for modeling the propagation of ultrashort optical pulses while it can describe a wide range of solutions its ultrashort pulse solutions with a few cycles which the conventional nonlinear schrodinger equation does not possess have drawn much attention in such a region the solutions can become quite close to a singularity and thus existing numerical methods cease to work stably or even if they do they require high computational cost in this paper we propose a robust numerical integration method which is obtained by combining a hodograph transformation and some structure preserving methods the resulting scheme successfully works even when the singularity occurs and thus gives a highly robust method for resolving near singular solutions of interest it is confirmed by numerical experiments and some new insights about derivative blow up are obtained

a numerical energy conserving method for the dnls equation an implicit numerical energy conserving method is developed for the derivative nonlinear schrodinger dnls equation for periodic boundary conditions we find no numerical high frequency modulational instabilities in addition to the modulational instability from a linear analysis around a nonlinear state for the dnls equation if the modulation is small and k0 a22 2 t k0 is the wavenumber and a the amplitude the numerical scheme is used to follow the nonlinear behavior of the dnls modulational instability the numerical code is also tested by the evolution for one soliton initial data these tests show that if the modulation is not small compared to the background wave amplitude new nonlinear numerical instabilities are introduced

a second order splitting method for the cahn hilliard equation a semi discrete finite element method requiring only continuous element is presented for the approximation of the solution of the evolutionary fourth order in space cahn hilliard equation optimal order error bounds are derived in various norms for an implementation which uses mass lumping the continuous problem has an energy based lyapunov functional it is proved that this property holds for the discrete problem

conservative numerical methods for x f x abstract two distinctly different numerical methods are developed for solving the conservative initial value problem x f x x o dot x o both methods conserve exactly the same total energy as does the given differential equation computer examples are described and

conservative numerical schemes for the ostrovsky equation the ostrovsky equation describes gravity waves under the influence of coriolis force it is known that solutions of this equation conserve the l 2 norm and an energy function that is determined non locally in this paper we propose four conservative numerical schemes for this equation a finite difference scheme and a pseudospectral scheme that conserve the norm and the same types of schemes that conserve the energy a numerical comparison of these schemes is also provided which indicates that the energy conservative schemes perform better than the norm conservative schemes

efficient local energy dissipation preserving algorithms for the cahn hilliard equation abstract in this paper we show that the cahn hilliard equation possesses a local energy dissipation law which is independent of boundary conditions and produces much more information of the original problem to inherit the intrinsic property we derive three novel local structure preserving algorithms for the 2d cahn hilliard equation by the concatenating method in particular when the nonlinear bulk potential f in the equation is chosen as the ginzburg landau double well potential the method discussed by zhang and qiao 2012 50 is a special case of our scheme ii thanks to the leibnitz rules and properties of operators the three schemes are rigorously proven to conserve the discrete local energy dissipation law in any local time space region under periodic boundary conditions the schemes are proven to possess the discrete mass conservation and total energy dissipation laws numerical experiments are conducted to show the performance of the proposed schemes

moderate deviations of minimum contrast estimators under contamination since statistical models are simplifications of reality it is important in estimation theory to study the behavior of estimators also under distributions slightly different from the proposed model in testing theory when dealing with test statistics where nuisance parameters are estimated knowledge of the behavior of the estimators of the nuisance parameters is needed under alternatives to evaluate the power in this paper the moderate deviation behavior of minimum contrast estimators is investigated not only under the supposed model but also under distributions close to the model a particular example is the multivariate maximum likelihood estimator determined within the proposed model the set up is quite general including for instance also discrete distributions the rate of convergence under alternatives is determined both when comparing the minimum contrast estimator with a natural parameter in the parameter space and when comparing it with the proposed true value in the parameter space it turns out that under the model the asymptotic optimality of the maximum likelihood estimator in the local sense continues to hold in the moderate deviation area

a second order accurate linearized difference scheme for the two dimensional cahn hilliard equation the cahn hilliard equation is a nonlinear evolutionary equation that is of fourth order in space in this paper a linearized finite difference scheme is derived by the method of reduction of order it is proved that the scheme is uniquely solvable and convergent with the convergence rate of order two in a discrete l2 norm the coefficient matrix of the difference system is symmetric and positive definite so many well known iterative methods e g gauss seidel sor can be used to solve the system

dissipative conservative galerkin method using discrete partial derivatives for nonlinear evolution equations a new method is proposed for designing galerkin schemes that retain the energy dissipation or conservation properties of nonlinear evolution equations such as the cahn hilliard equation the korteweg de vries equation or the nonlinear schrodinger equation in particular as a special case dissipative or conservative finite element schemes can be derived the key device there is the new concept of discrete partial derivatives as examples of the application of the present method dissipative or conservative galerkin schemes are presented for the three equations with some numerical experiments

notes on finite difference schemes to a parabolic elliptic system modelling chemotaxis we consider numerical solutions to a nonlinear parabolic elliptic system in mathematical biology an important property of solutions to the system is the conservation of l 1 norm applying the upwind technique we make a finite difference scheme satisfying the conservation of a discrete l 1 norm we also describe a remark on time discretization from the view point of lyapunov s property which is another important feature of the system several numerical examples clarify the nature of our scheme

numerical integrators that preserve symmetries and reversing symmetries we consider properties of flows the relationships between them and whether numerical integrators can be made to preserve these properties this is done in the context of automorphisms and antiautomorphisms of a certain group generated by maps associated to vector fields this new framework unifies several known constructions we also use the concept of covariance of a numerical method with respect to a group of coordinate transformations the main application is to explore the relationship between spatial symmetries reversing symmetries and time symmetry of flows and numerical integrators

the fourier collocation method for the cahn hilliard equation abstract in this paper a fourier collocation method for numerically solving cahn hilliard equations with periodic boundary conditions is developed we establish their semidiscrete and fully discrete schemes that inherit the energy dissipation property and mass conservation property from the associated continuous problem we prove existence and uniqueness of the numerical solution and derive the optimal error bounds we perform some numerical experiments which confirm our results

direct discretization of planar div curl problems a control volume method is proposed for planar div curl systems the method is independent of potential and least squares formulations and works directly with the div curl system the novelty of the technique lies in its use of a single local vector field component and two control volumes rather than the other way round a discrete vector field theory comes quite naturally from this idea and is developed in the paper error estimates are proved for the method and other ramifications investigated

data driven smooth tests when the hypothesis is composite in recent years several authors have recommended smooth tests for testing goodness of fit however the number of components in the smooth test statistic should be chosen well otherwise considerable loss of power may occur schwarz s selection rule provides one such good choice earlier results on simple null hypotheses are extended here to composite hypotheses which tend to be of more practical interest for general composite hypotheses consistency of the data driven smooth tests holds at essentially any alternative monte carlo experiments on testing exponentiality and normality show that the data driven version of neyman s test compares well to other even specialized tests over a wide range of alternatives

data driven smooth tests for composite hypotheses the classical problem of testing goodness of fit of a parametric family is reconsidered a new test for this problem is proposed and investigated the new test statistic is a combination of the smooth test statistic and schwarz s selection rule more precisely as the sample size increases an increasing family of exponential models describing departures from the null model is introduced and schwarz s selection rule is presented to select among them schwarz s rule provides the right dimension given by the data while the smooth test in the right dimension finishes the job theoretical properties of the selection rules are derived under null and alternative hypotheses they imply consistency of data driven smooth tests for composite hypotheses at essentially any alternative

a hamiltonian conserving galerkin scheme for the camassa holm equation a new hamiltonian conserving galerkin scheme for the camassa holm equation is presented the scheme has an additional welcome feature that in exact arithmetic it is unconditionally stable in the sense that the solution is always bounded numerical examples that confirm the theory and the effectiveness of the scheme are also given

analytical and numerical aspects of certain nonlinear evolution equations i analytical abstract nonlinear partial difference equations are obtained which have as limiting forms the nonlinear schrodinger korteweg devries and modified korteweg devries equations these difference equations have a number of special properties they are constructed by methods related to the inverse scattering transform they can be used as a basis for numerical schemes to the associated nonlinear evolution equations experiments have shown that they compare very favorably with other known numerical methods papers ii iii in paper ii the ablowitz ladik scheme for the nonlinear schrodinger equation is compared to other known numerical schemes and generally proved to be faster than all utilized finite difference schemes but somewhat slower than the finite fourier pseudospectral methods in paper iii a proposed scheme for the korteweg de vries equation proved to be faster than both the finite difference and finite fourier methods already considered

conservative finite difference schemes for the degasperis procesi equation we consider the numerical integration of the degasperis procesi equation which was recently introduced as a completely integrable shallow water equation for the equation we propose nonlinear and linear finite difference schemes that preserve two invariants associated with the bi hamiltonian form of the equation at the same time we also prove the unique solvability of the schemes and show some numerical examples

statistical inference for the intensity in a partially observed jump diffusion abstract this article focuses on a partially observed linear diffusion with jumps described by a poisson process precisely we study an inferential problem for the intensity of the poisson process by establishing a moderate deviation principle for the least squares estimator of the intensity

an extension of the discrete variational method to nonuniform grids the discrete variational method is a method used to derive finite difference schemes that inherit the conservation dissipation property of the original equations although this method has mainly been developed for uniform grids we extend this method to multidimensional nonuniform meshes

the discrete variational derivative method based on discrete differential forms as is well known for pdes that enjoy a conservation or dissipation property numerical schemes that inherit this property are often advantageous in that the schemes are fairly stable and give qualitatively better numerical solutions in practice lately furihata and matsuo have developed the so called discrete variational derivative method that automatically constructs energy preserving or dissipative finite difference schemes although this method was originally developed on uniform meshes the use of non uniform meshes is of importance for multi dimensional problems on the other hand the theories of discrete differential forms have received much attention recently these theories provide a discrete analogue of the vector calculus on general meshes in this paper we show that the discrete variational derivative method and the discrete differential forms by bochev and hyman can be combined applications to the cahn hilliard equation and the klein gordon equation on triangular meshes are provided as demonstrations we also show that the schemes for these equations are h1 stable under some assumptions in particular one for the nonlinear klein gordon equation is obtained by combination of the energy conservation property and the discrete poincare inequality which are the temporal and spacial structures that are preserved by the above methods

continuous finite element methods which preserve energy properties for nonlinear problems abstract continuous finite elements are used to develop numerical methods which preserve energy properties for nonlinear ordinary and partial differential equations schemes are constructed that give high orders of accuracy time discretization algorithms for the cahn hilliard klein gordon korteweg de vries and schrodinger equations are presented comparisons are made with existing schemes

symmetry preserving discretization of turbulent flow we propose to perform turbulent flow simulations in such manner that the difference operators do have the same symmetry properties as the underlying differential operators i e the convective operator is represented by a skew symmetric coefficient matrix and the diffusive operator is approximated by a symmetric positive definite matrix mimicking crucial properties of differential operators forms in itself a motivation for discretizing them in a certain manner we give it a concrete form by noting that a symmetry preserving discretization of the navier stokes equations is stable on any grid and conserves the total mass momentum and kinetic energy for the latter the physical dissipation is to be turned off of coarse being stable on any grid the choice of the grid may be based on the required accuracy solely and the main question becomes how accurate is a symmetry preserving discretization its accuracy is tested for a turbulent flow in a channel by comparing the results to those of physical experiments and previous numerical studies the comparison is carried out for a reynolds number of 5600 which is based on the channel width and the mean bulk velocity based on the channel half width and wall shear velocity the reynolds number becomes 180 the comparison shows that with a fourth order symmetry preserving method a 64 64 32 grid suffices to perform an accurate numerical simulation

an explicit finite difference scheme with exact conservation properties abstract a finite difference scheme for the numerical study of the korteweg de vries equation is constructed it is explicit and yet conserves exactly the energy of the computed solutions the underlying idea can also be applied to more general equations or systems numerical experiments are included

integrable discretization of time and its application on the fourier pseudospectral method to the korteweg de vries equation this paper presents a new integrable discretization of the korteweg de vries kdv equation different from other discrete analogues we discretize the variable time and obtain an integrable differential difference system this system has the original kdv equation as a standard limit when the step size tends to zero the main idea is based on hirota s bilinear method and backlund transformation and can be applied to other integrable systems by applying the fourier pseudospectral method to the space variable we derive a new numerical scheme for the kdv equation numerical results are found to agree with the exact solution and the first five conservation quantities are preserved quite well

a least squares finite element method for the numerical solution of the navier stokes cahn hilliard system modeling the motion of the contact line in this article we discuss the numerical solution of the navier stokes cahn hilliard system modeling the motion of the contact line separating two immiscible incompressible viscous fluids near a solid wall the method we employ combines a finite element space approximation with a time discretization by operator splitting to solve the cahn hilliard part of the problem we use a least squares conjugate gradient method we also show that the scheme has the total energy decaying in time property under certain conditions our numerical experiments indicate that the method discussed here is accurate stable and efficient

uniform stabilization of numerical schemes for the critical generalized korteweg de vries equation with damping this work is devoted to the analysis of a fully implicit numerical scheme for the critical generalized korteweg de vries equation gkdv with p 4 in a bounded domain with a localized damping term the damping is supported in a subset of the domain so that the solutions of the continuous model issuing from small data are globally defined and exponentially decreasing in the energy space based in this asymptotic behavior of the solution we introduce a finite difference scheme which despite being one of the first order has the good property to converge in l 4 strong combining this strong convergence with discrete multipliers and a contradiction argument we show that the smallness of the initial condition leads to the uniform with respect to the mesh size exponential decay of the energy associated to the scheme numerical experiments are provided to illustrate the performance of the method and to confirm the theoretical results

regular article finite difference schemes for u t x α δg δu that inherit energy conservation or dissipation property we propose a new procedure for designing by rote finite difference schemes that inherit energy conservation or dissipation property from nonlinear partial differential equations such as the korteweg de vries kdv equation and the cahn hilliard equation the most important feature of our procedure is a rigorous discretization of variational derivatives using summation by parts which implies that the inherited properties are satisfied exactly since the inherited properties are kept even if the time mesh size changes in the time evolution process we can use some appropriate time mesh adaptive methods to obtain numerical solutions through the derived schemes because of these properties the derived schemes are expected to be numerically stable and yield solutions converging to pde solutions and sufficiently flexible to treat the inheritance of the energy conservation and dissipation properties are verified numerically for the kdv equation and the cahn hilliard equation

on the gpgpu parallelization issues of finite element approximate inverse preconditioning during the past decades explicit finite element approximate inverse preconditioning methods have been extensively used for efficiently solving sparse linear systems on multiprocessor systems the effectiveness of explicit approximate inverse preconditioning schemes relies on the use of efficient preconditioners that are close approximants to the coefficient matrix and are fast to compute in parallel new parallel computational techniques are proposed for the parallelization of the optimized banded generalized approximate inverse finite element matrix obgaifem algorithm based on the concept of the fish bone computational approach and for the explicit preconditioned conjugate gradient type methods on a general purpose graphics processing unit gpgpu the proposed parallel methods have been implemented using compute unified device architecture cuda developed by nvidia finally numerical results for the performance of the finite element explicit approximate inverse preconditioning for solving characteristic two dimensional boundary value problems on a massive multiprocessor interface on a gpu are presented the cuda implementation issues of the proposed methods are also discussed

dissipative or conservative finite difference schemes for complex valued nonlinear partial differential equations abstract we propose a new procedure for designing finite difference schemes that inherit energy conservation or dissipation property from complex valued nonlinear partial differential equations pdes such as the nonlinear schrodinger equation the ginzburg landau equation and the newell whitehead equation the procedure is a complex version of the procedure that furihata has recently presented for real valued nonlinear pdes furthermore we show that the proposed procedure can be modified for designing linearly implicit finite difference schemes that inherit energy conservation or dissipation property

a general framework for deriving integral preserving numerical methods for pdes a general procedure for constructing conservative numerical integrators for time dependent partial differential equations is presented in particular linearly implicit methods preserving a time discretized version of the invariant are developed for systems of partial differential equations with polynomial nonlinearities the framework is rather general and allows for an arbitrary number of dependent and independent variables with derivatives of any order it is proved formally that second order convergence is obtained the procedure is applied to a test case and numerical experiments are provided

numerical analysis of a continuum model of phase transition a method for solution of the cahn hilliard equation is presented unlike previous work the discrete equations for the new method possess a lyapunov function this makes it possible to prove convergence of the approximate solutions without assumptions beyond those necessary for existence and uniqueness of the differential equation several consequences are explored

new conservative difference schemes for a coupled nonlinear schrodinger system abstract in this paper two conservative difference schemes for solving a coupled nonlinear schrodinger cnls system are numerically analyzed firstly a nonlinear implicit two level finite difference scheme for cnls system is studied then a linear three level difference scheme for cnls system is presented an induction argument and the discrete energy method are used to prove the second order convergence and unconditional stability of the linear scheme numerical examples show the efficiency of the new scheme

delta method in large deviations and moderate deviations for estimators the delta method is a popular and elementary tool for deriving limiting distributions of transformed statistics while applications of asymptotic distributions do not allow one to obtain desirable accuracy of approximation for tail probabilities the large and moderate deviation theory can achieve this goal motivated by the delta method in weak convergence a general delta method in large deviations is proposed the new method can be widely applied to driving the moderate deviations of estimators and is illustrated by examples including the wilcoxon statistic the kaplan meier estimator the empirical quantile processes and the empirical copula function we also improve the existing moderate deviations results for m estimators and l statistics by the new method some applications of moderate deviations to statistical hypothesis testing are provided

finite difference calculus invariant structure of a class of algorithms for the nonlinear klein gordon equation in a previous work the authors have presented a formalism for deriving systematically invariant symmetric finite difference algorithms for nonlinear evolution differential equations that admit conserved quantities this formalism is herein cast in the context of exact finite difference calculus the algorithms obtained from the proposed formalism are shown to derive exactly from discrete scalar potential functions using finite difference calculus in the same sense as that of the corresponding differential equation being derivable from its associated energy function a conserved quantity a clear ramification of this result is that the derived algorithms preserve certain discrete invariant quantities which are the consistent counterpart of the invariant quantities in the continuous case results on the nonlinear stability of a class of algorithms that are derived using the proposed formalism and that preserve energy or linear momentum are discussed in the context of finite difference calculus some

a moderate deviation principle for 2 d stochastic navier stokes equations driven by multiplicative levy noises abstract in this paper we establish a moderate deviation principle for two dimensional stochastic navier stokes equations driven by multiplicative levy noises the weak convergence method introduced by budhiraja dupuis and ganguly in 3 plays a key role

accuracy and conservation properties in numerical integration the case of the korteweg de vries equation when numerically integrating time dependent differential equations it is often recommended to employ methods that preserve some of the invariant quantities mass energy etc of the problem being considered this recommendation is usually justified on the grounds that conservation of invariant quantities may ensure that the numerical solution possesses some important qualitative features however there are cases where schemes that preserve invariants are also advantageous in that they possess favourable error propagation mechanisms that render them superior from a quantitative point of view in the present paper we consider the korteweg de vries equation as a case study we show rigorously that for soliton problems and at leading order the error of conservative schemes consists of a phase error that grows linearly with time plus a complementary term that is bounded in the h 1 norm uniformly in time for general nonconservative schemes the error involves a linearly growing amplitude error a quadratically growing phase error and a complementary term that grows linearly in the h 1 norm numerical experiments are presented

analytical and numerical aspects of certain nonlinear evolution equations 1v numerical modified korteweg de vries equation various numerical methods are used in order to approximate the korteweg de vries equation namely i zabusky kruskal scheme ii hopscotch method iii a scheme due to goda iv a proposed local scheme v a proposed global scheme vi a scheme suggested by kruskal vii split step fourier method by tappert viii an improved split step fourier method and ix pseudospectral method by fornberg and whitham comparisons between our proposed scheme which is developed using notions of the inverse scattering transform and the other utilized schemes are obtained

mimetic discretizations for maxwell s equations we have constructed reliable finite difference methods for approximating the solution to maxwell s equations using accurate discrete analogs of differential operators that satisfy the identities and theorems of vector and tensor calculus in discrete form the numerical approximation does not have spurious modes and mimics many fundamental properties of the underlying physical problem including conservation laws symmetries in the solution and the nondivergence of particular vector fields numerical examples demonstrate the high quality of the method when the medium is strongly discontinuous and for nonorthogonal nonsmooth computational grids