Integration::Integration(const QString& formula, const QString& var, ApplicationWindow *parent, Graph *g, double start, double end)
: Filter(parent, g),
d_formula(formula),
d_variable(var)
{
d_init_err = false;
d_n = 0;
d_from = start;
d_to = end;
if (d_to == d_from)
d_init_err = true;
MyParser parser;
double x = 0.0;
parser.DefineVar(d_variable.ascii(), &x);
parser.SetExpr(d_formula.ascii());
try {
parser.Eval();
} catch(mu::ParserError &e) {
QMessageBox::critical(parent, tr("QtiPlot - Input error"), QString::fromStdString(e.GetMsg()));
d_init_err = true;
}
setObjectName(tr("Integration"));
d_integrand = AnalyticalFunction;
d_method = 1;
d_max_iterations = 20;
d_sort_data = false;
}
double Integration::trapezf(int n)
{
MyParser parser;
double x = d_from;
parser.DefineVar(d_variable.ascii(), &x);
parser.SetExpr(d_formula.ascii());
static double s;
if (n == 1){
x = d_from;
double aux = parser.Eval();
x = d_to;
return (s = 0.5*(d_to - d_from)*(aux + parser.Eval()));
} else {
int it = 1;
for(int j=1; j < n-1; j++)
it<<=1;
double tnm = it;
double del = (d_to - d_from)/tnm;
x = d_from + 0.5*del;
double sum = 0.0;
for(int j=1; j <= it; j++, x += del)
sum += parser.Eval();
s = 0.5*(s + (d_to - d_from)*sum/tnm);
return s;
}
}
double ScaleDraw::transformValue(double value) const
{
if (!formula_string.isEmpty())
{
double lbl=0.0;
try
{
MyParser parser;
if (formula_string.contains("x"))
parser.DefineVar("x", &value);
else if (formula_string.contains("y"))
parser.DefineVar("y", &value);
parser.SetExpr(formula_string.ascii());
lbl=parser.Eval();
}
catch(mu::ParserError &)
{
return 0;
}
return lbl;
}
else
return value;
}
/** Checks to see if this axis has valid parameters
*
*/
bool AxisDetails::valid()
{
if (m_cmbAxisType->currentIndex() == ScaleDraw::Numeric)
{
if (m_chkShowFormula->isChecked())
{
QString formula = m_txtFormula->text().lower();
try
{
double value = 1.0;
MyParser parser;
if (formula.contains("x"))
{
parser.DefineVar("x", &value);
}
else if (formula.contains("y"))
{
parser.DefineVar("y", &value);
}
parser.SetExpr(formula.ascii());
parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(this, tr("MantidPlot - Formula input error"), QString::fromStdString(e.GetMsg())+"\n"+tr("Valid variables are 'x' for Top/Bottom axes and 'y' for Left/Right axes!"));
return false;
}
}
}
Table *w = m_app->table(m_cmbColName->currentText());
return m_initialised && m_graph && !((m_cmbAxisType->currentIndex() == ScaleDraw::Text || m_cmbAxisType->currentIndex() == ScaleDraw::ColHeader) && !w);
}
void DoubleSpinBox::interpretText()
{
bool ok = false;
QString s = text();
double value = locale().toDouble(s, &ok);
if (ok && value == d_value)
return;
if (!ok){
MyParser parser;
parser.setLocale(QLocale());
parser.addGSLConstants();
try {
parser.SetExpr(s.toAscii().constData());
value = parser.Eval();
} catch (mu::ParserError &e){
lineEdit()->setText(textFromValue(d_value));
return;
}
}
if (setValue(value))
emit valueChanged(d_value);
else
lineEdit()->setText(textFromValue(d_value));
}
void NonLinearFit::removeDataSingularities()
{
MyParser parser;
for (int i = 0; i < d_p; i++){
double param = gsl_vector_get(d_param_init, i);
parser.DefineVar(d_param_names[i].ascii(), ¶m);
}
QMapIterator<QString, double> it(d_constants);
while (it.hasNext()) {
it.next();
parser.DefineConst(it.key().ascii(), it.value());
}
double xvar;
parser.DefineVar("x", &xvar);
parser.SetExpr(d_formula.ascii());
for (int i = 0; i < d_n; i++){
xvar = d_x[i];
try {
parser.EvalRemoveSingularity(&xvar);
} catch(MyParser::Pole){
QApplication::restoreOverrideCursor();
QMessageBox::critical((ApplicationWindow *)parent(), QObject::tr("QtiPlot"),
QObject::tr("Ignored data point at x = %1.").arg(xvar));
removePole(i);
}
}
}
void NonLinearFit::calculateFitCurveData(double *X, double *Y)
{
MyParser parser;
for (int i=0; i<d_p; i++)
parser.DefineVar(d_param_names[i].ascii(), &d_results[i]);
QMapIterator<QString, double> i(d_constants);
while (i.hasNext()) {
i.next();
parser.DefineConst(i.key().ascii(), i.value());
}
double x;
parser.DefineVar("x", &x);
parser.SetExpr(d_formula.ascii());
if (d_gen_function){
double X0 = d_x[0];
double step = (d_x[d_n - 1] - X0)/(d_points - 1);
for (int i=0; i<d_points; i++){
x = X0 + i*step;
X[i] = x;
Y[i] = parser.EvalRemoveSingularity(&x, false);
}
} else {
for (int i=0; i<d_points; i++) {
x = d_x[i];
X[i] = x;
Y[i] = parser.EvalRemoveSingularity(&x, false);
}
}
}
bool NonLinearFit::setFormula(const QString& s, bool guess)
{
if (s.isEmpty()){
QMessageBox::critical((ApplicationWindow *)parent(), tr("QtiPlot - Input function error"),
tr("Please enter a valid non-empty expression! Operation aborted!"));
d_init_err = true;
return false;
}
if (d_formula == s)
return true;
if (guess)
setParametersList(guessParameters(s));
if (!d_p){
QMessageBox::critical((ApplicationWindow *)parent(), tr("QtiPlot - Fit Error"),
tr("There are no parameters specified for this fit operation. Please define a list of parameters first!"));
d_init_err = true;
return false;
}
try {
double *param = new double[d_p];
MyParser parser;
double xvar;
parser.DefineVar("x", &xvar);
for (int k = 0; k < (int)d_p; k++){
param[k] = gsl_vector_get(d_param_init, k);
parser.DefineVar(d_param_names[k].ascii(), ¶m[k]);
}
QMapIterator<QString, double> i(d_constants);
while (i.hasNext()) {
i.next();
parser.DefineConst(i.key().ascii(), i.value());
}
parser.SetExpr(s.ascii());
parser.Eval() ;
delete[] param;
} catch(mu::ParserError &e){
QMessageBox::critical((ApplicationWindow *)parent(), tr("QtiPlot - Input function error"), QString::fromStdString(e.GetMsg()));
d_init_err = true;
return false;
}
d_init_err = false;
d_formula = s;
return true;
}
double NonLinearFit::eval(double *par, double x)
{
MyParser parser;
for (int i=0; i<d_p; i++)
parser.DefineVar(d_param_names[i].ascii(), &par[i]);
QMapIterator<QString, double> i(d_constants);
while (i.hasNext()) {
i.next();
parser.DefineConst(i.key().ascii(), i.value());
}
parser.DefineVar("x", &x);
parser.SetExpr(d_formula.ascii());
return parser.EvalRemoveSingularity(&x, false);
}
double user_d(const gsl_vector * x, void *params) {
int n = ((struct FitData *)params)->n;
int p = ((struct FitData *)params)->p;
double *X = ((struct FitData *)params)->X;
double *Y = ((struct FitData *)params)->Y;
double *sigma = ((struct FitData *)params)->sigma;
NonLinearFit *fitter = (NonLinearFit *)((struct FitData *) params)->fitter;
const char *function = fitter->formula().ascii();
QStringList parNames = fitter->parameterNames();
double val=0;
MyParser parser;
try {
double *parameters = new double[p];
double xvar;
parser.DefineVar("x", &xvar);
for (int i=0; i < p; i++) {
parameters[i]=gsl_vector_get(x,i);
parser.DefineVar(parNames[i].ascii(), ¶meters[i]);
}
QMapIterator<QString, double> i(fitter->constants());
while (i.hasNext()){
i.next();
parser.DefineConst(i.key().ascii(), i.value());
}
parser.SetExpr(function);
for (int j = 0; j < n; j++) {
xvar = X[j];
double s = 1.0/sqrt(sigma[j]);
try {
double t = (parser.EvalRemoveSingularity(&xvar) - Y[j])/s;
val += t*t;
} catch (MyParser::Pole) {
return GSL_POSINF; //weird, I know. blame gsl.
}
}
delete[] parameters;
} catch (mu::ParserError &e) {
QMessageBox::critical(0,"QtiPlot - Input function error",QString::fromStdString(e.GetMsg()));
return GSL_EINVAL;
}
return val;
}
double UserFunction2D::operator()(double x, double y) {
if (d_formula.isEmpty())
return 0.0;
MyParser parser;
double result = 0.0;
try {
parser.DefineVar("x", &x);
parser.DefineVar("y", &y);
parser.SetExpr((const std::string)d_formula.toAscii().constData());
result = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(nullptr, "MantidPlot - Input function error",
QString::fromStdString(e.GetMsg()));
}
return result;
}
int user_df(const gsl_vector *x, void *params, gsl_matrix *J) {
int n = ((struct FitData *)params)->n;
int p = ((struct FitData *)params)->p;
double *X = ((struct FitData *)params)->X;
double *sigma = ((struct FitData *)params)->sigma;
NonLinearFit *fitter = (NonLinearFit *)((struct FitData *) params)->fitter;
const char *function = fitter->formula().ascii();
QStringList parNames = fitter->parameterNames();
try {
double *param = new double[p];
MyParser parser;
double xvar;
parser.DefineVar("x", &xvar);
for (int k=0; k<p; k++) {
param[k] = gsl_vector_get(x,k);
parser.DefineVar(parNames[k].ascii(), ¶m[k]);
}
QMapIterator<QString, double> i(fitter->constants());
while (i.hasNext()){
i.next();
parser.DefineConst(i.key().ascii(), i.value());
}
parser.SetExpr(function);
for (int i = 0; i < n; i++) {
xvar = X[i];
double s = 1.0/sqrt(sigma[i]);
for (int j = 0; j < p; j++)
try {
gsl_matrix_set (J, i, j, 1.0/s*parser.DiffRemoveSingularity(&xvar, ¶m[j], param[j]));
} catch (MyParser::Pole) {
return GSL_ESING;
}
}
delete[] param;
} catch (mu::ParserError &) {
return GSL_EINVAL;
}
return GSL_SUCCESS;
}
double ScaleDraw::transformValue(double value) const
{
if (!d_formula.isEmpty()){
double lbl=0.0;
try{
MyParser parser;
if (d_formula.contains("x", Qt::CaseInsensitive))
parser.DefineVar("x", &value);
else if (d_formula.contains("y", Qt::CaseInsensitive))
parser.DefineVar("y", &value);
parser.SetExpr(d_formula.lower().ascii());
lbl = parser.Eval();
}
catch(mu::ParserError &){
return 0;
}
return lbl;
} else
return value;
}
bool Plot3DDialog::updatePlot()
{
int axis=-1;
if (generalDialog->currentWidget()==(QWidget*)bars)
{
emit updateBars(boxBarsRad->text().toDouble());
}
if (generalDialog->currentWidget()==(QWidget*)points)
{
if (boxPointStyle->currentItem() == 0)
emit updatePoints(boxSize->text().toDouble(), boxSmooth->isChecked());
else if (boxPointStyle->currentItem() == 1)
emit updateCross(boxCrossRad->text().toDouble(), boxCrossLinewidth->text().toDouble(),
boxCrossSmooth->isChecked(), boxBoxed->isChecked());
else if (boxPointStyle->currentItem() == 2)
emit updateCones(boxConesRad->text().toDouble(), boxQuality->value());
}
if (generalDialog->currentWidget()==(QWidget*)title)
{
emit updateTitle(boxTitle->text(),titleColor,titleFont);
}
if (generalDialog->currentWidget()==(QWidget*)colors)
{
emit updateTransparency(boxTransparency->value()*0.01);
emit updateDataColors(fromColor,toColor);
emit updateColors(meshColor,axesColor,numColor,labelColor,bgColor,gridColor);
}
if (generalDialog->currentWidget()==(QWidget*)general)
{
emit showColorLegend(boxLegend->isChecked());
emit updateMeshLineWidth(boxMeshLineWidth->value());
emit adjustLabels(boxDistance->value());
emit updateResolution (boxResolution->value());
emit showColorLegend(boxLegend->isChecked());
emit setNumbersFont(numbersFont);
emit updateZoom(boxZoom->value()*0.01);
emit updateScaling(boxXScale->value()*0.01,boxYScale->value()*0.01,
boxZScale->value()*0.01);
}
if (generalDialog->currentWidget()==(QWidget*)scale)
{
axis=axesList->currentRow();
QString from=boxFrom->text().toLower();
QString to=boxTo->text().toLower();
double start,end;
bool error=false;
try
{
MyParser parser;
parser.SetExpr(from.toAscii().constData());
start=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(0,tr("Start limit error"), QString::fromStdString(e.GetMsg()));
boxFrom->setFocus();
error=true;
return false;
}
try
{
MyParser parser;
parser.SetExpr(to.toAscii().constData());
end=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(0,tr("End limit error"), QString::fromStdString(e.GetMsg()));
boxTo->setFocus();
error=true;
return false;
}
if (start>=end)
{
QMessageBox::critical(0,tr("Input error"),
tr("Please enter scale limits that satisfy: from < to!"));
boxTo->setFocus();
return false;
}
if (! error)
emit updateScale(axis,scaleOptions(axis, start, end,
boxMajors->text(), boxMinors->text()));
}
if (generalDialog->currentWidget()==(QWidget*)axes)
{
axis=axesList2->currentRow();
labels[axis] = boxLabel->text();
emit updateLabel(axis, boxLabel->text(),axisFont(axis));
emit updateTickLength(axis,boxMajorLength->text().toDouble(),
boxMinorLength->text().toDouble());
}
return true;
}
void FunctionDialog::acceptPolar() {
QString from = boxPolarFrom->text().toLower();
QString to = boxPolarTo->text().toLower();
QString points = boxPolarPoints->text().toLower();
double start, end;
try {
MyParser parser;
parser.SetExpr(from.toAscii().constData());
start = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(0, tr("MantidPlot - Start limit error"),
QString::fromStdString(e.GetMsg()));
boxPolarFrom->setFocus();
return;
}
try {
MyParser parser;
parser.SetExpr(to.toAscii().constData());
end = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(0, tr("MantidPlot - End limit error"),
QString::fromStdString(e.GetMsg()));
boxPolarTo->setFocus();
return;
}
if (start >= end) {
QMessageBox::critical(
0, tr("MantidPlot - Input error"),
tr("Please enter parameter limits that satisfy: from < end!"));
boxPolarTo->setFocus();
return;
}
QString rformula = boxPolarRadius->currentText();
QString tformula = boxPolarTheta->currentText();
bool error = false;
try {
MyParser parser;
double parameter = start;
;
parser.DefineVar((boxPolarParameter->text()).toAscii().constData(),
¶meter);
parser.SetExpr(rformula.toAscii().constData());
parser.Eval();
// cppcheck-suppress unreadVariable
parameter = end;
parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(0, tr("MantidPlot - Input function error"),
QString::fromStdString(e.GetMsg()));
boxPolarRadius->setFocus();
error = true;
}
try {
MyParser parser;
double parameter = start;
;
parser.DefineVar((boxPolarParameter->text()).toAscii().constData(),
¶meter);
parser.SetExpr(tformula.toAscii().constData());
parser.Eval();
// cppcheck-suppress unreadVariable
parameter = end;
parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(0, tr("MantidPlot - Input function error"),
QString::fromStdString(e.GetMsg()));
boxPolarTheta->setFocus();
error = true;
}
// Collecting all the information
int type = boxType->currentIndex();
QStringList formulas;
formulas += rformula;
formulas += tformula;
if (!error) {
d_app->updateFunctionLists(type, formulas);
if (!graph)
d_app->newFunctionPlot(formulas, start, end, boxPolarPoints->value(),
boxPolarParameter->text(), type);
else {
if (curveID >= 0)
graph->modifyFunctionCurve(curveID, type, formulas,
boxPolarParameter->text(), start, end,
boxPolarPoints->value());
else
graph->addFunction(formulas, start, end, boxPolarPoints->value(),
boxPolarParameter->text(), type);
}
}
}
bool FunctionCurve::loadData(int points, bool xLog10Scale)
{
if (!points)
points = dataSize();
double *X = (double *)malloc(points*sizeof(double));
if (!X){
QMessageBox::critical(0, QObject::tr("QtiPlot - Memory Allocation Error"),
QObject::tr("Not enough memory, operation aborted!"));
return false;
}
double *Y = (double *)malloc(points*sizeof(double));
if (!Y){
QMessageBox::critical(0, QObject::tr("QtiPlot - Memory Allocation Error"),
QObject::tr("Not enough memory, operation aborted!"));
free(X);
return false;
}
double step = (d_to - d_from)/(double)(points - 1.0);
if (d_function_type == Normal){
MyParser parser;
double x = d_from;
try {
parser.DefineVar(d_variable.ascii(), &x);
QMapIterator<QString, double> i(d_constants);
while (i.hasNext()){
i.next();
parser.DefineConst(i.key().ascii(), i.value());
}
parser.SetExpr(d_formulas[0].ascii());
int lastButOne = points - 1;
try {
double xl = x, xr;
double y = parser.EvalRemoveSingularity(&x, false);
bool wellDefinedFunction = true;
if (!gsl_finite(y)){// try to find a first well defined point (might help for some not really bad functions)
wellDefinedFunction = false;
for (int i = 0; i < lastButOne; i++){
xl = x;
x += step;
xr = x;
y = parser.Eval();
if (gsl_finite(y)){
wellDefinedFunction = true;
int iter = 0;
double x0 = x, y0 = y;
while(fabs(xr - xl)/step > 1e-15 && iter < points){
x = 0.5*(xl + xr);
y = parser.Eval();
if (gsl_finite(y)){
xr = x;
x0 = x;
y0 = y;
} else
xl = x;
iter++;
}
d_from = x0;
X[0] = x0;
Y[0] = y0;
step = (d_to - d_from)/(double)(lastButOne);
break;
}
}
if (!wellDefinedFunction){
QMessageBox::critical(0, QObject::tr("QtiPlot"),
QObject::tr("The function %1 is not defined in the specified interval!").arg(d_formulas[0]));
free(X); free(Y);
return false;
}
} else {
X[0] = d_from;
Y[0] = y;
}
} catch (MyParser::Pole) {}
ScaleEngine *sc_engine = 0;
if (plot())
sc_engine = (ScaleEngine *)plot()->axisScaleEngine(xAxis());
if (xLog10Scale || (d_from > 0 && d_to > 0 && sc_engine &&
sc_engine->type() == ScaleTransformation::Log10)){
step = log10(d_to/d_from)/(double)(points - 1);
for (int i = 1; i < lastButOne; i++ ){
x = d_from*pow(10, i*step);
X[i] = x;
try {
Y[i] = parser.EvalRemoveSingularity(&x, false);
} catch (MyParser::Pole){}
}
} else {
for (int i = 1; i < lastButOne; i++ ){
x += step;
X[i] = x;
try {
Y[i] = parser.EvalRemoveSingularity(&x, false);
} catch (MyParser::Pole){}
}
}
//the last point might be outside the interval, therefore we calculate it separately at its precise value
x = d_to;
X[lastButOne] = x;
try {
Y[lastButOne] = parser.EvalRemoveSingularity(&x, false);
} catch (MyParser::Pole){}
} catch(mu::ParserError &e) {}
} else if (d_function_type == Parametric || d_function_type == Polar) {
QStringList aux = d_formulas;
MyParser xparser;
MyParser yparser;
double par;
if (d_function_type == Polar) {
QString swap=aux[0];
aux[0]="("+swap+")*cos("+aux[1]+")";
aux[1]="("+swap+")*sin("+aux[1]+")";
}
try {
QMapIterator<QString, double> i(d_constants);
while (i.hasNext()){
i.next();
xparser.DefineConst(i.key().ascii(), i.value());
yparser.DefineConst(i.key().ascii(), i.value());
}
xparser.DefineVar(d_variable.ascii(), &par);
yparser.DefineVar(d_variable.ascii(), &par);
xparser.SetExpr(aux[0].ascii());
yparser.SetExpr(aux[1].ascii());
par = d_from;
for (int i = 0; i<points; i++ ){
X[i] = xparser.Eval();
Y[i] = yparser.Eval();
par += step;
}
} catch(mu::ParserError &) {}
}
if (curveType() == QwtPlotCurve::Yfx)
setData(X, Y, points);
else
setData(Y, X, points);
free(X); free(Y);
return true;
}
void SurfaceDialog::acceptFunction()
{
ApplicationWindow *app = static_cast<ApplicationWindow *>(this->parent());
QString Xfrom=boxXFrom->text().lower();
QString Xto=boxXTo->text().lower();
QString Yfrom=boxYFrom->text().lower();
QString Yto=boxYTo->text().lower();
QString Zfrom=boxZFrom->text().lower();
QString Zto=boxZTo->text().lower();
double fromX, toX, fromY,toY, fromZ,toZ;
try
{
MyParser parser;
parser.SetExpr(Xfrom.ascii());
fromX=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(app, tr("MantidPlot - X Start limit error"), QString::fromStdString(e.GetMsg()));
boxXFrom->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(Xto.ascii());
toX=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(app, tr("MantidPlot - X End limit error"), QString::fromStdString(e.GetMsg()));
boxXTo->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(Yfrom.ascii());
fromY=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(app, tr("MantidPlot - Y Start limit error"), QString::fromStdString(e.GetMsg()));
boxYFrom->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(Yto.ascii());
toY=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(app, tr("MantidPlot - Y End limit error"), QString::fromStdString(e.GetMsg()));
boxYTo->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(Zfrom.ascii());
fromZ=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(app, tr("MantidPlot - Z Start limit error"), QString::fromStdString(e.GetMsg()));
boxZFrom->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(Zto.ascii());
toZ=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(app, tr("MantidPlot - Z End limit error"), QString::fromStdString(e.GetMsg()));
boxZTo->setFocus();
return;
}
if (fromX >= toX || fromY >= toY || fromZ >= toZ)
{
QMessageBox::critical(app, tr("MantidPlot - Input error"),
tr("Please enter limits that satisfy: from < end!"));
boxXTo->setFocus();
return;
}
QString formula=boxFunction->currentText();
bool error=false;
try
{
MyParser parser;
double x=fromX;
double y=fromY;
parser.DefineVar("x", &x);
parser.DefineVar("y", &y);
parser.SetExpr(formula.ascii());
parser.Eval();
// cppcheck-suppress unreadVariable
x=toX;
// cppcheck-suppress unreadVariable
y=toY;
parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(0, tr("MantidPlot - Input function error"), QString::fromStdString(e.GetMsg()));
boxFunction->setFocus();
error=true;
}
if (!error){
QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
if (!d_graph){
app->plotSurface(boxFunction->currentText(),fromX, toX, fromY, toY, fromZ, toZ,
boxFuncColumns->value(), boxFuncRows->value());
} else
d_graph->addFunction(boxFunction->currentText(),fromX, toX, fromY, toY, fromZ, toZ,
boxFuncColumns->value(), boxFuncRows->value());
app->updateSurfaceFuncList(boxFunction->currentText());
QApplication::restoreOverrideCursor();
close();
}
}
void SurfaceDialog::acceptParametricSurface()
{
ApplicationWindow *app = static_cast<ApplicationWindow *>(this->parent());
MyParser parser;
double u = 1.0, v = 1.0;
parser.DefineVar("u", &u);
parser.DefineVar("v", &v);
int list_size = 15;
QString x_formula = boxX->text();
try {
parser.SetExpr(x_formula.ascii());
parser.Eval();
} catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - X Formula Error"), QString::fromStdString(e.GetMsg()));
boxX->setFocus();
return;
}
app->d_param_surface_func.remove(x_formula);
app->d_param_surface_func.push_front(x_formula);
while ((int)app->d_param_surface_func.size() > list_size)
app->d_param_surface_func.pop_back();
QString y_formula = boxY->text();
try {
parser.SetExpr(y_formula.ascii());
parser.Eval();
} catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - Y Formula Error"), QString::fromStdString(e.GetMsg()));
boxY->setFocus();
return;
}
app->d_param_surface_func.remove(y_formula);
app->d_param_surface_func.push_front(y_formula);
while ((int)app->d_param_surface_func.size() > list_size)
app->d_param_surface_func.pop_back();
QString z_formula = boxZ->text();
try {
parser.SetExpr(z_formula.ascii());
parser.Eval();
} catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - Z Formula Error"), QString::fromStdString(e.GetMsg()));
boxZ->setFocus();
return;
}
app->d_param_surface_func.remove(z_formula);
app->d_param_surface_func.push_front(z_formula);
while ((int)app->d_param_surface_func.size() > list_size)
app->d_param_surface_func.pop_back();
QString ufrom = boxUFrom->text().lower();
QString uto = boxUTo->text().lower();
QString vfrom = boxVFrom->text().lower();
QString vto = boxVTo->text().lower();
double ul, ur, vl, vr;
try{
parser.SetExpr(ufrom.ascii());
ul = parser.Eval();
}
catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - u start limit error"), QString::fromStdString(e.GetMsg()));
boxUFrom->setFocus();
return;
}
try{
parser.SetExpr(uto.ascii());
ur = parser.Eval();
}
catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - u end limit error"), QString::fromStdString(e.GetMsg()));
boxUTo->setFocus();
return;
}
try{
parser.SetExpr(vfrom.ascii());
vl = parser.Eval();
}
catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - v start limit error"), QString::fromStdString(e.GetMsg()));
boxVFrom->setFocus();
return;
}
try{
parser.SetExpr(vto.ascii());
vr = parser.Eval();
}
catch(mu::ParserError &e){
QMessageBox::critical(app, tr("MantidPlot - u end limit error"), QString::fromStdString(e.GetMsg()));
boxVTo->setFocus();
return;
}
QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
if (!d_graph)
app->plotParametricSurface(x_formula, y_formula, z_formula,
ul, ur, vl, vr, boxColumns->value(), boxRows->value(),
boxUPeriodic->isChecked(), boxVPeriodic->isChecked());
else
d_graph->addParametricSurface(x_formula, y_formula, z_formula,
ul, ur, vl, vr, boxColumns->value(), boxRows->value(),
boxUPeriodic->isChecked(), boxVPeriodic->isChecked());
QApplication::restoreOverrideCursor();
close();
}
bool Plot3DDialog::updatePlot()
{
if (!d_plot)
return false;
ApplicationWindow *app = (ApplicationWindow *)this->parent();
if (!app)
return false;
if (generalDialog->currentPage()==(QWidget*)bars){
d_plot->setBarRadius(boxBarsRad->text().toDouble());
d_plot->setBarStyle();
} else if (generalDialog->currentPage() == (QWidget*)points){
if (boxPointStyle->currentIndex() == 0) {
d_plot->setDotOptions(boxSize->text().toDouble(), boxSmooth->isChecked());
d_plot->setDotStyle();
} else if (boxPointStyle->currentIndex() == 1){
d_plot->setCrossOptions(boxCrossRad->text().toDouble(), boxCrossLinewidth->text().toDouble(),
boxCrossSmooth->isChecked(), boxBoxed->isChecked());
d_plot->setCrossStyle();
} else if (boxPointStyle->currentIndex() == 2) {
d_plot->setConeOptions(boxConesRad->text().toDouble(), boxQuality->value());
d_plot->setConeStyle();
}
app->custom3DActions(d_plot);
} else if (generalDialog->currentPage()==(QWidget*)title){
d_plot->setTitle(boxTitle->text().remove("\n"), btnTitleColor->color(), titleFont);
} else if (generalDialog->currentPage()==(QWidget*)colors){
d_plot->changeTransparency(boxTransparency->value()*0.01);
d_plot->setDataColors(btnFromColor->color(), btnToColor->color());
d_plot->setMeshColor(btnMesh->color());
d_plot->setAxesColor(btnAxes->color());
d_plot->setNumbersColor(btnNumbers->color());
d_plot->setLabelsColor(btnLabels->color());
d_plot->setBackgroundColor(btnBackground->color());
d_plot->setGridColor(btnGrid->color());
} else if (generalDialog->currentPage()==(QWidget*)general){
d_plot->showColorLegend(boxLegend->isChecked());
d_plot->setResolution(boxResolution->value());
d_plot->setMeshLineWidth(boxMeshLineWidth->value());
d_plot->setLabelsDistance(boxDistance->value());
d_plot->setNumbersFont(numbersFont);
d_plot->setZoom(boxZoom->value()*0.01);
d_plot->setScale(boxXScale->value()*0.01, boxYScale->value()*0.01, boxZScale->value()*0.01);
} else if (generalDialog->currentPage()==(QWidget*)scale){
int axis = axesList->currentRow();
QString from=boxFrom->text().lower();
QString to=boxTo->text().lower();
double start, end;
try {
MyParser parser;
parser.SetExpr(from.ascii());
start = parser.Eval();
} catch(mu::ParserError &e){
QMessageBox::critical(0,tr("QtiPlot - Start limit error"), QString::fromStdString(e.GetMsg()));
boxFrom->setFocus();
return false;
}
try {
MyParser parser;
parser.SetExpr(to.ascii());
end = parser.Eval();
} catch(mu::ParserError &e){
QMessageBox::critical(0,tr("QtiPlot - End limit error"), QString::fromStdString(e.GetMsg()));
boxTo->setFocus();
return false;
}
d_plot->updateScale(axis, scaleOptions(axis, start, end, boxMajors->text(), boxMinors->text()));
} else if (generalDialog->currentPage()==(QWidget*)axes){
int axis = axesList2->currentRow();
labels[axis] = boxLabel->text();
switch(axis)
{
case 0:
d_plot->setXAxisLabel(boxLabel->text().remove("\n"));
d_plot->setXAxisLabelFont(axisFont(axis));
d_plot->setXAxisTickLength(boxMajorLength->text().toDouble(), boxMinorLength->text().toDouble());
break;
case 1:
d_plot->setYAxisLabel(boxLabel->text().remove("\n"));
d_plot->setYAxisLabelFont(axisFont(axis));
d_plot->setYAxisTickLength(boxMajorLength->text().toDouble(), boxMinorLength->text().toDouble());
break;
case 2:
d_plot->setZAxisLabel(boxLabel->text().remove("\n"));
d_plot->setZAxisLabelFont(axisFont(axis));
d_plot->setZAxisTickLength(boxMajorLength->text().toDouble(), boxMinorLength->text().toDouble());
break;
}
}
d_plot->update();
app->modifiedProject(d_plot);
return true;
}
void FilterDialog::filter() {
double from = 0.0, to = 0.0;
try {
MyParser parser;
parser.SetExpr(boxStart->text().replace(",", ".").toAscii().constData());
from = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(this, tr("MantidPlot - Frequency input error"),
QString::fromStdString(e.GetMsg()));
boxStart->setFocus();
return;
}
if (from < 0) {
QMessageBox::critical(this, tr("MantidPlot - Frequency input error"),
tr("Please enter positive frequency values!"));
boxStart->setFocus();
return;
}
if (filter_type >= FFTFilter::BandPass) {
try {
MyParser parser;
parser.SetExpr(boxEnd->text().replace(",", ".").toAscii().constData());
to = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(this, tr("MantidPlot - High Frequency input error"),
QString::fromStdString(e.GetMsg()));
boxEnd->setFocus();
return;
}
if (to < 0) {
QMessageBox::critical(this, tr("MantidPlot - High Frequency input error"),
tr("Please enter positive frequency values!"));
boxEnd->setFocus();
return;
}
if (from >= to) {
QMessageBox::critical(
this, tr("MantidPlot - Frequency input error"),
tr("Please enter frequency limits that satisfy: Low < High !"));
boxEnd->setFocus();
return;
}
}
FFTFilter *f = new FFTFilter(dynamic_cast<ApplicationWindow *>(parent()),
graph, boxName->currentText(), filter_type);
if (filter_type == FFTFilter::BandPass) {
f->setBand(from, to);
f->enableOffset(boxOffset->isChecked());
} else if (filter_type == FFTFilter::BandBlock) {
f->setBand(from, to);
f->enableOffset(!boxOffset->isChecked());
} else
f->setCutoff(from);
f->setColor(boxColor->currentIndex());
f->run();
delete f;
}
void FitDialog::accept()
{
QString curve = boxCurve->currentText();
QStringList curvesList = d_graph->curvesList();
if (curvesList.contains(curve) <= 0)
{
QMessageBox::critical(this,tr("Warning"),
tr("The curve <b> %1 </b> doesn't exist anymore! Operation aborted!").arg(curve));
boxCurve->clear();
boxCurve->addItems(curvesList);
return;
}
if (!validInitialValues())
return;
QString from=boxFrom->text().toLower();
QString to=boxTo->text().toLower();
QString tolerance=boxTolerance->text().toLower();
double start, end, eps;
try
{
MyParser parser;
parser.SetExpr(CONFS(from).toAscii().constData());
start=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(this, tr("Start limit error"),QString::fromStdString(e.GetMsg()));
boxFrom->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(CONFS(to).toAscii().constData());
end=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(this, tr("End limit error"),QString::fromStdString(e.GetMsg()));
boxTo->setFocus();
return;
}
if (start>=end)
{
QMessageBox::critical(0, tr("Input error"),
tr("Please enter x limits that satisfy: from < end!"));
boxTo->setFocus();
return;
}
try
{
MyParser parser;
parser.SetExpr(CONFS(tolerance).toAscii().constData());
eps=parser.Eval();
}
catch(mu::ParserError &e)
{
QMessageBox::critical(0, tr("Tolerance input error"),QString::fromStdString(e.GetMsg()));
boxTolerance->setFocus();
return;
}
if (eps<0 || eps>=1)
{
QMessageBox::critical(0, tr("Tolerance input error"),
tr("The tolerance value must be positive and less than 1!"));
boxTolerance->setFocus();
return;
}
int i, n=0, rows=boxParams->rowCount();
if (!boxParams->isColumnHidden(2))
{
for (i=0;i<rows;i++)
{//count the non-constant parameters
QCheckBox *cb = (QCheckBox*)boxParams->cellWidget(i, 2);
if (!cb->isChecked())
n++;
}
}
else
n=rows;
QStringList parameters;
double *paramsInit = new double[n];
QString formula;
// recursively define variables for user functions used in formula
bool found_uf;
do {
found_uf = false;
for (i=0; i<d_user_function_names.count(); i++)
if (boxFunction->text().contains(d_user_function_names[i])) {
QStringList l = d_user_functions[i].split("=");
formula += QString("%1=%2\n")
.arg(d_user_function_names[i])
.arg(l[1]);
found_uf = true;
}
} while (found_uf);
formula += boxFunction->text();
// define variables for builtin functions used in formula
for (i=0; i<d_built_in_function_names.count(); i++)
if (formula.contains(d_built_in_function_names[i]))
formula.prepend(QString("%1=%2\n")
.arg(d_built_in_function_names[i])
.arg(d_built_in_functions[i]));
if (!boxParams->isColumnHidden(2))
{
int j = 0;
for (i=0;i<rows;i++)
{
QCheckBox *cb = (QCheckBox*)boxParams->cellWidget(i, 2);
if (!cb->isChecked())
{
paramsInit[j] = QLocale().toDouble(boxParams->item(i,1)->text());
parameters << boxParams->item(i,0)->text();
j++;
}
else
formula.prepend(QString("%1=%2\n")
.arg(boxParams->item(i,0)->text())
.arg(CONFS(boxParams->item(i,1)->text())));
}
}
else
{
for (i=0;i<n;i++)
{
paramsInit[i] = QLocale().toDouble(boxParams->item(i,1)->text());
parameters << boxParams->item(i,0)->text();
}
}
ApplicationWindow *app = (ApplicationWindow *)this->parent();
if (d_fitter)
{
delete d_fitter;
d_fitter = 0;
}
if (boxUseBuiltIn->isChecked() && categoryBox->currentRow() == 1)
fitBuiltInFunction(funcBox->currentItem()->text(), paramsInit);
else if (boxUseBuiltIn->isChecked() && categoryBox->currentRow() == 3)
{
d_fitter = new PluginFit(app, d_graph);
if (!((PluginFit*)d_fitter)->load(d_plugin_files_list[funcBox->currentRow()])){
d_fitter = 0;
return;}
d_fitter->setInitialGuesses(paramsInit);
}
else
{
d_fitter = new NonLinearFit(app, d_graph);
((NonLinearFit*)d_fitter)->setParametersList(parameters);
((NonLinearFit*)d_fitter)->setFormula(formula);
d_fitter->setInitialGuesses(paramsInit);
}
delete[] paramsInit;
if (!d_fitter->setDataFromCurve(curve, start, end) ||
!d_fitter->setYErrorSource ((Fit::ErrorSource)boxYErrorSource->currentIndex(),
tableNamesBox->currentText()+"_"+colNamesBox->currentText()))
{
delete d_fitter;
d_fitter = 0;
return;
}
d_fitter->setTolerance (eps);
d_fitter->setAlgorithm((Fit::Algorithm)boxAlgorithm->currentIndex());
d_fitter->setColor(boxColor->currentIndex());
d_fitter->generateFunction(generatePointsBtn->isChecked(), generatePointsBox->value());
d_fitter->setMaximumIterations(boxPoints->value());
d_fitter->scaleErrors(scaleErrorsBox->isChecked());
if (d_fitter->name() == tr("MultiPeak") && ((MultiPeakFit *)d_fitter)->peaks() > 1)
{
((MultiPeakFit *)d_fitter)->enablePeakCurves(app->generatePeakCurves);
((MultiPeakFit *)d_fitter)->setPeakCurvesColor(app->peakCurvesColor);
}
d_fitter->fit();
double *res = d_fitter->results();
if (!boxParams->isColumnHidden(2))
{
int j = 0;
for (i=0;i<rows;i++)
{
QCheckBox *cb = (QCheckBox*)boxParams->cellWidget(i, 2);
if (!cb->isChecked())
boxParams->item(i, 1)->setText(QLocale().toString(res[j++], 'g', boxPrecision->value()));
}
}
else
{
for (i=0;i<rows;i++)
boxParams->item(i, 1)->setText(QLocale().toString(res[i], 'g', boxPrecision->value()));
}
}
bool Plot3DDialog::updatePlot() {
if (!d_plot)
return false;
ApplicationWindow *app = static_cast<ApplicationWindow *>(this->parent());
if (!app)
return false;
if (generalDialog->currentWidget() == static_cast<QWidget *>(bars)) {
d_plot->setBarRadius(boxBarsRad->text().toDouble());
d_plot->setBarStyle();
} else if (generalDialog->currentWidget() == static_cast<QWidget *>(points)) {
if (boxPointStyle->currentIndex() == 0) {
d_plot->setDotOptions(boxSize->text().toDouble(), boxSmooth->isChecked());
d_plot->setDotStyle();
} else if (boxPointStyle->currentIndex() == 1) {
d_plot->setCrossOptions(
boxCrossRad->text().toDouble(), boxCrossLinewidth->text().toDouble(),
boxCrossSmooth->isChecked(), boxBoxed->isChecked());
d_plot->setCrossStyle();
} else if (boxPointStyle->currentIndex() == 2) {
d_plot->setConeOptions(boxConesRad->text().toDouble(),
boxQuality->value());
d_plot->setConeStyle();
}
app->custom3DActions(d_plot);
} else if (generalDialog->currentWidget() == static_cast<QWidget *>(title)) {
d_plot->setTitle(boxTitle->toPlainText().remove("\n"),
btnTitleColor->color(), titleFont);
} else if (generalDialog->currentWidget() == static_cast<QWidget *>(colors)) {
d_plot->changeTransparency(boxTransparency->value() * 0.01);
d_plot->setDataColors(btnFromColor->color(), btnToColor->color());
d_plot->setMeshColor(btnMesh->color());
d_plot->setAxesColor(btnAxes->color());
d_plot->setNumbersColor(btnNumbers->color());
d_plot->setLabelsColor(btnLabels->color());
d_plot->setBackgroundColor(btnBackground->color());
d_plot->setGridColor(btnGrid->color());
} else if (generalDialog->currentWidget() ==
static_cast<QWidget *>(general)) {
d_plot->showColorLegend(boxLegend->isChecked());
d_plot->setResolution(boxResolution->value());
d_plot->setMeshLineWidth(boxMeshLineWidth->value());
d_plot->setLabelsDistance(boxDistance->value());
d_plot->setNumbersFont(numbersFont);
d_plot->setZoom(zoom * boxZoom->value() * 0.01);
d_plot->setScale(xScale * boxXScale->value() * 0.01,
yScale * boxYScale->value() * 0.01,
zScale * boxZScale->value() * 0.01);
} else if (generalDialog->currentWidget() == static_cast<QWidget *>(scale)) {
int axis = axesList->currentRow();
QString from = boxFrom->text().toLower();
QString to = boxTo->text().toLower();
double start, end;
try {
MyParser parser;
parser.SetExpr(from.toAscii().constData());
start = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(0, tr("MantidPlot - Start limit error"),
QString::fromStdString(e.GetMsg()));
boxFrom->setFocus();
return false;
}
try {
MyParser parser;
parser.SetExpr(to.toAscii().constData());
end = parser.Eval();
} catch (mu::ParserError &e) {
QMessageBox::critical(0, tr("MantidPlot - End limit error"),
QString::fromStdString(e.GetMsg()));
boxTo->setFocus();
return false;
}
/*double xsc = d_plot->xScale();
double ysc = d_plot->yScale();
double zsc = d_plot->zScale();
if (axis == 2)
{
double start0 = scales[0].toDouble();
double end0 = scales[1].toDouble();
zsc *= (end0 - start0)/(end - start);
QMessageBox::information(this,"OK","OK");
}*/
d_plot->updateScale(axis, scaleOptions(axis, start, end, boxMajors->text(),
boxMinors->text()));
// d_plot->setScale(xsc,ysc,zsc*0.1);
} else if (generalDialog->currentWidget() == static_cast<QWidget *>(axes)) {
int axis = axesList2->currentRow();
labels[axis] = boxLabel->toPlainText();
switch (axis) {
case 0:
d_plot->setXAxisLabel(boxLabel->toPlainText().remove("\n"));
d_plot->setXAxisLabelFont(axisFont(axis));
d_plot->setXAxisTickLength(boxMajorLength->text().toDouble(),
boxMinorLength->text().toDouble());
break;
case 1:
d_plot->setYAxisLabel(boxLabel->toPlainText().remove("\n"));
d_plot->setYAxisLabelFont(axisFont(axis));
d_plot->setYAxisTickLength(boxMajorLength->text().toDouble(),
boxMinorLength->text().toDouble());
break;
case 2:
d_plot->setZAxisLabel(boxLabel->toPlainText().remove("\n"));
d_plot->setZAxisLabelFont(axisFont(axis));
d_plot->setZAxisTickLength(boxMajorLength->text().toDouble(),
boxMinorLength->text().toDouble());
break;
}
}
d_plot->update();
app->modifiedProject(d_plot);
return true;
}
void IntDialog::accept()
{
QString curveName = boxName->currentText();
QwtPlotCurve *c = graph->curve(curveName);
QStringList curvesList = graph->analysableCurvesList();
if (!c || !curvesList.contains(curveName))
{
QMessageBox::critical((ApplicationWindow *)parent(), tr("SciDAVis") +" - "+ tr("Warning"),
tr("The curve <b> %1 </b> doesn't exist anymore! Operation aborted!").arg(curveName));
boxName->clear();
boxName->insertStringList(curvesList);
return;
}
double start = 0, stop = 0;
double minx = c->minXValue();
double maxx = c->maxXValue();
// Check the Xmin
QString from = boxStart->text().toLower();
if(from=="min")
{
boxStart->setText(QString::number(minx));
return;
}
else if(from=="max")
{
boxStart->setText(QString::number(maxx));
return;
}
else
{
try
{
MyParser parser;
parser.SetExpr((boxStart->text()).toAscii().constData());
start=parser.Eval();
if(start<minx)
{
QMessageBox::warning((ApplicationWindow *)parent(), tr("Input error"),
tr("Please give a number larger or equal to the minimum value of X, for the lower limit.\n If you do not know that value, type min in the box."));
boxStart->clear();
boxStart->setFocus();
return;
}
if(start > maxx)
{
QMessageBox::warning((ApplicationWindow *)parent(), tr("Input error"),
tr("Please give a number smaller or equal to the maximum value of X, for the lower limit.\n If you do not know that value, type max in the box."));
boxStart->clear();
boxStart->setFocus();
return;
}
}
catch(mu::ParserError &e)
{
QMessageBox::critical((ApplicationWindow *)parent(),tr("Start limit error"),QString::fromStdString(e.GetMsg()));
boxStart->clear();
boxStart->setFocus();
return;
}
}
// Check Xmax
QString end=boxEnd->text().toLower();
if(end=="min")
{
boxEnd->setText(QString::number(minx));
return;
}
else if(end=="max")
{
boxEnd->setText(QString::number(maxx));
return;
}
else
{
try
{
MyParser parser;
parser.SetExpr((boxEnd->text()).toAscii().constData());
stop = parser.Eval();
if(stop > maxx)
{
//FIXME: I don't understand why this doesn't work for FunctionCurves!!(Ion)
/*QMessageBox::warning((ApplicationWindow *)parent(), tr("Input error"),
tr("Please give a number smaller or equal to the maximum value of X, for the upper limit.\n If you do not know that value, type max in the box."));
boxEnd->clear();
boxEnd->setFocus();
return;
*/
boxEnd->setText(QString::number(maxx));
}
if(stop < minx)
{
QMessageBox::warning((ApplicationWindow *)parent(), tr("Input error"),
tr("Please give a number larger or equal to the minimum value of X, for the upper limit.\n If you do not know that value, type min in the box."));
boxEnd->clear();
boxEnd->setFocus();
return;
}
}
catch(mu::ParserError &e)
{
QMessageBox::critical((ApplicationWindow *)parent(), tr("End limit error"),QString::fromStdString(e.GetMsg()));
boxEnd->clear();
boxEnd->setFocus();
return;
}
}
Integration *i = new Integration((ApplicationWindow *)this->parent(), graph, curveName,
boxStart->text().toDouble(), boxEnd->text().toDouble());
i->setMethod((Integration::InterpolationMethod)boxMethod->currentIndex());
i->run();
delete i;
}