the adiabatic mode parabolic equation an efficient method is developed for solving problems involving three dimensional acoustic propagation in a thin spherical shell the method is based on the adiabatic mode approximation which assumes that energy coupling between modes is negligible and the parabolic equation pe method which efficiently handles spatial variations in the acoustic properties the pe method is used to solve two dimensional wave equations for the adiabatic mode coefficients over latitude and longitude due to this reduction in dimension the adiabatic mode pe should be practical for solving global scale ocean acoustics problems these types of problems which are currently of great interest heaney et al j acoust soc am 90 2586 2594 1991 baggeroer and munk phys today 45 9 22 30 1992 are too large to solve with other existing three dimensional models

an energy conserving parabolic equation for elastic media the accuracy of the parabolic equation pe method is improved for range dependent elastic media the energy conserving elastic pe is a generalization of the energy conserving acoustic pe a range dependent waveguide is approximated by a sequence of range independent regions a linear approximation is derived for the nonlinear condition corresponding to conservation of compressional energy flux at the vertical interfaces between range independent regions the rotated elastic pe is applied to generate reference solutions for demonstrating accuracy the energy conserving elastic pe is applied to solve problems involving the propagation of an interface wave along a sloping ocean bottom and a problem involving mode cutoff and coupling into shear wave beams in the ocean bottom the latter problem is analogous to a cutoff problem that was first solved with the acoustic pe f b jensen and w a kuperman j acoust soc am 67 1564 1566 1980

the lagrangian relaxation method for solving integer programming problems this article originally appeared in management science january 1981 volume 27 number 1 pp 1 18 published by the institute of management sciences one of the most computationally useful ideas of the 1970s is the observation that many hard integer programming problems can be viewed as easy problems complicated by a relatively small set of side constraints dualizing the side constraints produces a lagrangian problem that is easy to solve and whose optimal value is a lower bound for minimization problems on the optimal value of the original problem the lagrangian problem can thus be used in place of a linear programming relaxation to provide bounds in a branch and bound algorithm this approach has led to dramatically improved algorithms for a number of important problems in the areas of routing location scheduling assignment and set covering this paper is a review of lagrangian relaxation based on what has been learned in the last decade

projected newton methods and optimization of multicommodity flows a superlinearly convergent newton like method for linearly constrained optimization problems is adapted for solution of multicommodity network flow problems of the type arising in communication and transportation networks we show that the method can be implemented approximately by making use of conjugate gradient iterations without the need to compute explicitly the hessian matrix preliminary computational results suggest that this type of method is capable of yielding highly accurate solutions of nonlinear multicommodity flow problems far more efficiently than any of the methods available at present

a mapping approach for handling sloping interfaces a mapping approach for handling sloping interfaces in parabolic equation solutions is developed and tested at each range the medium is rigidly translated vertically so that a sloping interface becomes horizontal to simplify the approach the slope is assumed to be small and the extra terms that arise in the wave equation under the mapping are neglected the effects of these terms can be approximately accounted for by applying a leading order correction to the phase the mapping introduces variations in topography which are relatively easy to handle for the case of a pressure release boundary condition the accuracy of the approach is demonstrated for problems involving fluid sediments the approach should also be accurate for problems involving elastic sediments and should be useful for solving three dimensional problems involving variable topography

making and evaluating point forecasts typically point forecasting methods are compared and assessed by means of an error measure or scoring function with the absolute error and the squared error being key examples the individual scores are averaged over forecast cases to result in a summary measure of the predictive performance such as the mean absolute error or the mean squared error i demonstrate that this common practice can lead to grossly misguided inferences unless the scoring function and the forecasting task are carefully matched effective point forecasting requires that the scoring function be specified ex ante or that the forecaster receives a directive in the form of a statistical functional such as the mean or a quantile of the predictive distribution if the scoring function is specified ex ante the forecaster can issue the optimal point forecast namely the bayes rule if the forecaster receives a directive in the form of a functional it is critical that the scoring function be consistent for it in the sense that t

finite difference solution to the parabolic wave equation in this article an implicit finite difference ifd scheme is introduced for solving underwater acoustic wave propagation problems this scheme is applied to the parabolic wave equation model for the acoustic problem with a free surface and arbitrary bottom and bottom boundary conditions the mathematical theory of the implicit finite difference method as applied to the parabolic equation is developed and a computer model to implement the ifd is introduced representative problems are presented to demonstrate the capability of the ifd model advantages and limitations of the ifd are discussed

a two way parabolic equation for acoustic backscattering in the ocean the parabolic equation pe method is generalized to handle backscattered acoustic energy in the ocean the two way pe is based on the single scattering approximation and the approach of two way coupled modes in which range dependent environments are approximated by a sequence of range independent regions at the vertical boundaries between regions the solution of the two way pe is required to satisfy two continuity conditions the range derivative in one of the conditions is replaced by a higher order pe depth operator the reflected and transmitted fields that satisfy these conditions are computed with an efficient iteration scheme the outgoing and incoming fields are propagated by two way range marching the two way pe which is presently implemented for two dimensional problems is a practical method for solving large scale reverberation problems the accuracy of the two way pe is demonstrated by comparisons with reference solutions the two way pe is applied to simulate the localization of a source of backscattering using the method of back propagation

