integral_transform.cpp

Go to the documentation of this file.

#include "integral_transform.h"
#include <QFormLayout>
#include <QDoubleValidator>

IntegralTransform::IntegralTransform()
  : _widget(new QWidget()), ui(new Ui::IntegralTransform), _dT(0.0)
{
  ui->setupUi(_widget);
  ui->lineEditCustom->setValidator(
      new QDoubleValidator(0.0001, 1000, 4, ui->lineEditCustom));

  connect(ui->buttonCompute, &QPushButton::clicked, this,
          &IntegralTransform::on_buttonCompute_clicked);

  connect(ui->lineEditCustom, &QLineEdit::editingFinished, this, [=]() {
    _dT = ui->lineEditCustom->text().toDouble();
    emit parametersChanged();
  });

  connect(ui->radioActual, &QRadioButton::toggled, this, [=](bool toggled) {
    if (toggled)
    {
      _dT = 0.0;
      emit parametersChanged();
    }
  });

  connect(ui->radioCustom, &QRadioButton::toggled, this, [=](bool toggled) {
    if (toggled)
    {
      _dT = ui->lineEditCustom->text().toDouble();
      emit parametersChanged();
    }
  });
}

IntegralTransform::~IntegralTransform()
{
  delete ui;
  delete _widget;
}

std::optional<PlotData::Point> IntegralTransform::calculateNextPoint(size_t index)
{
  if (index == 0)
  {
    return {};
  }

  const auto& prev = dataSource()->at(index - 1);
  const auto& p = dataSource()->at(index);

  double dt = (_dT == 0.0) ? (p.x - prev.x) : _dT;

  if (dt <= 0)
  {
    return {};
  }

  double val = (p.y + prev.y) * dt / (2.0);
  _accumulated_value += val;
  PlotData::Point out = { p.x, _accumulated_value };
  return out;
}

QWidget* IntegralTransform::optionsWidget()
{
  const size_t data_size = dataSource()->size();

  if (!dataSource() || data_size < 2)
  {
    _widget->setEnabled(false);
  }
  return _widget;
}

void IntegralTransform::reset()
{
  _accumulated_value = 0.0;
  TransformFunction_SISO::reset();
}

bool IntegralTransform::xmlSaveState(QDomDocument& doc, QDomElement& parent_element) const
{
  QDomElement widget_el = doc.createElement("options");

  if (ui->radioActual->isChecked())
  {
    widget_el.setAttribute("radioChecked", "radioActual");
  }
  else
  {
    widget_el.setAttribute("radioChecked", "radioCustom");
  }
  widget_el.setAttribute("lineEdit", ui->lineEditCustom->text());

  parent_element.appendChild(widget_el);
  return true;
}

bool IntegralTransform::xmlLoadState(const QDomElement& parent_element)
{
  QDomElement widget_el = parent_element.firstChildElement("options");

  ui->lineEditCustom->setText(widget_el.attribute("lineEdit"));

  if (widget_el.attribute("radioChecked") == "radioActual")
  {
    ui->radioActual->setChecked(true);
  }
  else
  {
    ui->radioCustom->setChecked(true);
  }
  return true;
}

void IntegralTransform::on_buttonCompute_clicked()
{
  if (!dataSource() || dataSource()->size() < 2)
  {
    return;
  }

  const size_t data_size = dataSource()->size();

  // calculate automatic diff
  std::vector<double> diff;
  diff.reserve(data_size - 1);
  double prev_t = dataSource()->at(0).x;
  for (size_t i = 1; i < data_size; i++)
  {
    double t = dataSource()->at(i).x;
    double delta = t - prev_t;
    prev_t = t;
    diff.push_back(delta);
  }

  size_t first = 0;
  size_t last = diff.size();
  if (data_size > 10)
  {
    std::sort(diff.begin(), diff.end());
    first = last / 5;
    last = (last * 4) / 5;
  }
  double total = 0;
  for (size_t i = first; i < last; i++)
  {
    total += diff[i];
  }
  double estimated_dt = total / static_cast<double>(last - first);
  ui->lineEditCustom->setText(QString::number(estimated_dt, 'g', 4));

  if (ui->radioCustom->isChecked())
  {
    _dT = estimated_dt;
    emit parametersChanged();
  }
}