a signal processing approach to symmetry detection we present an algorithm that detects rotational and reflectional symmetries of two dimensional objects both symmetry types are effectively detected and analyzed using the angular correlation ac which measures the correlation between images in the angular direction the ac is accurately computed using the pseudopolar fourier transform which rapidly computes the fourier transform of an image on a near polar grid we prove that the ac of symmetric images is a periodic signal whose frequency is related to the order of the symmetry this frequency is recovered via spectrum estimation which is a proven technique in signal processing with a variety of efficient solutions we also provide a novel approach for finding the center of symmetry and demonstrate the applicability of our scheme to the analysis of real images

frequency and spatial sampling strategies for crosshole seismic waveform spectral inversion experiments analyses of synthetic frequency domain acoustic waveform data provide new insights into the design and imaging capability of crosshole surveys the full complex fourier spectral data offer significantly more information than other data representations such as the amplitude phase or hartley spectrum extensive eigenvalue analyses are used for further inspection of the information content offered by the seismic data the goodness of different experimental configurations is investigated by varying the choice of 1 the frequencies 2 the source and receiver spacings along the boreholes and 3 the borehole separation with only a few carefully chosen frequencies a similar amount of information can be extracted from the seismic data as can be extracted with a much larger suite of equally spaced frequencies optimized data sets should include at least one very low frequencycomponent the remaining frequencies should be chosen fromthe upper end of the spectrum available this strategy proved to be applica

3d finite difference modeling of elastic wave propagation in the laplace fourier domain abstractwith the recent interest in the laplace fourier domain full waveform inversion we have developed new heterogeneous 3d fourth and second order staggered grid finite difference schemes for modeling seismic wave propagation in the laplace fourier domain our approach is based on the integro interpolation technique for the velocity stress formulation in the cartesian coordinate system five averaging elastic coefficients and three averaging densities are necessary to describe the heterogeneous medium with harmonic averaging of the bulk and shear moduli and arithmetic averaging of density in the fourth order approximation we improved the accuracy of the scheme using a combination of integral identities for two elementary volumes small and large around spatial gridpoints where the wave variables are defined two solution approaches are provided both of which are solved with iterative krylov methods in the first approach the stress variables are eliminated and a linear system for the veloc

a mass redistributed finite element method mr fem for acoustic problems using triangular mesh the accuracy of numerical results using standard finite element method fem in acoustic problems will deteriorate with increasing frequency due to the dispersion error such dispersion error depends on the balance between the stiffness and mass of discretization equation systems this paper reports an improved finite element method fem for solving acoustic problems by re distributing the mass in the mass matrix to tune the balance aiming to minimize the dispersion errors this is done by shifting the integration point locations when computing the entries of the mass matrix while ensuring the mass conservation the new method is verified through the detailed numerical error analysis and a strategy is also proposed for the best mass redistribution in terms of minimizing dispersion error the relative dispersion error of present mass redistributed finite element method mr fem is found to be much smaller than the fem solution in both theoretical prediction and numerical examination the present mr fem works well by using the linear triangular elements that can be generated automatically which enables automation in computation and saving computational cost in mesh generation numerical examples demonstrate the advantages of mr fem in comparison with the standard fem using the same triangular meshes and quadrilateral meshes

accurate computation of zernike moments in polar coordinates an algorithm for high precision numerical computation of zernike moments is presented the algorithm based on the introduced polar pixel tiling scheme does not exhibit the geometric error and numerical integration error which are inherent in conventional methods based on cartesian coordinates this yields a dramatic improvement of the zernike moments accuracy in terms of their reconstruction and invariance properties the introduced image tiling requires an interpolation algorithm which turns out to be of the second order importance compared to the discretization error various comparisons are made between the accuracy of the proposed method and that of commonly used techniques the results reveal the great advantage of our approach

simultaneous multifrequency inversion of full waveform seismic data we present a simultaneous multifrequency inversion approach for seismic data interpretation this algorithm inverts all frequency data components simultaneously a data weighting scheme balances the contributions from different frequency data components so the inversion process does not become dominated by high frequency data components which produce a velocity image with many artifacts a gauss newton minimization approach achieves a high convergence rate and an accurate reconstructed velocity image by introducing a modified adjoint formulation we can calculate the jacobian matrix efficiently allowing the material properties in the perfectly matched layers pmls to be updated automatically during the inversion process this feature ensures the correct behavior of the inversion and implies that the algorithm is appropriate for realistic applications where a priori information of the background medium is unavailable two different regularization schemes an l2 norm and a weighted l2 norm function a

a dispersion minimizing compact finite difference scheme for the 2d helmholtz equation in this paper we present a dispersion minimizing compact finite difference scheme for solving the 2d helmholtz equation which is a fourth order scheme the error between the numerical wavenumber and the exact wavenumber is analyzed and a refined choice strategy based on minimizing the numerical dispersion is proposed for choosing weight parameters numerical results are presented to demonstrate the efficiency and accuracy of the compact finite difference scheme with refined parameters

a consistent test for the functional form of a regression based on a difference of variance estimators in this paper we study the problem of testing the functional form of a given regression model a consistent test is proposed which is based on the difference of the least squares variance estimator in the assumed regression model and a nonparametric variance estimator the corresponding test statistic can be shown to be asymptotically normal under the null hypothesis and under fixed alternatives with different rates of convergence corresponding to both cases this provides a simple asymptotic test where the asymptotic results can also be used for the calculation of the type ii error of the procedure at any particular point of the alternative and for the construction of tests for precise hypotheses finally the finite sample performance of the new test is investigated in a detailed simulation study which also contains a comparison with the commonly used tests

on difference based variance estimation in nonparametric regression when the covariate is high dimensional summary we consider the problem of estimating the noise variance in homoscedastic nonparametric regression models for low dimensional covariates t d d 1 2 difference based estimators have been investigated in a series of papers for a given length of such an estimator difference schemes which minimize the asymptotic mean squared error can be computed for d 1 and d 2 however from numerical studies it is known that for finite sample sizes the performance of these estimators may be deficient owing to a large finite sample bias we provide theoretical support for these findings in particular we show that with increasing dimension d this becomes more drastic if d4 these estimators even fail to be consistent a different class of estimators is discussed which allow better control of the bias and remain consistent when d4 these estimators are compared numerically with kernel type estimators which are asymptotically efficient and some guidance is given about when their use becomes necessary

high accuracy wave simulation revised derivation numerical analysis and testing of a nearly analytic integration discrete method for solving acoustic wave equation the nearly analytic integration discrete naid method for solving the two dimensional acoustic wave equation has been fully mathematically revised analyzed and tested the naid method is an alternative numerical modeling method for generating synthetic seismograms the acoustic wave equation is first transformed into a system of first order ordinary differential equations odes with respect to time variable t and then directly integrated at a small time interval of tn tn 1 to obtain semi discrete ordinary differential equations the integral kernel is expanded into a truncated taylor series to which the integration operator is explicitly applied the high order temporal derivatives involved in the integral kernel are replaced by high order spatial derivatives which then are approximately calculated as a weighted linear combination of the displacement the particle velocity and their spatial gradients in this article we investigate the theoretical properties of the revised naid method including the discrete error and the stability criteria numerical results for constant and layered velocity models show that comparing to the lax wendroff correction lwc scheme and the staggered grid finite difference method the naid method can effectively suppress the numerical dispersion and source noises caused by the discretization of the acoustic wave equation with too coarse spatial grids or when models have strong velocity contrasts between adjacent layers the proposed naid method has been applied in computing the acoustic wavefields for two heterogeneous models the corner edge model and the marmousi model promising numerical results illustrate that the naid method provides an encouraging tool for large scale and complex wave simulation and inversion problems based on the acoustic equation

