gemseo_calibration.post.data_versus_model.post

# Copyright 2021 IRT Saint Exupéry, https://www.irt-saintexupery.com
#
# This program is free software; you can redistribute it and/or
# modify it under the terms of the GNU Lesser General Public
# License version 3 as published by the Free Software Foundation.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
# Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with this program; if not, write to the Free Software Foundation,
# Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
"""Plot the model data versus the reference data."""
from __future__ import annotations

import logging

from gemseo.core.dataset import Dataset

from gemseo_calibration.post.multiple_scatter import MultipleScatter
from gemseo_calibration.post_processor import CalibrationPostProcessor

LOGGER = logging.getLogger(__name__)


[docs]class DataVersusModel(CalibrationPostProcessor):
    """Scatter plot of the model data versus the reference ones."""

    def _plot(
        self,
        output: str,
    ) -> None:
        opt_name = f"Opt[{output}]"
        init_name = f"Init[{output}]"
        ref_name = f"Ref[{output}]"

        dataset = Dataset()
        dataset.add_variable(
            opt_name,
            self._posterior_model_data.get_data_by_names(output, False).mean(1)[
                :, None
            ],
        )
        dataset.add_variable(
            init_name,
            self._prior_model_data.get_data_by_names(output, False).mean(1)[:, None],
        )
        dataset.add_variable(
            ref_name,
            self._reference_data.get_data_by_names(output, False).mean(1)[:, None],
        )
        plot = MultipleScatter(dataset, x=ref_name, y=[init_name, opt_name])
        plot.color = ["blue", "red"]
        plot.xlabel = "Reference"
        plot.ylabel = "Model"
        plot.labels = {opt_name: "After calibration", init_name: "Before calibration"}
        figures = plot.execute(save=False, show=False)
        for figure in figures:
            self._add_figure(figure)
Source code for gemseo_calibration.post.data_versus_model.post
