A tutorial on pseudospectral methods for computational optimal control

[English] This paper is a tutorial introduction to pseudospectral optimal control. With pseudospectral methods, a function is approximated as a linear combination of smooth basis functions, which are often chosen to be Legendre or Chebyshev polynomials. Collocation of the differential-algebraic equations is performed at orthogonal collocation points, which are selected to yield interpolation of high accuracy. Pseudospectral methods directly discretize the original optimal control problem to recast it into a nonlinear programming format. A numerical optimizer is then employed to find approximate local optimal solutions. The paper also briefly describes the functionality and implementation of PSOPT, an open source software package written in C++ that employs pseudospectral discretization methods to solve multi-phase optimal control problems. The software implements the Legendre and Chebyshev pseudospectral methods, and it has useful features such as automatic differentiation, sparsity detection, and automatic scaling. The use of pseudospectral methods is illustrated in two problems taken from the literature on computational optimal control. [Portuguese] Este artigo e um tutorial introdutorio sobre controle otimo pseudo-espectral. Em metodos pseudo-espectrais, uma funcao e aproximada como uma combinacao linear de funcoes de base suaves, tipicamente escolhidas como polinomios de Legendre ou Chebyshev. A colocacao de equacoes algebrico-diferenciais e realizada em pontos de colocacao ortogonal, que sao selecionados de modo a minimizar o erro de interpolacao. Metodos pseudoespectrais discretizam o problema de controle otimo original de modo a converte-lo em um problema de programa cao nao-linear. Um otimizador numerico e entao empregado para obter solucoes localmente otimas. Este artigo tambem descreve sucintamente a funcionalidade e a implementacao de um pacote computacional de codigo aberto escrito em C++ chamado PSOPT. Tal pacote emprega metodos de discretizacao pseudo-spectrais para resolver problemas de controle otimo com multiplas fase. O PSOPT permite a utilizacao de metodos de Legendre ou Chebyshev, e possui caractersticas uteis tais como diferenciacao automatica, deteccao de esparsidade e escalonamento automatico. O uso de metodos pseudo-espectrais e ilustrado em dois problemas retirados da literatura de controle otimo computacional.

Granular flow in the marginal ice zone

The region of sea ice near the edge of the sea ice pack is known as the marginal ice zone (MIZ), and its dynamics are complicated by ocean wave interaction with the ice cover, strong gradients in the atmosphere and ocean and variations in sea ice rheology. This paper focuses on the role of sea ice rheology in determining the dynamics of the MIZ. Here, sea ice is treated as a granular material with a composite rheology describing collisional ice floe interaction and plastic interaction. The collisional component of sea ice rheology depends upon the granular temperature, a measure of the kinetic energy of flow fluctuations. A simplified model of the MIZ is introduced consisting of the along and across momentum balance of the sea ice and the balance equation of fluctuation kinetic energy. The steady solution of these equations is found to leading order using elementary methods. This reveals a concentrated region of rapid ice flow parallel to the ice edge, which is in accordance with field observations, and previously called the ice jet. Previous explanations of the ice jet relied upon the existence of ocean currents beneath the ice cover. We show that an ice jet results as a natural consequence of the granular nature of sea ice.

Computational aspects of harmonic wavelet Galerkin methods and an application to a precipitation front propagation model

This article is dedicated to harmonic wavelet Galerkin methods for the solution of partial differential equations. Several variants of the method are proposed and analyzed, using the Burgers equation as a test model. The computational complexity can be reduced when the localization properties of the wavelets and restricted interactions between different scales are exploited. The resulting variants of the method have computational complexities ranging from O(N(3)) to O(N) (N being the space dimension) per time step. A pseudo-spectral wavelet scheme is also described and compared to the methods based on connection coefficients. The harmonic wavelet Galerkin scheme is applied to a nonlinear model for the propagation of precipitation fronts, with the front locations being exposed in the sizes of the localized wavelet coefficients. (C) 2011 Elsevier Ltd. All rights reserved.

A new re-ranking method using enhanced pseudo-relevance feedback for content-based medical image retrieval

We propose a novel re-ranking method for content-based medical image retrieval based on the idea of pseudo-relevance feedback (PRF). Since the highest ranked images in original retrieval results are not always relevant, a naive PRF based re-ranking approach is not capable of producing a satisfactory result. We employ a two-step approach to address this issue. In step 1, a Pearson's correlation coefficient based similarity update method is used to re-rank the high ranked images. In step 2, after estimating a relevance probability for each of the highest ranked images, a fuzzy SVM ensemble based approach is adopted to re-rank the images. The experiments demonstrate that the proposed method outperforms two other re-ranking methods.