comparison of source independent methods of elastic waveform inversion in this paper we investigate several source independent methods of nonlinear full waveform inversion of multicomponent elastic wave data this includes iterative estimation of source signature ies standard trace normalization stn and average trace normalization atn inversion methods all are based on the finite element method in the frequency domain one synthetic elastic crosshole model is used to compare the recovered images with all these methods as well as the known source signature kss inversion method the numerical experiments show that the ies method is superior to both stn and atn methods in two component elastic wave inversion in the frequency domain when the source signature is unknown the stn and atn methods have limitations associated with near zero amplitudes or polarity reversals in traces from one of the components which destroy the energy balance in the normalized traces and cause a loss of frequency information but the atn method is somewhat superior to the stn method in suppressing random noise and improving stability as the developed formulas and the numerical experiments show we suggest the ies method as a practical procedure for multicomponent seismic inversion

a bayesian approach to estimate uncertainty for full waveform inversion using a priori information from depth migration abstractfull waveform inversion fwi enables us to obtain high resolution subsurface images however estimating model uncertainties associated with this technique is still a challenging problem we have used a bayesian inference framework to estimate model uncertainties associated with fwi the uncertainties were assessed based on an a posteriori covariance operator evaluated at the maximum a posteriori model for the prior distribution we have used a spatially nonstationary covariance operator based on a plane wave construction with local dips measured from migrated images preconditioned frequency domain fwi was used to estimate the maximum a posteriori model efficient manipulation of the posterior covariance was based on a low rank approximation of the data misfit hessian preconditioned by the prior covariance operator the strong decay of the singular values indicated that data were mostly informative about a low dimensional subspace of model parameters to reduce computational cost of the random

projection based polygonality measurement measuring the degree to which a shape resembles a polygon referred to as polygonality is a difficult problem due to the intrinsic diversity in the form and distortion of shapes caused by digitization and similarity transformation this paper proposes a generic approach for this problem by performing the measurement in the projection space where the radon image of some primitive shapes which compose the shape becomes apparent the obtained measures take value in 0 1 with 1 corresponding to analytical polygons they are robust to additive noise and boundary distortion and is invariant to similarity transformation the new framework generalizes existing polygonal measures such as triangularity and quadrangularity in addition the ability to estimate polygon s geometric quantities in the projection space allows approximating a shape by analytical polygons the efficiency of the proposed approach is demonstrated through a number of experiments on both synthetic and real data sets

source independent time domain waveform inversion using convolved wavefields application to the encoded multisource waveform inversion full waveform inversion requires a good estimation of the source wavelet to improve our chances of a successful inversion this is especially true for an encoded multisource time domain implementation which conventionally requires separate source modeling as well as the fourier transform of wavefields as an alternative we have developed the source independent time domain waveform inversion using convolved wavefields specifically the misfit function consists of the convolution of the observed wavefields with a reference trace from the modeled wavefield plus the convolution of the modeled wavefields with a reference trace from the observed wavefield in this case the source wavelet of the observed and the modeled wavefields are equally convolved with both terms in the misfit function and thus the effects of the source wavelets are eliminated furthermore because the modeled wavefields play a role of low pass filtering the observed wavefields in the misfit function the frequency selection stra

applying finite element analysis to the memory variable formulation of wave propagation in anelastic media finite element methods are applied to solution of seismic wave motion in linear viscoelastic media using the memory variable formalism the displacements are represented as a superposition of a set of basis functions it is shown that if memory variables are represented using the spatial derivatives of those basis functions rather than the basis functions themselves the equations to be solved are simpler and require less computer memory using this formulation results for sh waves in one and two dimensions are calculated using a simple explicit finite element in space finite difference in time scheme these results agree with those found with a method of lines solution results in a homogeneous medium also agree with the frequency domain solutions of kjartansson s constant q method

a frequency space 2 d scalar wave extrapolator using extended 25 point finite difference operator finite difference frequency domain modeling for the generation of synthetic seismograms and crosshole tomography has been an active field of research since the 1980s the generation of synthetic seismograms with the time domain finite difference technique has achieved considerable success for waveform crosshole tomography and for wider applications in seismic reverse time migration this became possible with the rapid development of high performance computers however the space frequency x finite difference modeling technique is still beyond the capability of the modern supercomputer in terms of both cost and computer memory therefore finite difference time domain modeling is much more popular among exploration geophysicists a limitation of the space frequency domain is that the recently developed nine point scheme still requires that g the number of grid points per wavelength be 5 this value is greater than for most other numerical modeling techniques for example the pseudospectral scheme

mimetic seismic wave modeling including topography on deformed staggered grids abstractfinite difference methods for modeling seismic waves are known to be inaccurate when including a realistic topography due to the large dispersion errors that appear in the modelled surface waves and the scattering introduced by the staircase approximation to the topography as a consequence alternatives to finite difference methods have been proposed to circumvent these issues we present a new numerical scheme for 3d elastic wave propagation in the presence of strong topography this finite difference scheme is based upon a staggered grid of the lebedev type or fully staggered grid fsg it uses a grid deformation strategy to make a regular cartesian grid conform to a topographic surface in addition the scheme uses a mimetic approach to accurately solve the free surface condition and hence allows for a less restrictive grid spacing criterion in the computations the scheme can use high order operators for the spatial derivatives and obtain low dispersion results with as few as six points pe

grid dispersion and stability criteria of some common finite element methods for acoustic and elastic wave equations purely numerical methods based on finite element approximation of the acoustic or elastic wave equation are becoming increasingly popular for the generation of synthetic seismograms we present formulas for the grid dispersion and stability criteria for some popular finite element methods fem for wave propagation namely classical and spectral fem we develop an approach based on a generalized eigenvalue formulation to analyze the dispersive behavior of these fems for acoustic or elastic wave propagation that overcomes difficulties caused by irregular node spacing within the element and the use of high order polynomials as is the case for spectral fem analysis reveals that for spectral fem of order four or greater dispersion is less than 0 2 at four to five nodes per wavelength and dispersion is not angle dependent new results can be compared with grid dispersion results of some classical finite difference methods fdm used for acoustic or elastic wave propagation analysis reveals that fdm and classical fem require a larger sampling ratio than a spectral fem to obtain results with the same degree of accuracy the staggered grid fdm is an efficient scheme but the dispersion is angle dependent with larger values along the grid axes on the other hand spectral fem of order four or greater is isotropic with small dispersion making it attractive for simulations with long propagation times