scalar electromagnetic propagation modelling using parabolic equations and the split step pade approximation this paper presents a solution of 2d electromagnetic wave propagation problems in complicated terrestrial domains scalar 2d parabolic approximations are derived from maxwell s equations for both vertical and horizontal polarization the parabolic equations are then solved using a new technique involving pade rational function approximations of the macroscopic operator this method allows for larger than normal range stepping speeding up computational time significantly the pade approximation and the numerical implementation are fully discussed the discontinuity at the earth s surface is handled directly by using classical continuity conditions and deriving exact interface conditions for linking the fields in the atmosphere to those in the terrain the interface conditions are then implemented using the concept of virtual points preliminary benchmark tests show the interface treatment to work well finally several example runs are presented illustrating results

a higher order parabolic equation for wave propagation in an ocean overlying an elastic bottom a higher order elastic parabolic equation hepe is derived for wave propagation in depth dependent and weakly range dependent fluid solid media galerkin s method is used to discretize the depth operators in the hepe within layers in which depth variations in the lame constants and density are continuous discontinuities in material properties are handled with centered differences for the interface conditions between layers the numerical solution of the hepe also involves the method of alternating directions and crank nicolson integration the hepe is applied to underwater wave propagation problems involving a water column over an elastic bottom including a weakly range dependent problem the accuracy of the hepe is demonstrated with benchmark calculations

of quantiles and expectiles consistent scoring functions choquet representations and forecast rankings summary in the practice of point prediction it is desirable that forecasters receive a directive in the form of a statistical functional for example forecasters might be asked to report the mean or a quantile of their predictive distributions when evaluating and comparing competing forecasts it is then critical that the scoring function used for these purposes be consistent for the functional at hand in the sense that the expected score is minimized when following the directive we show that any scoring function that is consistent for a quantile or an expectile functional can be represented as a mixture of elementary or extremal scoring functions that form a linearly parameterized family scoring functions for the mean value and probability forecasts of binary events constitute important examples the extremal scoring functions admit appealing economic interpretations of quantiles and expectiles in the context of betting and investment problems the choquet type mixture representations give rise to simple checks of whether a forecast dominates another in the sense that it is preferable under any consistent scoring function in empirical settings it suffices to compare the average scores for only a finite number of extremal elements plots of the average scores with respect to the extremal scoring functions which we call murphy diagrams permit detailed comparisons of the relative merits of competing forecasts

on the effect of error correlation on matched field geoacoustic inversion the effect of correlated data errors on matched field geoacoustic inversion for vertical array data is examined the correlated errors stem from our inability to model the inhomogeneities in the environment resulting in an additional error term beyond ambient noise simulated data with these correlated errors are generated and then inverted with or without using the proper covariance matrix results show that the correlated error has a negative impact on geoacoustic parameter estimation if not accounted for properly

bias corrected score decomposition for generalized quantiles summarydecompositions of the score of a forecast represent useful tools for assessing its performance we consider local score decompositions permitting detailed forecast assessments across a spectrum of conditions of interest we derive corrections to the bias of the decomposition components in the framework of point forecasts of quantile type functionals and illustrate their performance by simulation related bias corrections have thus far only been known for squared error criteria

a single scattering correction for the seismo acoustic parabolic equation an efficient single scattering correction that does not require iterations is derived and tested for the seismo acoustic parabolic equation the approach is applicable to problems involving gradual range dependence in a waveguide with fluid and solid layers including the key case of a sloping fluid solid interface the single scattering correction is asymptotically equivalent to a special case of a single scattering correction for problems that only have solid layers kusel et al j acoust soc am 121 808 813 2007 the single scattering correction has a simple interpretation conservation of interface conditions in an average sense that facilitated its generalization to problems involving fluid layers promising results are obtained for problems in which the ocean bottom interface has a small slope

rational operators for filtering the parabolic equation pe method has been applied to approximate the wave equation construct radiation conditions solve scattering problems construct initial conditions and derive energy flux conditions a new application of the pe method filters based on rational approximations of depth operators is described and tested the rational approximation is designed to act as the identity operator on the desired part of the spectrum of the depth separated wave equation and to annihilate other parts of the spectrum this approach does not require explicit knowledge of the spectrum the applications of the filter include directly solving eigenvalue problems annihilating components of the wave number spectrum generating initial conditions at the source range and eliminating gibbs oscillations that arise in energy conserving pe solutions