Dynamic regime marginal structural mean models for estimation of optimal dynamic treatment regimes, Part I: main content

Dynamic treatment regimes are set rules for sequential decision making based on patient covariate history. Observational studies are well suited for the investigation of the effects of dynamic treatment regimes because of the variability in treatment decisions found in them. This variability exists because different physicians make different decisions in the face of similar patient histories. In this article we describe an approach to estimate the optimal dynamic treatment regime among a set of enforceable regimes. This set is comprised by regimes defined by simple rules based on a subset of past information. The regimes in the set are indexed by a Euclidean vector. The optimal regime is the one that maximizes the expected counterfactual utility over all regimes in the set. We discuss assumptions under which it is possible to identify the optimal regime from observational longitudinal data. Murphy et al. (2001) developed efficient augmented inverse probability weighted estimators of the expected utility of one fixed regime. Our methods are based on an extension of the marginal structural mean model of Robins (1998, 1999) which incorporate the estimation ideas of Murphy et al. (2001). Our models, which we call dynamic regime marginal structural mean models, are specially suitable for estimating the optimal treatment regime in a moderately small class of enforceable regimes of interest. We consider both parametric and semiparametric dynamic regime marginal structural models. We discuss locally efficient, double-robust estimation of the model parameters and of the index of the optimal treatment regime in the set. In a companion paper in this issue of the journal we provide proofs of the main results.

Estimação de regressões aditivas via backfitting e integração marginal: performance em pequenas amostras

Nesta dissertação realizou-se um experimento de Monte Carlo para re- velar algumas características das distribuições em amostras finitas dos estimadores Backfitting (B) e de Integração Marginal(MI) para uma regressão aditiva bivariada. Está-se particularmente interessado em fornecer alguma evidência de como os diferentes métodos de seleção da janela hn, tais co- mo os métodos plug-in, impactam as propriedades em pequenas amostras dos estimadores. Está-se interessado, também, em fornecer evidência do comportamento de diferentes estimadores de hn relativamente a seqüência ótima de hn que minimiza uma função perda escolhida. O impacto de ignorar a dependência entre os regressores na estimação da janela é tam- bém investigado. Esta é uma prática comum e deve ter impacto sobre o desempenho dos estimadores. Além disso, não há nenhuma rotina atual- mente disponível nos pacotes estatísticos/econométricos para a estimação de regressões aditivas via os métodos de Backfitting e Integração Marginal. É um dos objetivos a criação de rotinas em Gauss para a implementação prática destes estimadores. Por fim, diferentemente do que ocorre atual- mente, quando a utilização dos estimadores-B e MI é feita de maneira completamente ad-hoc, há o objetivo de fornecer a usuários informação que permita uma escolha mais objetiva de qual estimador usar quando se está trabalhando com uma amostra finita.

Bounds on policy relevant parameters with discrete policy variation

When estimating policy parameters, also known as treatment effects, the assignment to treatment mechanism almost always causes endogeneity and thus bias many of these policy parameters estimates. Additionally, heterogeneity in program impacts is more likely to be the norm than the exception for most social programs. In situations where these issues are present, the Marginal Treatment Effect (MTE) parameter estimation makes use of an instrument to avoid assignment bias and simultaneously to account for heterogeneous effects throughout individuals. Although this parameter is point identified in the literature, the assumptions required for identification may be strong. Given that, we use weaker assumptions in order to partially identify the MTE, i.e. to stablish a methodology for MTE bounds estimation, implementing it computationally and showing results from Monte Carlo simulations. The partial identification we perfom requires the MTE to be a monotone function over the propensity score, which is a reasonable assumption on several economics' examples, and the simulation results shows it is possible to get informative even in restricted cases where point identification is lost. Additionally, in situations where estimated bounds are not informative and the traditional point identification is lost, we suggest a more generic method to point estimate MTE using the Moore-Penrose Pseudo-Invese Matrix, achieving better results than traditional methods.

IMPROVEMENT TO THE MARGINAL COPING FIT of COMMERCIALLY PURE TITANIUM CAST IN PHOSPHATE-BONDED INVESTMENT BY USING A SIMPLE PATTERN COATING TECHNIQUE

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Accuracy of Digital Images in the Detection of Marginal Micro leakage: An In Vitro Study