accurate and efficient seismic modeling in random media the optimum method for seismic modeling in random media must 1 be highly accurate to be sensitive to subtle effects of wave propagation 2 allow coarse sampling to model media that are large compared to the scale lengths and wave propagation distances which are long compared to the wavelengths this is necessary to obtain statistically meaningful overall attributes of wavefields high order staggered grid finite difference algorithms and the pseudospectral method combine high accuracy in time and space with coarse sampling investigations for random media reveal that both methods lead to nearly identical wavefields the small differences can be attributed mainly to differences in the numerical dispersion this result is important because it shows that errors of the numerical differentiation which are caused by poor polynomial interpolation near discontinuities do not accumulate but cancel in a random medium where discontinuities are numerous the differentiator can be longer than the medium scale leng

robust source independent elastic full waveform inversion in the time domain abstractaccurate estimation of source wavelet is crucial in a successful full waveform inversion fwi however it cannot be guaranteed in the case of real seismic data we have developed time domain source independent elastic fwi using the convolution based objective function that was originally developed for acoustic fwi we have applied a new time window on the reference traces in the objective function to suppress the noises induced by the convolution and crosscorrelation operations also we have adopted l1 huber and hybrid norm objective functions to improve the antinoise ability of our algorithm we implemented a multiscale inversion strategy to conduct the tests with a quasi newton limited memory broyden fletcher goldfarb shanno l bfgs method to reduce the sensitivity to initial models and to improve the quality of inversion results synthetic tests verified that the new added time window can not only improve the inversion results but also accelerate the convergence rate our method can be

a central limit theorem for generalized quadratic forms random variables of the form w n mathop sum limits 1 leqq i leqq n mathop sum limits text 1 leqq j leqq n w ijn x i x j are considered with x i independent not necessarily identically distributed and w ijn borel functions such that w ijn x i x j is square integrable and has vanishing conditional expectations e w ijn x i x j x i e w ijn x i x j x j 0 text a text s text a central limit theorem is proved under the condition that the normed fourth moment tends to 3 under some restrictions the condition is also necessary finally conditions on the individual tails of w ijn x i x j and an eigenvalue condition are given that ensure asymptotic normality of w n

source independent full waveform inversion of seismic data a set of seismic trace data is collected in an input data set that is first fourier transformed in its entirety into the frequency domain a normalized wavefield is obtained for each trace of the input data set in the frequency domain normalization is done with respect to the frequency response of a reference trace selected from the set of seismic trace data the normalized wavefield is source independent complex and dimensionless the normalized wavefield is shown to be uniquely defined as the normalized impulse response provided that a certain condition is met for the source this property allows construction of the inversion algorithm disclosed herein without any source or source coupling information the algorithm minimizes the error between data normalized wavefield and the model normalized wavefield the methodology is applicable to any 3 d seismic problem and damping may be easily included in the process

frequency domain classification of cyclic and dihedral symmetries of finite 2 d patterns this paper presents a very simple and effective algorithm for classifying cyclic and dihedral symmetries from images of finite two dimensional patterns it consists in an extension of a frequency domain algorithm proposed by the author for estimating two dimensional rotation of rigid objects within this framework the estimation of rotation is conveniently turned into the problem of detecting two orthogonal lines within the locus of the zero crossings of a certain function represented in orthogonal cartesian coordinates the main contribution of this paper is the use and interpretation of the ambiguities that arise in such a locus in the form of additional pairs of orthogonal lines when applied to symmetric patterns and the development of other algorithmic steps which allow a simple and fast discrimination between cyclic and dihedral symmetries unlike any other available method this algorithm does not require any conversion from cartesian to polar image representations besides classification reliability and consistency observed through several dozen testing experiments the nice features of this new method also include robustness to noise and ease of implementation with fast fourier transform algorithms which makes it amenable to almost real time pattern classifications applications several examples representative of the performance of the algorithm through extensive testing are reported and discussed in the paper

global superheating field for superconductors in a large bounded interval abstract in continuation of our study of the ginzburg landau equations we study the asymptotic behavior of the superheating field for a film of width 2d in the regime small d large where is the ginzburg landau parameter using new estimates on the inner magnetic potential we improve our conditions on the parameter d and we prove that the local superheating field analyzed in our preceding paper c bolley f foucher b helffer superheating field for the ginzburg landau equations in the case of a large bounded interval j math phys 41 11 2000 7263 7289 is a global one in the regime and d large we also give accurate estimates on the superheating field

invariant image recognition by zernike moments the problem of rotation scale and translation invariant recognition of images is discussed a set of rotation invariant features are introduced they are the magnitudes of a set of orthogonal complex moments of the image known as zernike moments scale and translation invariance are obtained by first normalizing the image with respect to these parameters using its regular geometrical moments a systematic reconstruction based method for deciding the highest order zernike moments required in a classification problem is developed the quality of the reconstructed image is examined through its comparison to the original one the orthogonality property of the zernike moments which simplifies the process of image reconstruction make the suggest feature selection approach practical features of each order can also be weighted according to their contribution to the reconstruction process the superiority of zernike moment features over regular moments and moment invariants was experimentally verified

source receiver compression scheme for full waveform seismic inversion we have developed a source receiver compression approach for reducing the computational time and memory usage of the acoustic and elastic full waveform inversions by detecting and quantifying the extent of redundancy in the data we assembled a reduced set of simultaneous sources and receivers that are weighted sums of the physical sources and receivers used in the survey because the numbers of these simultaneous sources and receivers could be significantly less than those of the physical sources and receivers the computational time and memory usage of any gradient type inversion method such as steepest descent nonlinear conjugate gradient contrast source inversion and quasi newton methods could be reduced the scheme is based on decomposing the data into their principal components using a singular value decomposition approach and the data reduction is done through the elimination of the small eigenvalues consequently this would suppress the effect of noise in the data moreover taking advantage

full wave directional illumination analysis in the frequency domain directional illumination analysis based on one way wave equations has been studied extensively however its inherent limitations e g one way propagation wide angle error and amplitude inaccuracy can severely hinder its applications for accurate survey design and true reflection imaging corrections in complex media we have analyzed the illumination in the frequency domain using full two way wave propagators considering the extensive computation and huge storage required for time domain methods and the fact that the illumination is frequency dependent this full wave analysis can provide frequency dependent full angle true amplitude illumination not only for the downgoing waves but also for the upgoing waves including turning waves and reflected waves two methods were considered to decompose the full wavefield into the local angle domain a direct full dimensional decomposition and more efficient split step decomposition composed of three lower dimensional decompositions the results of illumination analysis demonstrated the advantages of this method the two decomposition methods produced similar results

on symmetry detection a straight line is an axis ofsymmetry of a planar figure if the figure is invariant to reflection with respect to that line the purpose of this correspondence is to describe an o n log n time algorithm for enumerating all the axes of symmetry of a planar figure which is made up of possibly intersecting segments circles points etc the solution involves a reduction of the problem to a combinatorial question on words our algorithm is optimal since we can establish an n log n time lower bound for this problem

