Scikit-ReduceModel is a Python package to build reduced order models. This code gives the standard reduced-basis framework and provides an efficient and accurate solution for model building. Also, an gives an extension of it, implementing the hp-greedy refinement strategy, an enhancement approach for reduced-basis model building. This approach uses an automatic parameter space partitioning, where there is a local reduced basis on each partition. The procedure splits spaces in two, in a recursive way.
Surrogate models can be built, which are comprised of three sequential steps of data processing that extract the most relevant information and use it to learn the patterns of the solutions and build easy to evaluate expressions: building a reduced basis, applying the empirical interpolation method and using a machine learning algorithm to learn the behavior of the solutions to obtain the surrogate model.
This package is built with the philosophy and idea of usability of scikit-learn modules. For usage examples, see the documentation.
In science and engineering is known that obtaining numerical simulations by solving differential equations can be more computational demanding than desired. For example, in the field of general relativity, to obtain expressions of gravitational waves could cost months using supercomputers. Furthermore, there are studies as parameter estimation that can require up to millions of sequential estimations, dominating the computational expense of the problem. In the last years, these problems were addresed building surrogate models from high cost simulations, taking advantage of the redundancy of the solutions with respect to the parameter space, which can build solutions in real time with no loss of accuracy.
To install the latest stable version of Scikit-ReducedModel from PyPI:
pip install scikit-reducedmodel
To install the developer version (may be unstable):
git clone https://github.com/francocerino/scikit-reducedmodel
cd scikit-reducedmodel
pip install .
In order to build a surrogate model, we need to be familiar with a set of functions parameterized by $λ$, denoted as $f_λ(x)$.
The known functions at given parameters are named as the training set (training_set
).
The associated parameters to training_set
are given by parameters
.
We need also a distretization of $x$ , given by x_set
.
Then, we can first build the reduced basis, in this case, we use the default parameters.
from skreducedmodel.reducedbasis import ReducedBasis
rb = ReducedBasis()
rb.fit(training_set = training_set,
parameters = parameters
physical_points = x_set)
In the second step, with the reduced basis built, the empirical interpolation method is applied.
from skreducedmodel.empiricalinterpolation import EmpiricalInterpolation
eim = EmpiricalInterpolation(rb)
eim.fit()
Finally, we build the reduced model from our eim object
from skreducedmodel.surrogate import Surrogate
surrogate = Surrogate(eim)
surrogate.fit()
In case we are interested in studying only Surrogate
objects, the package has a function that automates the whole process.
from skreducedmodel.mksurrogate import mksurrogate
surrogate = mksurrogate(parameters = param,
training_set = training_set,
physical_points = x_set,
)
To obtain an expression of a function $f´_{λ´}$ with any of the three methods, we must use:
rb.transform(f_test,parameter_test)
eim.transform(f_test,parameter_test)
surrogate.predict(parameter_test)
We encourage users to contribute with ideas, code, or by reporting bugs. To report bugs or issues, users should create an issue in the project repository. To contribute with code, please submit a pull request. We suggest that you contact us at francocerino@gmail.com prior to undertaking any significant improvement that requires substantial effort to address technical and design aspects before beginning development.
Franco Cerino [francocerino@gmail.com](francocerino@gmail.com) (FaMAF-UNC).
Agustín Rodríguez-Medrano (IATE-OAC-CONICET, FaMAF-UNC).