Purpose: To evaluate the accuracy of Image Tool Software 3.0 (ITS 3.0) to detect marginal microleakage using the stereomicroscope as the validation criterion and ITS 3.0 as the tool under study.Materials and Methods: Class V cavities were prepared at the cementoenamel junction of 61 bovine incisors, and 53 halves of them were used. Using the stereomicroscope, microleakage was classified dichotomously: presence or absence. Next, ITS 3.0 was used to obtain measurements of the microleakage, so that 0.75 was taken as the cut-off point, and values equal to or greater than 0.75 indicated its presence, while values between 0.00 and 0.75 indicated its absence. Sensitivity and specificity were calculated by point and given as 95% confidence interval (95% CI).Results: The accuracy of the ITS 3.0 was verified with a sensitivity of 0.95 (95% CI: 0.89 to 1.00) and a specificity of 0.92 (95% CI: 0.84 to 0.99).Conclusion: Digital diagnosis of marginal microleakage using ITS 3.0 was sensitive and specific.

THE EFFECT of IMMEDIATE DENTIN SEALING on THE MARGINAL ADAPTATION and BOND STRENGTHS of TOTAL-ETCH and SELF-ETCH ADHESIVES

Statement of problem. Sealing ability and bond strengths of total-etch and self-etch dentin adhesives used for immediate dentin sealing have not been assessed and established.Purpose. The purpose of this study was to determine the effectiveness of immediate dentin sealing (IDS) using total-etch or self-etch dentin adhesives on microleakage and microtensile bond strength.Material and methods. Twenty recently extracted molars were selected, and standard MOD inlay preparations were made with the gingival margins located below the cemento-enamel unction. The teeth were assigned to 4 experimental groups (n=5) according to the indirect composite restoration cementation technique used: (1) immediate dentin sealing with Adper Single Bond (TEBI); (2) conventional adhesive cementation technique using Adper Single Bond (TEAI); (3) immediate dentin sealing using Adper Prompt L-Pop (SEBI); or (4) conventional adhesive cementation technique using Adper Prompt L-Pop (SEAI). The restored teeth were thermal cycled 1,000 times between 5 degrees and 55 degrees C and then immersed in 50% ammoniacal silver nitrate. Three specimens per restoration were evaluated for microleakage, according to predefined scores, and submitted to Friedman's test (alpha-.05). The specimens were then sectioned to obtain 0.8 +/- 0.2-mm-thick sticks (with n ranging from 32 to 57 specimens) and submitted to microtensile bond strength (mu TBS) testing. The obtained data were submitted to 2-way ANOVA test (alpha=.05).Results. None of the experimental groups demonstrated complete elimination of marginal microleakage. There were significant differences in microleakage of the tested adhesives (P>.001). IDS microleakage scores were similar to those obtained using the conventional cementation technique (CCT) for both adhesives. The highest mean bond strengths were obtained with TEBI (51.1 MPa), whereas SEAI showed the lowest mean bond strengths (1.7 MPa). IDS resulted in significantly higher bond strengths than CCT (P<.001).Conclusions. Total-etch and self-etch adhesives have a significant effect on IDS. IDS resulted in high bond strengths for both adhesives; however, the microleakage was similar to that obtained with CCT (J Prosthet Dent 2009;102:1-9)

Influence of the mineral content of the seston on tropical cladocerans of a marginal lake.

Bioassays were carried out seasonally to evaluate individual growth and reproduction of cladocerans, from a marginal lake, with the addition of nitrogen (N), phosphorus (P), and both N and P to natural seston; Methods: Cohorts originated from cultivated females were submitted to the following treatments: 1) lake seston, 2) lake seston + P, 3) lake seston + N, and 4) lake seston + NP; Results: The sestonic C:P and C:N molar ratios were always high and limiting, according to threshold ratios estimated for temperate lakes. P addition to seston enhanced the growth rates of one species, D. birgei. A significant higher growth rate of B. longirostris was found in the seston enriched with N compared to natural seston, as well as a higher fecundity of M. minuta. The fecundity of D. birgei was significantly higher in the seston enriched with both N and P. C, N, and P body content of cladocerans was similar to that of temperate counterparts; Conclusion: Energy limitation related to carbon content or food quality seems to be most important in controlling cladocerans' populations in the lake than food mineral content.

Occurrence of Pseudosida ramosa Daday, 1904 (Crustacea: Cladocera: Ctenopoda) in marginal lakes of the Paranapanema River, with comments on its distribution.