intrinsic manova for riemannian manifolds with an application to kendall s space of planar shapes we propose an intrinsic multifactorial model for data on riemannian manifolds that typically occur in the statistical analysis of shape due to the lack of a linear structure linear models cannot be defined in general to date only one way manova is available for a general multifactorial model we assume that variation not explained by the model is concentrated near elements defining the effects by determining the asymptotic distributions of respective sample covariances under parallel transport we show that they can be compared by standard manova often in applications manifolds are only implicitly given as quotients where the bottom space parallel transport can be expressed through a differential equation for kendall s space of planar shapes we provide an explicit solution we illustrate our method by an intrinsic two way manova for a set of leaf shapes while biologists can identify genotype effects by sight we can detect height effects that are otherwise not identifiable

validation of linear regression models a new test is proposed in order to verify that a regression function say g has a prescribed linear parametric form this procedure is based on the large sample behavior of an empirical l 2 distance between g and the subspace u spanned by the regression functions to be verified the asymptotic distribution of the test statistic is shown to be normal with parameters depending only on the variance of the observations and the l 2 distance between the regression function g and the model space u based on this result a test is proposed for the hypothesis that g is not in a preassigned l 2 neighborhood of u whichallows the verification of the model u at a controlled type i error rate the suggested procedure is very easy to apply because of its asymptotic normal law and the simple form of the test statistic in particular it does not require nonparametric estimators of the regression function and hence the test does not depend on the subjective choice of smoothing parameters

a recursive green s function technique for acoustic scattering from heterogeneous objects a fast efficient algorithm for computing acoustic fields scattered by inhomogeneous objects in an otherwise homogeneous space is presented the algorithm called the recursive green s function method rgfm constructs the domain green s function by recursively combining known green s functions from smaller subdomains the fields on the scatterer surface are then computed using a boundary integral formulation proper implementation of the rgfm results in a storage requirement of o n and computational costs of o n3 2 and o n2 for two and three dimensional problems respectively where n is the total number of discrete points in the inhomogeneous region results are compared with those obtained from exact solutions to show the accuracy of the method

a new 3 d numerical approach to the solution of pdes with optimal accuracy on irregular domains and cartesian meshes abstract a new 3 d numerical approach for the time dependent wave and heat equations as well as for the time independent laplace equation on irregular domains with the dirichlet boundary conditions has been developed trivial cartesian meshes and simple 27 point uniform and nonuniform stencil equations are used for 3 d irregular domains the calculation of the coefficients of the stencil equations is based on the minimization of the local truncation error of the stencil equations and yields the optimal order of accuracy of the new technique very small distances 0 1 h 1 0 9 h where h is the grid size between the grid points of a cartesian mesh and the boundary do not worsen the accuracy of the new technique at similar 27 point stencils the accuracy of the new approach is two orders higher than that for the linear finite elements the numerical results for irregular domains show that at the same number of degrees of freedom the new approach is even much more accurate than the high order up to the fifth order tetrahedral finite elements with much wider stencils the wave and heat equations can be uniformly treated with the new approach the order of the time derivative in these equations does not affect the coefficients of the stencil equations of the semi discrete systems the new approach can be directly applied to other partial differential equations

2d frequency domain elastic full waveform inversion using time domain modeling and a multistep length gradient approach abstractto decouple the parameters in elastic full waveform inversion fwi we evaluated a new multistep length gradient approach to assign individual weights separately for each parameter gradient and search for an optimal step length along the composite gradient direction to perform wavefield extrapolations for the inversion we used parallelized high precision finite element fe modeling in the time domain the inversion was implemented in the frequency domain the data were obtained at every subsurface grid point using the discrete fourier transform at each time domain extrapolation step we also used frequency selection to reduce cycle skipping time windowing to remove the artifacts associated with different source spatial patterns between the test and predicted data and source wavelet estimation at the receivers over the full frequency spectrum by using a fast fourier transform in the inversion the velocity and density reconstructions behaved differently as a low wavenumber tomography for v

accurate calculation of zernike moments zernike moments zms are very effective global image descriptors which are used in many digital image processing applications the digitization process compromises the accuracy of the moments and therefore several of its properties are affected there are two major discretization errors namely the geometric error and numerical integration error in this paper we propose two new algorithms which eliminate these errors the first algorithm performs the exact computation of geometric moments gms over a unit disk and then uses gms to zms relationship to compute the latter this algorithm is computationally more expensive and it becomes numerically instable for higher order moments therefore we develop a second algorithm based on gaussian quadrature numerical integration the second algorithm reduces both the errors simultaneously and its accuracy increases as the degree of gaussian quadrature numerical integration increases the proposed algorithms are observed to provide very accurate zms which result in improved image reconstruction reduction in reconstruction error and improvement in rotation and scale invariance exhaustive experiments are provided to support improved accuracy of zms and time complexity analysis is performed for the existing and the proposed methods

statistical shape analysis clustering learning and testing using a differential geometric treatment of planar shapes we present tools for 1 hierarchical clustering of imaged objects according to the shapes of their boundaries 2 learning of probability models for clusters of shapes and 3 testing of newly observed shapes under competing probability models clustering at any level of hierarchy is performed using a minimum variance type criterion and a markov process statistical means of clusters provide shapes to be clustered at the next higher level thus building a hierarchy of shapes using finite dimensional approximations of spaces tangent to the shape space at sample means we implicitly impose probability models on the shape space and results are illustrated via random sampling and classification hypothesis testing together hierarchical clustering and hypothesis testing provide an efficient framework for shape retrieval examples are presented using shapes and images from eth surrey and amcom databases

on the recovery of a function on a circular domain we consider the problem of estimating a function f x y on the unit disk f x y x sup 2 y sup 2 spl les 1 given a discrete and noisy data recorded on a regular square grid an estimate of f x y based on a class of orthogonal and complete functions over the unit disk is proposed this class of functions has a distinctive property of being invariant to rotation of axes about the origin of coordinates yielding therefore a rotationally invariant estimate for radial functions the orthogonal set has a particularly simple form being related to the classical legendre polynomials we give the statistical accuracy analysis of the proposed estimate of f x y in the sense of the l sub 2 metric it is found that there is an inherent limitation in the precision of the estimate due to the geometric nature of a circular domain this is explained by relating the accuracy issue to the celebrated problem in the analytic number theory called the lattice points of a circle in fact the obtained bounds for the mean integrated squared error are determined by the best known result so far on the problem of lattice points within the circular domain

a consistent nonparametric test of symmetry in linear regression models abstract this article proposes a consistent nonparametric test of the hypothesis that the disturbance in a linear regression model is distributed symmetrically around zero simulation results show that the test has good size and power properties for sample sizes as small as 50 we illustrate the use of the test in a cross country model of inflation and monetary growth

