An object-oriented optimization framework for large-scale inverse problems

2021 ◽  
pp. 104790
Author(s):  
Ettore Biondi ◽  
Guillaume Barnier ◽  
Robert G. Clapp ◽  
Francesco Picetti ◽  
Stuart Farris
2021 ◽  
Vol 47 (2) ◽  
pp. 1-34
Author(s):  
Umberto Villa ◽  
Noemi Petra ◽  
Omar Ghattas

We present an extensible software framework, hIPPYlib, for solution of large-scale deterministic and Bayesian inverse problems governed by partial differential equations (PDEs) with (possibly) infinite-dimensional parameter fields (which are high-dimensional after discretization). hIPPYlib overcomes the prohibitively expensive nature of Bayesian inversion for this class of problems by implementing state-of-the-art scalable algorithms for PDE-based inverse problems that exploit the structure of the underlying operators, notably the Hessian of the log-posterior. The key property of the algorithms implemented in hIPPYlib is that the solution of the inverse problem is computed at a cost, measured in linearized forward PDE solves, that is independent of the parameter dimension. The mean of the posterior is approximated by the MAP point, which is found by minimizing the negative log-posterior with an inexact matrix-free Newton-CG method. The posterior covariance is approximated by the inverse of the Hessian of the negative log posterior evaluated at the MAP point. The construction of the posterior covariance is made tractable by invoking a low-rank approximation of the Hessian of the log-likelihood. Scalable tools for sample generation are also discussed. hIPPYlib makes all of these advanced algorithms easily accessible to domain scientists and provides an environment that expedites the development of new algorithms.


1993 ◽  
Vol 2 (4) ◽  
pp. 133-144 ◽  
Author(s):  
Jon B. Weissman ◽  
Andrew S. Grimshaw ◽  
R.D. Ferraro

The conventional wisdom in the scientific computing community is that the best way to solve large-scale numerically intensive scientific problems on today's parallel MIMD computers is to use Fortran or C programmed in a data-parallel style using low-level message-passing primitives. This approach inevitably leads to nonportable codes and extensive development time, and restricts parallel programming to the domain of the expert programmer. We believe that these problems are not inherent to parallel computing but are the result of the programming tools used. We will show that comparable performance can be achieved with little effort if better tools that present higher level abstractions are used. The vehicle for our demonstration is a 2D electromagnetic finite element scattering code we have implemented in Mentat, an object-oriented parallel processing system. We briefly describe the application. Mentat, the implementation, and present performance results for both a Mentat and a hand-coded parallel Fortran version.


OOIS’97 ◽  
1998 ◽  
pp. 493-503 ◽  
Author(s):  
Arnd G. Grosse ◽  
Dietmar A. Kottmann ◽  
Jörn Hartroth
Keyword(s):  

Author(s):  
Zhao Sun ◽  
Yifu Wang ◽  
Lei Pan ◽  
Yunhong Xie ◽  
Bo Zhang ◽  
...  

AbstractPine wilt disease (PWD) is currently one of the main causes of large-scale forest destruction. To control the spread of PWD, it is essential to detect affected pine trees quickly. This study investigated the feasibility of using the object-oriented multi-scale segmentation algorithm to identify trees discolored by PWD. We used an unmanned aerial vehicle (UAV) platform equipped with an RGB digital camera to obtain high spatial resolution images, and multi-scale segmentation was applied to delineate the tree crown, coupling the use of object-oriented classification to classify trees discolored by PWD. Then, the optimal segmentation scale was implemented using the estimation of scale parameter (ESP2) plug-in. The feature space of the segmentation results was optimized, and appropriate features were selected for classification. The results showed that the optimal scale, shape, and compactness values of the tree crown segmentation algorithm were 56, 0.5, and 0.8, respectively. The producer’s accuracy (PA), user’s accuracy (UA), and F1 score were 0.722, 0.605, and 0.658, respectively. There were no significant classification errors in the final classification results, and the low accuracy was attributed to the low number of objects count caused by incorrect segmentation. The multi-scale segmentation and object-oriented classification method could accurately identify trees discolored by PWD with a straightforward and rapid processing. This study provides a technical method for monitoring the occurrence of PWD and identifying the discolored trees of disease using UAV-based high-resolution images.


2022 ◽  
Vol 14 (2) ◽  
pp. 273
Author(s):  
Mengyao Li ◽  
Rui Zhang ◽  
Hongxia Luo ◽  
Songwei Gu ◽  
Zili Qin

In recent years, the scale of rural land transfer has gradually expanded, and the phenomenon of non-grain-oriented cultivated land has emerged. Obtaining crop planting information is of the utmost importance to guaranteeing national food security; however, the acquisition of the spatial distribution of crops in large-scale areas often has the disadvantages of excessive calculation and low accuracy. Therefore, the IO-Growth method, which takes the growth stage every 10 days as the index and combines the spectral features of crops to refine the effective interval of conventional wavebands for object-oriented classification, was proposed. The results were as follows: (1) the IO-Growth method obtained classification results with an overall accuracy and F1 score of 0.92, and both values increased by 6.98% compared to the method applied without growth stages; (2) the IO-Growth method reduced 288 features to only 5 features, namely Sentinel-2: Red Edge1, normalized difference vegetation index, Red, short-wave infrared2, and Aerosols, on the 261st to 270th days, which greatly improved the utilization rate of the wavebands; (3) the rise of geographic data processing platforms makes it simple to complete computations with massive data in a short time. The results showed that the IO-Growth method is suitable for large-scale vegetation mapping.


2021 ◽  
Author(s):  
E. Biondi ◽  
G. Barnier ◽  
R. Clapp ◽  
F. Picetti ◽  
S. Farris

2012 ◽  
Author(s):  
Piotr J. Bandyk ◽  
Justin Freimuth ◽  
George Hazen

Object-oriented programming offers a natural approach to solving complex problems by focusing on individual aspects, or objects, and describing the ways in which they interact using interfaces. Modularity, extensibility, and code re-use often make OOP more appealing than its procedural counterpart. Code can be implemented in a more intuitive way and often mirrors the theory it derives from. Two examples are given in the form of real programs: a 3D panel code solver and a system-of-systems model for seabasing and environment sensing. Both are examples of large-scale frameworks and leverage the benefits offered by the object-oriented paradigm.


Author(s):  
A. Andrade-Campos

The use of optimization methods in engineering is increasing. Process and product optimization, inverse problems, shape optimization, and topology optimization are frequent problems both in industry and science communities. In this paper, an optimization framework for engineering inverse problems such as the parameter identification and the shape optimization problems is presented. It inherits the large experience gain in such problems by the SiDoLo code and adds the latest developments in direct search optimization algorithms. User subroutines in Sdl allow the program to be customized for particular applications. Several applications in parameter identification and shape optimization topics using Sdl Lab are presented. The use of commercial and non-commercial (in-house) Finite Element Method codes to evaluate the objective function can be achieved using the interfaces pre-developed in Sdl Lab. The shape optimization problem of the determination of the initial geometry of a blank on a deep drawing square cup problem is analysed and discussed. The main goal of this problem is to determine the optimum shape of the initial blank in order to save latter trimming operations and costs.


Sign in / Sign up

Export Citation Format

Share Document