Aim: Contribute to the understanding of the distribution of a littoral cladoceran, Pseudosida ramosa Daday, 1904, in Brazil, with special reference for the Paranapanema River basin (Southeast Brazil); Methods: The specimens were sampled in two floodplain lakes located in the middle Paranapanema River basin, in the boundary of São Paulo and Parana States, using a standard plankton net (50 [mu]m mesh size). Trinocular microscopic with camera lucida and milimetric scale was used for the observation and measurements of the specimens; Results: Only three individuals were found (in winter and summer), despite of an extensive sampling program during 17 consecutive months of observations in both lakes. The largest body length was 1,146 [mu]m. The taxonomical identification followed Korovchinsky (1992); Conclusions: This species have a wide geographical distribution in Brazil, but it is restricted to zones colonized by aquatic macrophytes and seems to occur in low population abundance. The scarcity of registers of P. ramosa demonstrates the necessity of deeper investigations of littoral habits as well as the important role of this kind of environment for the conservation of the rivers/reservoirs biota, including some possible rare species.

Marginal and Internal Discrepancies Related to Margin Design of Ceramic Crowns Fabricated by a CAD/CAM System

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Vertical marginal discrepancy of ceramic copings with different ceramic materials finish lines and luting agents: An in vitro evaluation

Statement of problem. Prosthetic restorations that fit poorly may affect periodontal health and occlusion. Studies that have evaluated the accuracy of fit of ceramic restorations before and after cementation assessed primarily intracoronal restorations.Purpose. This in vitro study evaluated the effect of different finish lines, ceramic manufacturing techniques, and luting agents on the vertical discrepancy of ceramic copings.Material and methods. Two stainless steel molars were prepared for complete crowns with 2 different finish lines (heavy chamfer and rounded shoulder); each molar was duplicated to fabricate 90 copings. A total of 180 copings generated 18 groups (n=10 for each finish line-coping material-luting agent combination). Luting agents tested included zinc phosphate, resin-modified glass ionomer (Fuji Plus), and resin composite cements (Panavia F). A metal frame was developed on which to screw the stainless steel model and a ceramic coping; the distance (mum) between 2 predetermined points was measured before and after cementation by a profile projector under a torquing force. A 4-way ANOVA with repeated measurements was performed to assess the influence of each factor in the vertical marginal discrepancy: 3 between-coping factors (finish line-coping material-luting agent) and 1 within-coping factor (before and after cementation) (alpha=.05).Results. Procera copings presented the lowest mean values (P<.05) of vertical marginal discrepancy before and after cementation (25/44 mum) when compared to Empress 2 (68/110 mum) and InCeram Alumina copings (57/117 mum), regardless of any combinations among all finish lines and luting agents tested.Conclusion. Considering each factor separately, the ceramic manufacturing technique appeared to be the most important factor tested for the definitive vertical discrepancy of all-ceramic copings, with lower mean values for Procera copings.

Deproteinized bovine bone mineral in marginal defects at implants installed immediately into extraction sockets: an experimental study in dogs

Aim: To evaluate the influence of deproteinized bovine bone mineral (DBBM) particles concomitant with the placement of a collagen membrane on alveolar ridge preservation and on osseointegration of implants placed into alveolar sockets immediately after tooth extraction. Material and methods: The pulp tissue of the mesial roots of 3P3 was removed in six Labrador dogs and the root canals were filled. Flaps were elevated in the right side of the mandible, and the buccal and lingual alveolar bony plates were exposed. The third premolar was hemi-sectioned and the distal root was removed. A recipient site was prepared and an implant was placed lingually. After implant installation, defects of about 0.6mm wide and 3.1mm depth resulted at the buccal aspects of the implant, both at the test and at the control sites. The same surgical procedures and measurements were performed on the left side of the mandible. However, DBBM particles with a size of 0.25-1mm were placed into the remaining defect concomitant with the placement of a collagen membrane. Results: All implants were integrated into mature bone. No residual DBBM particles were detected at the test sites after 4 months of healing. Both the test and the control sites showed buccal alveolar bone resorption, 1.8 +/- 1.1 and 2.1 +/- 1mm, respectively. The most coronal bone-to-implant contact at the buccal aspect was 2 +/- 1.1 an 2.8 +/- 1.3mm, at the test and the control sites, respectively. This difference in the distance was statistically significant. Conclusion: The application of DBBM concomitant with a collagen membrane to fill the marginal defects around implants placed into the alveolus immediately after tooth extraction contributed to improved bone regeneration in the defects. However, with regard to buccal bony crest preservation, a limited contribution of DBBM particles was achieved.