dispersion minimizing quadrature rules for c1 quadratic isogeometric analysis abstract we develop quadrature rules for the isogeometric analysis of wave propagation and structural vibrations that minimize the discrete dispersion error of the approximation the rules are optimal in the sense that they only require two quadrature points per element to minimize the dispersion error barto et al 2017 1 and they are equivalent to the optimized blending rules we recently described our approach further simplifies the numerical integration instead of blending two three point standard quadrature rules we construct directly a single two point quadrature rule that reduces the dispersion error to the same order for uniform meshes with periodic boundary conditions also we present a 2 5 point rule for both uniform and non uniform meshes with arbitrary boundary conditions consequently we reduce the computational cost by using the proposed quadrature rules various numerical examples demonstrate the performance of these quadrature rules

robust rotation angle estimator the conventional method of estimating the rotation angle of a pattern using the principal axes is not suitable for circular symmetric patterns since their eigenvalues are similar in both directions in the paper a robust method of estimating a rotation angle using the phase information of zernike moments is presented the experimental results show that the proposed method estimates the rotation angle of the circular symmetric patterns more accurately than the principal axes method even in the presence of noise

accuracy of finite difference and finite element modeling of the scalar and elastic wave equations numerical solutions of the scalar and elastic wave equations have greatly aided geophysicists in both forward modeling and migration of seismic wave fields in complicated geologic media and they promise to be invaluable in solving the full inverse problem this paper quantitatively compares finite difference and finite element solutions of the scalar and elastic hyperbolic wave equations for the most popular implicit and explicit time domain and frequency domain techniques in addition to versatility and ease of implementation it is imperative that one choose the most cost effective solution technique for a fixed degree of accuracy to be of value a solution technique must be able to minimize 1 numerical attenuation or amplification 2 polarization errors 3 numerical anisotropy 4 errors in phase and group velocities 5 extraneous numerical parasitic modes 6 numerical diffraction and scattering and 7 errors and transmission coefficients this paper shows that in homogeneous media the

on the accuracy of zernike moments for image analysis we give a detailed analysis of the accuracy of zernike moments in terms of their discretization errors and the reconstruction power it is found that there is an inherent limitation in the precision of computing the zernike moments due to the geometric nature of a circular domain this is explained by relating the accuracy issue to a celebrated problem in analytic number theory of evaluating the lattice points within a circle

comparing nonparametric versus parametric regression fits in general there will be visible differences between a parametric and a nonparametric curve estimate it is therefore quite natural to compare these in order to decide whether the parametric model could be justified an asymptotic quantification is the distribution of the integrated squared difference between these curves we show that the standard way of bootstrapping this statistic fails we use and analyse a different form of bootstrapping for this task we call this method the wild bootstrap and apply it to fitting engel curves in expenditure data analysis

iterative finite difference solution analysis of acoustic wave equation in the laplace fourier domain we have investigated numerical characteristics of iterative solutions to the acoustic wave equation in the laplacefourier lf domain we transformed the time domain acoustic wave equation into the lf domain the transformed equation was discretized with finite differences and was solved with iterative methods finite difference modeling experiments demonstrate that iterative methods require an infinitesimal stopping tolerance to accurately compute the pressure field especially at long offsets to understand the requirement for such infinitesimal tolerance values we analyzed the evolution of intermediate solution vectors residual vectors and search direction vectors during the iteration the analysis showed that the requirement arises from the fact that in the solution space the amplitude of the pressure field varies more than sixty orders of magnitude on the common log scale accordingly we propose a rule of thumb for choosing a proper stopping tolerance value we also examined numerical dispersion errors in terms of the grid sampling resolutions per skin depth and wavelength we found that despite the similarity of the form of the acoustic wave and electromagnetic diffusion equations the former is different from the latter due to the fact that in the lf domain the skin depth of the acoustic wave equation is decoupled from its wavelength this aspect requires that in the lf domain its grid size be determined by considering the minimum grid sampling resolutions based not only the wavelength but also the skin depth

dispersion minimized mass for isogeometric analysis abstract we introduce the dispersion minimized mass for isogeometric analysis to approximate the structural vibration which we model as a second order differential eigenvalue problem the dispersion minimized mass reduces the eigenvalue error significantly from the optimum order of 2 p to the superconvergence order of 2 p 2 for the p th order isogeometric elements with maximum continuity which in return leads to a more accurate method we first establish the dispersion error where the leading error term is explicitly written in terms of the stiffness and mass entries for arbitrary polynomial order isogeometric elements we derive the dispersion minimized mass in one dimension by solving a p dimensional local matrix problem for the p th order approximation and then extend it to multiple dimensions on tensor product grids we show that the dispersion minimized mass can also be obtained by approximating the mass matrix using optimally blended quadratures we generalize the lower order quadrature blending results to arbitrary polynomial order isogeometric approximations as well as to arbitrary quadrature rules various numerical examples validate the eigenvalue and eigenfunction error estimates

waveform inversion using a logarithmic wavefield although waveform inversion has been studied extensively since its beginning 20 years ago applications to seismic field data have been limited and most of those applications have been for global seismology or engineering seismology scale problems not for exploration scale data as an alternative to classical waveform inversion we propose the use of a new objective function constructed by taking the logarithm of wavefields allowing consideration of three types of objective function namely amplitude only phase only or both in our waveform inversion we estimate the source signature as well as the velocity structure by including functions of amplitudes and phases of the source signature in the objective function we compute the steepest descent directions by using a matrix formalism derived from a frequency domain finite element finite difference modeling technique our numerical algorithms are similar to those of reverse time migration and waveform inversion based on the adjoint state of the wave equation in order to demonstrate the practical applicability of our algorithm we use a synthetic data set from the marmousi model and seismic data collected from the korean continental shelf for noisefree synthetic data the velocity structure produced by our inversion algorithm is closer to the true velocity structure than that obtained with conventional waveform inversion when random noise is added the inverted velocity model is also close to the true marmousi model but when frequencies below 5 hz are removed from the data the velocity structure is not as good as those for the noise free and noisy data for field data we compare the time domain synthetic seismograms generated for the velocity model inverted by our algorithm with real seismograms and find that the results show that our inversion algorithm reveals short period features of the subsurface although we use wrapped phases in our examples we still obtain reasonable results we expect that if we were to use correctly unwrapped phases in the inversion algorithm we would obtain better results

time space domain mesh free finite difference based on least squares for 2d acoustic wave modeling abstractfinite difference fd methods approximate derivatives through a weighted summation of function values from neighboring nodes traditionally these neighboring nodes are assumed to be distributed regularly such as in square or rectangular lattices to improve geometric flexibility one option is to develop fd in a mesh free discretization in which scattered nodes can be placed suitably with respect to irregular boundaries or arbitrarily shaped anomalies without coordinate transformation or forming any triangles or tetrahedra etc these mesh free fds have had successful applications especially in computational geoscience however they are all space domain fd schemes in which fd coefficients are derived by approximating spatial derivatives individually in the space domain for acoustic wave modeling it has been proven that space domain fd methods normally have higher dispersion error than time space domain fd methods which determine fd coefficients by approximating the time space domain disp

