1 // Copyright (C) 2005 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License.
9 // This library is distributed in the hope that it will be useful
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
18 // See http://www.salome-platform.org/
21 #include "StdMeshersGUI_DistrPreview.h"
22 #include "CASCatch.hxx"
24 #include <Expr_NamedUnknown.hxx>
25 #include <Expr_GeneralExpression.hxx>
27 StdMeshersGUI_DistrPreview::StdMeshersGUI_DistrPreview( QWidget* p, StdMeshers::StdMeshers_NumberOfSegments_ptr h )
33 myConv( CUT_NEGATIVE ),
37 myHypo = StdMeshers::StdMeshers_NumberOfSegments::_duplicate( h );
38 myVars.ChangeValue( 1 ) = new Expr_NamedUnknown( "t" );
39 myDensity = insertCurve( QString() );
40 myDistr = insertCurve( QString() );
41 myMsg = insertMarker( new QwtPlotMarker( this ) );
42 setMarkerPos( myMsg, 0.5, 0.5 );
43 setMarkerLabelPen( myMsg, QPen( Qt::red, 1 ) );
44 QFont f = markerFont( myMsg );
47 setMarkerFont( myMsg, f );
48 setCurvePen( myDensity, QPen( Qt::red, 1 ) );
51 setCurvePen( myDistr, QPen( dc, 1 ) );
52 setCurveSymbol( myDistr, QwtSymbol( QwtSymbol::XCross, QBrush( dc ), QPen( dc ), QSize( 5, 5 ) ) );
53 setAutoLegend( true );
55 setLegendPos( Qwt::Bottom );
56 setCurveTitle( myDensity, tr( "SMESH_DENSITY_FUNC" ) );
57 setCurveTitle( myDistr, tr( "SMESH_DISTR" ) );
60 StdMeshersGUI_DistrPreview::~StdMeshersGUI_DistrPreview()
64 bool StdMeshersGUI_DistrPreview::isTableFunc() const
69 void StdMeshersGUI_DistrPreview::tableFunc( SMESH::double_array& f ) const
74 QString StdMeshersGUI_DistrPreview::function() const
79 int StdMeshersGUI_DistrPreview::nbSeg() const
84 int StdMeshersGUI_DistrPreview::pointsCount() const
89 void StdMeshersGUI_DistrPreview::setConversion( Conversion conv, const bool upd )
96 bool StdMeshersGUI_DistrPreview::setParams( const QString& func, const int nbSeg, const int points, const bool upd )
99 myTableFunc = SMESH::double_array();
100 myFunction = func.isEmpty() ? "0" : func;
101 myPoints = points>0 ? points : 2;
102 myNbSeg = nbSeg>0 ? nbSeg : 1;
103 bool res = init( func );
109 bool StdMeshersGUI_DistrPreview::setParams( const SMESH::double_array& f, const int nbSeg, const bool upd )
113 if( myTableFunc.length()%2==1 )
114 myTableFunc.length( myTableFunc.length()-1 );
117 myPoints = myTableFunc.length()/2;
118 myNbSeg = nbSeg>0 ? nbSeg : 1;
123 return myTableFunc.length()>0;
126 bool StdMeshersGUI_DistrPreview::createTable( SMESH::double_array& func )
128 if( myExpr.IsNull() )
134 const double xmin = 0.0, xmax = 1.0;
136 double d = (xmax-xmin)/double(myPoints-1);
137 func.length( 2*myPoints );
139 for( int i=0, j=0; i<myPoints; j++ )
142 double t = xmin + d*j, f = funcValue( t, ok );
152 func.length( func.length()-2*err );
156 void StdMeshersGUI_DistrPreview::update()
158 SMESH::double_array graph, distr;
165 myIsDone = createTable( graph );
167 if( graph.length()>=2 )
169 StdMeshers::StdMeshers_NumberOfSegments_var h =
170 StdMeshers::StdMeshers_NumberOfSegments::_narrow( myHypo );
172 if( !CORBA::is_nil( h.in() ) )
174 SMESH::double_array* arr = 0;
176 arr = h->BuildDistributionTab( myTableFunc, myNbSeg, ( int )myConv );
178 arr = h->BuildDistributionExpr( myFunction.latin1(), myNbSeg, ( int )myConv );
187 bool correct = graph.length()>=2 && distr.length()>=2;
194 setMarkerLabel( myMsg, QString() );
196 int size = graph.length()/2;
197 double* x = new double[size], *y = new double[size];
198 double min_x, max_x, min_y, max_y;
199 for( int i=0; i<size; i++ )
203 if( !convert( y[i] ) )
205 min_x = 0.0; max_x = 1.0; min_y = 0.0; max_y = 1.0;
206 delete[] x; delete[] y;
211 if( i==0 || y[i]<min_y )
213 if( i==0 || y[i]>max_y )
215 if( i==0 || x[i]<min_x )
217 if( i==0 || x[i]>max_x )
221 setAxisScale( curveXAxis( myDensity ), min_x, max_x );
222 setAxisScale( curveYAxis( myDensity ), min( 0.0, min_y ), max( 0.0, max_y ) );
223 setCurveData( myDensity, x, y, size );
230 size = distr.length();
231 x = new double[size];
232 y = new double[size];
233 for( int i=0; i<size; i++ )
238 setCurveData( myDistr, x, y, size );
247 CASCatch_CATCH(Standard_Failure)
249 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
253 void StdMeshersGUI_DistrPreview::showError()
255 setAxisScale( curveXAxis( myDensity ), 0.0, 1.0 );
256 setAxisScale( curveYAxis( myDensity ), 0.0, 1.0 );
257 setCurveData( myDensity, 0, 0, 0 );
258 setCurveData( myDistr, 0, 0, 0 );
259 setMarkerLabel( myMsg, tr( "SMESH_INVALID_FUNCTION" ) );
263 bool isCorrectArg( const Handle( Expr_GeneralExpression )& expr )
265 Handle( Expr_NamedUnknown ) sub = Handle( Expr_NamedUnknown )::DownCast( expr );
267 return sub->GetName()=="t";
270 for( int i=1, n=expr->NbSubExpressions(); i<=n && res; i++ )
272 Handle( Expr_GeneralExpression ) sub = expr->SubExpression( i );
273 Handle( Expr_NamedUnknown ) name = Handle( Expr_NamedUnknown )::DownCast( sub );
276 if( name->GetName()!="t" )
280 res = isCorrectArg( sub );
285 bool StdMeshersGUI_DistrPreview::init( const QString& str )
287 bool parsed_ok = true;
290 myExpr = ExprIntrp_GenExp::Create();
291 myExpr->Process( ( Standard_CString ) str.latin1() );
293 CASCatch_CATCH(Standard_Failure)
295 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
299 bool syntax = false, args = false;
300 if( parsed_ok && myExpr->IsDone() )
303 args = isCorrectArg( myExpr->Expression() );
306 bool res = parsed_ok && syntax && args;
312 double StdMeshersGUI_DistrPreview::funcValue( const double t, bool& ok )
314 if( myExpr.IsNull() )
317 myValues.ChangeValue( 1 ) = t;
320 double res = calc( ok );
325 double StdMeshersGUI_DistrPreview::calc( bool& ok )
331 res = myExpr->Expression()->Evaluate( myVars, myValues );
333 CASCatch_CATCH(Standard_Failure) {
334 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
342 bool StdMeshersGUI_DistrPreview::isDone() const
347 bool StdMeshersGUI_DistrPreview::convert( double& v ) const
356 // in StdMeshers_NumberOfSegments.cc
357 // const double PRECISION = 1e-7;
362 CASCatch_CATCH(Standard_Failure)
364 Handle(Standard_Failure) aFail = Standard_Failure::Caught();