a logical topology and discrete capacity assignment algorithm for reconfigurable networks a reconfigurable network is a circuit switched network where the effective topology and capacities can be dynamically adapted to changes in the traffic requirements or to changes in the structure of the network due to failures the application of reconfigurable networks considered in this paper is data communications we formulated the joint topology capacity and routing problem in a reconfigurable network as a nonlinear mixed integer programming problem to solve this problem we present a partial branch and bound algorithm the reduced gradient method is used to solve the problem with the integrality constraints relaxed lagrangian relaxation is applied to obtain lower bounds in the computational experiments the algorithm found good solutions in a few minutes of cpu time in addition the reduction in the minimal delay due to the reconfiguration capability can be as large as 62

complex pade approximants for wide angle acoustic propagators procedures for constructing wide angle propagators that attenuate evanescent modes in parabolic equation pe models of underwater sound propagation are examined in this context two new split step pade approximations to exponentiated square root operators are proposed the first method involves a propagator derived from a pade approximant of a square root operator that has been rotated in the complex plane while the propagator of the second method is obtained from a complex coefficient rational approximation of the square root operator numerical evidence confirms that these proposed methods yield improved behavior compared to previously suggested techniques in selected areas of the spectral domain

a split step pade solution for the parabolic equation method a split step pade solution is derived for the parabolic equation pe method higher order pade approximations are used to reduce both numerical errors and asymptotic errors e g phase errors due to wide angle propagation this approach is approximately two orders of magnitude faster than solutions based on pade approximations that account for asymptotic errors but not numerical errors in contrast to the split step fourier solution which achieves similar efficiency for some problems the split step pade solution is valid for problems involving very wide propagation angles large depth variations in the properties of the waveguide and elastic ocean bottoms the split step pade solution is practical for global scale problems

asymmetric least squares estimation and testing this paper considers estimation and hypothesis tests for coefficients of linear regression models where the coefficient estimates are based on location measures defined by an asymmetric least squares criterion function these asymmetric least squares estimators have properties which are analogous to regression quantile estimators but are much simpler to calculate as are the corresponding test statistics the coefficient estimators can be used to construct test statistics for homoskedasticity and conditional symmetry of the error distribution and we find these tests compare quite favorably with other commonly used tests of these null hypotheses in terms of local relative efficiency consequently asymmetric least squares estimation provides a convenient and relatively efficient method of summarizing the conditional distribution of a dependent variable given the regressors and a means of testing whether a linear model is an adequate characterization of the typical value for this conditional distribution

energy conserving and reciprocal solutions for higher order parabolic equations conservation of sound energy in nonabsorbing media and acoustical reciprocity in stationary fluids at rest are the fundamental properties of the acoustic field preserving these properties in parabolic approximations is essential for accurate modeling of sound propagation in shallow water environments recently godin proposed a family of wide angle parabolic equations pe s for motionless fluid which are energy conserving reciprocal and allow improved description of mode coupling without adversely affecting phase accuracy for moving media a generalized claerbout pe was also derived its solutions preserve the acoustic energy in the absence of dissipation and satisfy the flow reversal theorem frt which is an extension of reciprocity principle to moving media the report presents a numerical scheme for solving the proposed pe s the algorithm preserves the conservational properties of the original pseudo differential equations implicit finite difference solutions satisfy discrete analogs of the e

probabilistic forecasts calibration and sharpness summary probabilistic forecasts of continuous variables take the form of predictive densities or predictive cumulative distribution functions we propose a diagnostic approach to the evaluation of predictive performance that is based on the paradigm of maximizing the sharpness of the predictive distributions subject to calibration calibration refers to the statistical consistency between the distributional forecasts and the observations and is a joint property of the predictions and the events that materialize sharpness refers to the concentration of the predictive distributions and is a property of the forecasts only a simple theoretical framework allows us to distinguish between probabilistic calibration exceedance calibration and marginal calibration we propose and study tools for checking calibration and sharpness among them the probability integral transform histogram marginal calibration plots the sharpness diagram and proper scoring rules the diagnostic approach is illustrated by an assessment and ranking of probabilistic forecasts of wind speed at the stateline wind energy centre in the us pacific northwest in combination with cross validation or in the time series context our proposal provides very general nonparametric alternatives to the use of information criteria for model diagnostics and model selection

a higher order energy conserving parabolic equqation for range dependent ocean depth sound speed and density outgoing solutions of the wave equation including parabolic equation pe and normal mode solutions are usually formulated so that pressure is continuous with range for range dependent problems the accuracy of normal mode solutions has been improved by conserving energy rather than maintaining continuity of pressure porter et al the problem of energy conservation in one way equations j acoust soc am 89 1058 1067 1991 this approach is applied to derive a higher order energy conserving pe that provides improved accuracy for problems involving large ocean bottom slopes and large range and depth variations in sound speed and density a special numerical approach and complex pade coefficients are applied to suppress gibbs oscillations the back propagated half space field an improved pe starter is applied to handle wide propagation angles reference solutions generated with a complex ray model and with the rotated pe are used for comparison