efficient laplace domain modeling and inversion using an axis transformation technique abstractwe tested an axis transformation technique for modeling wave propagation in the laplace domain using a finite difference method this technique enables us to use small grids near the surface and large grids at depth accordingly we can reduce the number of grids and attain computational efficiency in modeling and inversion in the laplace domain we used a dispersion analysis and comparisons between modeled wavefields obtained on the regular and transformed axes we demonstrated in a synthetic laplace domain inversion technique shows that this method is efficient and yields a result comparable to that of a laplace domain inversion using a regular grid in a synthetic inversion example the memory usage reduced to less than 33 and the computation time reduced to 39 of those required for the regular grid case using a logarithmic transformation function

velocity model building by 3d frequency domain full waveform inversion of wide aperture seismic data we assessed 3d frequency domain fd acoustic fullwaveforminversionfwidataasatooltodevelophigh resolution velocity models from low frequency global offset data the inverse problem was posed as a classic leastsquares optimization problem solved with a steepest descent method inversion was applied to a few discrete frequencies allowing management of a limited subset of the 3d data volume the forward problem was solved with a finite difference frequency domain method based on a massively paralleldirectsolver allowingefficientmultiple shotsimulations the inversion code was fully parallelized for distributedmemory platforms taking advantage of a domain decomposition of the modeled wavefields performed by the direct solver after validation on simple synthetic tests fwi was applied to two targets channel and thrust system of the 3d seg eageoverthrustmodel correspondingto3ddomains of 78 752 25 km and 13 513 54 65 km respectively the maximum inverted frequencies are 15 and 7 hz for the two applications amaximum of 30 dual core biprocessor nodes with 8 gb of shared memory per node were used for the second target the main structures were imaged successfully at a resolution scale consistent with the inverted frequencies our study confirms the feasibility of 3d frequency domain fwi of global offset data on large distributed memory platforms to develop high resolution velocity models these high velocity models may provide accurate macromodelsforwave equationprestackdepthmigration

perturbational and nonperturbational inversion of rayleigh wave velocities the inversion of rayleigh wave dispersion curves is a classic geophysical inverse problem we have developed a set of matlab codes that performs forward modeling and inversion of rayleigh wave phase or group velocity measurements we describe two different methods of inversion a perturbational method based on finite elements and a nonperturbational method based on the recently developed dix type relation for rayleigh waves in practice the nonperturbational method can be used to provide a good starting model that can be iteratively improved with the perturbational method although the perturbational method is well known we solve the forward problem using an eigenvalue eigenvector solver instead of the conventional approach of root finding features of the codes include the ability to handle any mix of phase or group velocity measurements combinations of modes of any order the presence of a surface water layer computation of partial derivatives due to changes in material properties and layer boundaries and the implementation of an automatic grid of layers that is optimally suited for the depth sensitivity of rayleigh waves

seismic waveform inversion in the frequency domain part 1 theory and verification in a physical scale model seismic waveforms contain much information that is ignored under standard processing schemes seismic waveform inversion seeks to use the full information content of the recorded wavefield in this paper i present apply and evaluate a frequency space domain approach to waveform inversion the method is a local descent algorithm that proceeds from a starting model to refine the model in order to reduce the waveform misfit between observed and model data the model data are computed using a full wave equation viscoacoustic frequency domain finite difference method ray asymptotics are avoided and higher order effects such as diffractions and multiple scattering are accounted for automatically the theory of frequency domain waveform wavefield inversion can be expressed compactly using a matrix formalism that uses finite difference finite element frequency domain modeling equations expressions for fast local descent inversion using back propagation techniques then follow naturally implementation of

wave propagation in heterogeneous media effects of fine scale heterogeneity we study the multiple scattering effects caused by fine scale heterogeneity for this purpose it is unreasonable to rely on a linearization of the dependency of the wavefield in the parameters that describe the medium therefore the only tools that correctly model wave propagation are based on the numerical solution of the two way wave equation by finite differences or finite elements a fine grid provides a straightforward approach to account for fine scale heterogeneity in this situation there is no need for high order schemes variational methods allow us to exhibit a numerical scheme that accounts for heterogeneity and that is by construction stable provided a stability condition is fulfilled this condition is a sufficient stability condition contrary to the classical necessary stability conditions in addition general mathematical results prove the finite difference solution is close to the solution of the wave equation when the grid is fine enough the multiple scattering effects caused by fine scale heterogeneity are very important in particular we observe that imaging the so computed synthetic data by standard migration techniques that assume a linearization of the above mentioned dependency shows a strongly noise corrupted image this illustrates the importance of preprocessing data to remove the effect of multiple scattering we try to improve the signal to noise ratio by removing multiples related to the free surface although significant noise reduction is achieved even more sophisticated preprocessing is required to obtain a clear image of the subsurface

testing for lack of fit in inverse regression with applications to biophotonic imaging summary we propose two test statistics for use in inverse regression problems y d k c where k is a given linear operator which cannot be continuously inverted thus only noisy indirect observationsy for the function are available both test statistics have a counterpart in classical hypothesis testing where they are called the order selection test and the data driven neyman smooth test we also introduce two model selection criteria which extend the classical akaike information criterion and bayes information criterion to inverse regression problems in a simulation study we show that the inverse order selection and neyman smooth tests outperform their direct counterparts in many cases the theory is motivated by data arising in confocal fluorescence microscopy here images are observed with blurring modelled as convolution and stochastic error at subsequent times the aim is then to reduce the signal to noise ratio by averaging over the distinct images in this context it is relevant to decide whether the images are still equal or have changed by outside influences such as moving of the object table

3d frequency domain finite difference viscoelastic wave modeling using weighted average 27 point operators with optimal coefficients abstractexperimental data suggest that the viscoelastic behavior of rocks is more easily and accurately described in the frequency domain than in the time domain supporting the idea of simulating seismic wave propagation in the frequency domain we evaluated weighted averaged 27 point finite difference operators for 3d viscoelastic wave modeling in the frequency domain within the proposed framework we developed general equations for normalized phase velocities that can be used with arbitrary finite difference operators three sets of weighting coefficients for second order central finite difference operators that minimize the numerical dispersion for up to five grid points per wavelength were found using a damped least squares ls criterion as well as a global optimization scheme based on l1 and l2 norm criteria the three sets produced very similar dispersion curves and improvement provided by global optimization appeared marginal in this respect we also evaluated a discrete form for the heterogen

on image analysis by moments research has been performed investigating the use of moments for pattern recognition in recent years the basic problem of the influence of discretization and noise on moment accuracy as object descriptors has been barely investigated in this paper the detailed error analysis involved in the moment method is discussed several new techniques to increase the accuracy and efficiency of moment descriptor are proposed we utilize these results for the problem of image reconstruction from the orthogonal legendre moments computed from discrete and noisy data the automatic selection of an optimal number of moments is also discussed

testing for image symmetries with application to confocal microscopy statistical tests are introduced for checking whether an image function f x y defined on the unit disc d x y x 2 y 2 les 1 is invariant under certain symmetry transformations of d given that discrete and noisy data are observed invariance under reflections or under rotations by rational angles is considered as well as rotational invariance these symmetry relations can be naturally expressed as restrictions for the zernike moments of f x y therefore the test statistics are based on the l 2 distance between zernike series estimates of the image function itself and its version obtained after applying the symmetry transformation the asymptotic distribution of the test statistics under both the hypothesis of symmetry as well as under fixed alternatives is derived furthermore the quality of the asymptotic approximations via simulation studies is investigated the usefulness of our theory is verified by examining an important problem in confocal microscopy i e possible imprecise alignments in the optical path of the microscope are investigated for optical systems with rotational symmetry the theoretical point spread function psf is reflection symmetric with respect to two orthogonal axes and rotationally invariant if the detector plane matches the optical plane of the microscope the tests are used to investigate whether the required symmetries can indeed be detected in the empirical psf

numerical analysis of a narrow angle one way elastic wave equation and extension to curvilinear coordinates in this paper we review the finite difference implementation of a narrow angle one way vector wave equation for elastic 3d media extrapolation is performed in the frequency domain where the second order accurate lateral spatial difference operators are sufficiently accurate for narrow angle propagation we perform a numerical analysis of the finite difference scheme to highlight the stability and dispersion characteristics the von neumann stability criterion indicates that extracting a reference phase during the extrapolation step noticeably improves the forward marching scheme and dispersion analysis shows that numerical grid anisotropy is minimal for the propagation path lengths source pulse spectral content and angular range of forward propagation of interest although the algorithm is reasonable its computational efficiency is limited by the second order accurate extrapolation step therefore the extrapolation scheme can be improved we extend the cartesian narrow angle formulation to curvilinear coordinates where the computational grid tracks the true wavefront in a reference medium and the wavefield derivative normal to the reference wavefront is evaluated locally using the cartesian propagator an example of curvilinear extrapolation for a simple model consisting of a high velocity sphere within a homogeneous background velocity structure shows that the narrow angle propagator is capable of modeling frequency dependent geometric spreading and diffraction effects in curvilinear coordinates

crosshole seismic inversion with normalized full waveform amplitude data we investigate a simple scheme for full waveform amplitude spectrum inversion of crosshole seismic data with an unknown source wavelet the method is based on our 2d 2 5d finite element method of acoustic wave modeling the normalized amplitude data defined as the spectral ratio of the original trace amplitude to the average amplitude for the entire common shot gather are used for full waveform inversion in the frequency domain in essence the normalization of amplitudes removes the source wavelet and is easily carried out in the time domain or frequency domain from crosshole seismic surveying two synthetic models simulating respectively a dipping fracture model and a complicated sedimentary structure are examined with the inversion scheme the numerical results show that clear images of the targets can be obtained by using the single and multiple frequency data for comparison two other amplitude inversion schemes in which the source signatures are known and estimated with an iterative procedure were carried out for the same models the results show that the iteratively source estimated procedure also produces satisfactory images of the velocity structure and yields an approximate amplitude of the source wavelet with the multifrequency data three such schemes yield very competitive results

dispersion reducing finite elements for transient acoustics this paper focuses on increasing the accuracy of low order four node quadrilateral finite elements for the transient analysis of wave propagation modified integration rules originally proposed for time harmonic problems provide the basis for the proposed technique the modified integration rules shift the integration points to locations away from the conventional gauss or gauss lobatto integration points with the goal of reducing the discretization errors specifically the dispersion error presented here is an extension of the idea to time dependent analysis using implicit as well as explicit time stepping schemes the locations of the stiffness integration points remain unchanged from those in time harmonic case on the other hand the locations of the integration points for the mass matrix depend on the time stepping scheme and the step size furthermore the central difference method needs to be modified from its conventional form to facilitate fully explicit computation the superior performance of the proposed algorithms is illustrated with the help of several numerical examples

on the detection of the axes of symmetry of symmetric and almost symmetric planar images a simple algorithm is presented for finding all the axes of symmetry of symmetric and almost symmetric planar images having nonuniform gray level intensity images the technique which can also be used in conjunction with planar curves is based on the identification of the centroids of the given image and other related sets of points followed by a maximization of a specially defined coefficient of symmetry owing to the nature of the required procedures which are strictly based on the computation of mean values the method is not particularly affected by the presence of noise and is especially suited for application in conjunction with digitized figures

optimal blended spectral element operators for acoustic wave modeling spectral element methods based on high order polynomials are among the most commonly used techniques for computing accurate simulations of wave propagation phenomena in complex media however to retain computational efficiency very high order polynomials cannot be used and errors such as numerical dispersion and numerical anisotropy cannot be totally avoided in the present work we devise an approach for reducing such errors by considering modified discrete wave operators we analyze consistent and lumped operators together with blended operators weighted averages of consistent and lumped operators furthermore using the operator blending approach and a novel dispersion analysis method we develop optimal spectral element operators that have increased numerical accuracy without resorting to very high order operators the new operators are faster and computationally more efficient than consistent operators our approach is based on the tensor product decomposition of the element matrices into 1d f

speed of estimation in positron emission tomography and related inverse problems we study a continuous idealization of the image reconstruction problem in positron emission tomography considered as an example of bivariate density estimation based on indirect observations given a large sample of indirect observations we consider the size of the equivalent sample of observations whose original exact positions would allow equally accurate estimation of the image of interest both for indirect and for direct observations we establish exact minimax rates of convergence of estimation for all possible estimators over suitable smoothness classes of functions the approach of the paper is applicable to a wide class of linear inverse problems

a comparative study of finite element and spectral element methods in seismic wavefield modeling the computational accuracy and efficiency of finite element method and spectral element method sem are investigated thoroughly in time domain elastic wavefield modeling the diagonal mass matrices of the fem and sem free from matrix inversion are compared comprehensively by making full use of the mass lumped technique and quadrature rules we investigate the fem and sem based respectively on quadrilateral with the polynomials of degrees one and two and on triangular grids with the polynomials of degrees one and three generally the numerical solutions based on quadrilateral grids have a higher precision than those computed on triangular grids when the same order of polynomials is used the fem has a comparable accuracy to the sem with the same number of interpolant points in view of the triangular and quadrilateral sems the former suffers from larger computational costs and relatively lower accuracy compared with the latter furthermore the convergence study proves that the triangular sem produces consistently larger errors than the quadrilateral sem for any order and element sizes however the triangular sem can adapt to arbitrary complex geometries effectively in terms of efficiency the fem has an efficiency comparable with the sem on condition that the order of interpolation polynomials is identical in addition a perfectly matched layer pml boundary condition in variational form is deduced by introducing four intermediate variables in frequency domain the pml avoids convolution calculation and obtains an exact solution through inverse fourier transform in time domain the numerical examples verify the validity and effectiveness of the pml

fast generic polar harmonic transforms generic polar harmonic transforms have recently been proposed to extract rotation invariant features from images and their usefulness has been demonstrated in a number of pattern recognition problems however direct computation of these transforms from their definition is inefficient and is usually slower than some efficient computation strategies that have been proposed for other methods this paper presents a number of novel computation strategies to compute these transforms rapidly the proposed methods are based on the inherent recurrence relations among complex exponential and trigonometric functions used in the definition of the radial and angular kernels of these transforms the employment of these relations leads to recursive and addition chain based strategies for fast computation of harmonic function based kernels experimental results show that the proposed method is about 10 times faster than direct computation and 5 times faster than fast computation of zernike moments using the q recursive strategy thus among all existing rotation invariant feature extraction methods polar harmonic transforms are the fastest

orthogonal moment features for use with parametric and non parametric classifiers this research examines a variety of approaches for using two dimensional orthogonal polynomials for the recognition of handwritten arabic numerals it also makes use of parametric and non parametric statistical and neural network classifiers polynomials including legendre zernike and pseudo zernike are used to generate moment based features which are invariant to location size and optionally rotation an efficient method for computing the moments via geometric moments is presented a side effect of this method also yields scale invariance a new approach to location invariance using a minimum bounding circle is presented and a detailed analysis of the rotational properties of the moments is given data partitioning tests are performed to evaluate the various feature types and classifiers for rotational invariant character recognition the highest percentage of correctly classified characters was 91 7 and for non rotational invariant recognition it was 97 6 this compares with a previous effort using the same data and test conditions of 94 8 the techniques developed here should also be applicable to other areas of shape recognition

3d finite difference frequency domain modeling of visco acoustic wave propagation using a massively parallel direct solver a feasibility study we present a finite difference frequency domain method for 3d visco acoustic wave propagation modeling in the frequency domain the underlying numerical problem is the resolution of a large sparse system of linear equations whose right hand side term is the source this system is solved with a massively parallel direct solver we first present an optimal 3d finite difference stencil for frequency domain modeling the method is based on a parsimonious staggered grid method differential operators are discretized with second order accurate staggered grid stencils on different rotated coordinate systems to mitigate numerical anisotropy an antilumped mass strategy is implemented to minimize numerical dispersion the stencil incorporates 27 grid points and spans two grid intervals dispersion analysis shows that four grid points per wavelength provide accurate simulations in the 3d domain to assess the feasibility of the method for frequency domain full waveform inversion we computed simulations in the 3d seg eage overthrust model for frequencies 5 7 and 10 hz results confirm the huge memory requirement of the factorization several hundred figabytes but also the cpu efficiency of the resolution phase few seconds per shot heuristic scalability analysis suggests that the memory complexity of the factorization is o 35n 4 for a n 3 grid our method may provide a suitable tool to perform frequency domain full waveform inversion using a large distributed memory platform further investigation is still necessary to assess more quantitatively the respective merits and drawbacks of time and frequency domain modeling of wave propagation to perform 3d full waveform inversion

an optimal nearly analytic discrete weighted runge kutta discontinuous galerkin hybrid method for acoustic wavefield modeling abstractthe newly developed optimal nearly analytic discrete onad and the weighted runge kutta discontinuous galerkin wrkdg methods can effectively suppress the numerical dispersion caused by discretizing wave equations but it is difficult for onad to implement on flexible meshes whereas the wrkdg has high computational cost for wavefield simulations we have developed a new hybrid algorithm by combining the onad method with the wrkdg method in this hybrid algorithm the computational domain was split into several subdomains in which the subdomain for the onad method used regular cartesian grids whereas the subdomain for the wrkdg method used triangular grids the hybrid method was at least third order spatially accurate we have applied the proposed method to simulate the scalar wavefields for different models including a homogeneous model a rough topography model a fracture model and a cave model the numerical results found that the hybrid method can deal with complicated geometrical struc

on the solutions of the one dimensional ginzburg landau equations for superconductivity this paper gives a complete description of the solutions of the one dimensional ginzburg landau equations which model superconductivity phenomena in infinite slabs we investigate this problem over the entire range of physically important parameters a the size of the slab the ginzburg landau parameter and h0 the exterior magnetic field we do extensive numerical computations using the software auto and determine the number symmetry and stability of solutions for all values of the parameters in particular our experiments reveal the existence of two key points in parameter space which play a central role in the formation of the complicated patterns by means of bifurcation phenomena our global description also allows us to separate the various physically important regimes to classify previous results in each regime according to the values of the parameters and to derive new open problems in addition our investigation provides new insight into the problem of differentiating between the types of superconductors in terms of the parameters

optimally blended spectral finite element scheme for wave propagation and nonstandard reduced integration we study the dispersion and dissipation of the numerical scheme obtained by taking a weighted averaging of the consistent finite element mass matrix and lumped spectral element mass matrix for the small wave number limit we find and prove that for the optimum blending the resulting scheme a provides 2p 4 order accuracy for p th order method two orders more accurate compared with finite and spectral element schemes b has an absolute accuracy which is mathcal o p 3 and mathcal o p 2 times better than that of the pure finite and spectral element schemes respectively c tends to exhibit phase lag moreover we show that the optimally blended scheme can be efficiently implemented merely by replacing the usual gaussian quadrature rule used to assemble the mass and stiffness matrices by novel nonstandard quadrature rules which are also derived

a high order numerical approach with cartesian meshes for modeling of wave propagation and heat transfer on irregular domains with inhomogeneous materials abstract recently we have developed a new numerical approach for pdes with constant coefficients on irregular domains and cartesian meshes in this paper we extend it to a much more general case of pdes with variable coefficients that have a lot of applications e g the modeling of functionally graded materials the inhomogeneous materials obtained by 3 d printing and many others here we consider the 2 d wave and heat equations for isotropic and anisotropic inhomogeneous materials the idea of the extension to the case of pdes with variable coefficients is based on the representation of the stencil coefficients as functions of the mesh size this leads to the increase in the size of the local system of algebraic equations solved for each grid point of the new approach however this does not change the size of the global system of semidiscrete equations and practically does not increase the computational costs of the proposed technique similar to our previous technique the new 2 d approach with compact 9 point stencils uses trivial cartesian meshes for complex irregular domains and provides the fourth order of accuracy for the wave and heat equations with variable coefficients the calculation of the coefficients of the stencil equations is based on the minimization of the local truncation error of the stencil equations and yields the optimal order of accuracy of the new technique at similar 9 point stencils the accuracy of the new approach is much higher than that for the linear finite elements the numerical results for irregular domains show that at the same number of degrees of freedom the new approach is even much more accurate than the high order up to the third order finite elements with much wider stencils the wave and heat equations are uniformly treated with the new approach

a new family of finite difference schemes to solve the heterogeneous acoustic wave equation equivalent staggered grid scheme esg is a new family of schemes based on the finite difference method fdm the method is applied to acoustic wave propagation in variable density media and the results are compared with those from some classic fdm approaches the main feature of this new family is that it is designed to generate results numerically equivalent to those using the standard staggered grid formulations ssg but with the same memory requirements of simple grid schemes hence it results in a reduction of memory usage by 33 in 2d and 50 in 3d problems compared to the memory usage of ssg the first order system of equations in terms of pressure and velocity is not used here instead the formulation is based on applying new central difference operators to the secondorder acoustic wave equation in terms of pressure obtaining the same level of accuracy and stability as the ssg schemes the equivalence between the esg and ssg is mathematically demonstrated and issues concerning the application of seismic sources and the boundary conditions are addressed