1 // Copyright (c) 1997-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
4 // This file is part of Open CASCADE Technology software library.
6 // This library is free software; you can redistribute it and / or modify it
7 // under the terms of the GNU Lesser General Public version 2.1 as published
8 // by the Free Software Foundation, with special exception defined in the file
9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
10 // distribution for complete text of the license and disclaimer of any warranty.
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
18 #include <ExprIntrp_yaccintrf.hxx>
19 #include <ExprIntrp_yaccanal.hxx>
20 #include <Expr_GeneralExpression.hxx>
21 #include <Expr_NamedExpression.hxx>
22 #include <Expr_NamedUnknown.hxx>
23 #include <Expr_NamedConstant.hxx>
24 #include <Expr_NamedFunction.hxx>
25 #include <Expr_NumericValue.hxx>
26 #include <Expr_UnaryFunction.hxx>
27 #include <Expr_BinaryFunction.hxx>
28 #include <Expr_PolyFunction.hxx>
29 #include <Expr_Exponentiate.hxx>
30 #include <Expr_Absolute.hxx>
31 #include <Expr_ArcCosine.hxx>
32 #include <Expr_ArcSine.hxx>
33 #include <Expr_ArcTangent.hxx>
34 #include <Expr_ArgCosh.hxx>
35 #include <Expr_ArgSinh.hxx>
36 #include <Expr_ArgTanh.hxx>
37 #include <Expr_Cosh.hxx>
38 #include <Expr_Cosine.hxx>
39 #include <Expr_Exponential.hxx>
40 #include <Expr_LogOf10.hxx>
41 #include <Expr_LogOfe.hxx>
42 #include <Expr_Sign.hxx>
43 #include <Expr_Sine.hxx>
44 #include <Expr_Sinh.hxx>
45 #include <Expr_SquareRoot.hxx>
46 #include <Expr_Tangent.hxx>
47 #include <Expr_Tanh.hxx>
48 #include <Expr_Equal.hxx>
49 #include <Expr_SystemRelation.hxx>
50 #include <Expr_UnknownIterator.hxx>
51 #include <Expr_FunctionDerivative.hxx>
52 #include <Expr.hxx> // ATTENTION POUR PROTECTION BUG STACK
54 #include <Expr_SequenceOfGeneralExpression.hxx>
55 #include <Expr_Operators.hxx>
56 #include <ExprIntrp_SyntaxError.hxx>
57 #include <Expr_Array1OfNamedUnknown.hxx>
58 #include <Expr_Array1OfGeneralExpression.hxx>
65 static char ExprIntrp_assname[30];
66 TCollection_AsciiString ExprIntrp_funcdefname;
67 static Standard_Integer ExprIntrp_nbargs;
68 static Standard_Integer ExprIntrp_nbdiff;
71 extern "C" void ExprIntrp_StartFunction()
74 ExprIntrp_GetResult(funcname);
75 TCollection_AsciiString name(funcname);
76 ExprIntrp_Recept.PushName(name);
80 extern "C" void ExprIntrp_StartDerivate()
83 ExprIntrp_GetResult(funcname);
84 TCollection_AsciiString name(funcname);
85 ExprIntrp_Recept.PushName(name);
88 extern "C" void ExprIntrp_EndDerivate()
91 degree = ExprIntrp_GetDegree();
92 ExprIntrp_Recept.PushValue(degree);
95 extern "C" void ExprIntrp_Derivation()
97 ExprIntrp_Recept.PushValue(1);
99 ExprIntrp_GetResult(name);
100 TCollection_AsciiString thename(name);
101 Handle(Expr_NamedExpression) namexp = ExprIntrp_Recept.GetNamed(thename);
102 if (namexp.IsNull()) {
103 namexp = new Expr_NamedUnknown(thename);
105 if (!namexp->IsKind(STANDARD_TYPE(Expr_NamedUnknown))) {
106 ExprIntrp_SyntaxError::Raise();
108 ExprIntrp_Recept.Push(namexp);
111 extern "C" void ExprIntrp_DerivationValue()
114 ExprIntrp_GetResult(num);
115 Standard_Integer degree = ExprIntrp_Recept.PopValue();
117 ExprIntrp_Recept.PushValue(degree);
120 extern "C" void ExprIntrp_EndDerivation()
122 Standard_Integer degree = ExprIntrp_Recept.PopValue();
123 Handle(Expr_NamedUnknown) var = Handle(Expr_NamedUnknown)::DownCast(ExprIntrp_Recept.Pop());
124 Handle(Expr_GeneralExpression) exp = ExprIntrp_Recept.Pop();
125 exp = exp->NDerivative(var,degree);
126 ExprIntrp_Recept.Push(exp);
129 extern "C" void ExprIntrp_StartDifferential()
131 ExprIntrp_StartDerivate();
132 ExprIntrp_nbdiff = 0;
135 extern "C" void ExprIntrp_DiffDegreeVar()
138 ExprIntrp_GetResult(name);
140 if (strncasecmp(name,"X",1)) {
142 if ( name[ 0 ] != 'X' && name[ 0 ] != 'x' ) {
144 ExprIntrp_SyntaxError::Raise();
148 Standard_Integer rank = atoi(s);
149 ExprIntrp_Recept.PushValue(rank);
153 extern "C" void ExprIntrp_DiffVar()
155 ExprIntrp_Recept.PushValue(1);
156 ExprIntrp_DiffDegreeVar();
159 extern "C" void ExprIntrp_DiffDegree()
162 ExprIntrp_GetResult(name);
163 Standard_Integer deg = atoi(name);
164 ExprIntrp_Recept.PushValue(deg);
167 extern "C" void ExprIntrp_VerDiffDegree()
170 ExprIntrp_GetResult(name);
171 Standard_Integer deg = atoi(name);
172 Standard_Integer thedeg = ExprIntrp_Recept.PopValue();
174 ExprIntrp_SyntaxError::Raise();
176 ExprIntrp_Recept.PushValue(deg);
179 extern "C" void ExprIntrp_EndDifferential()
181 TCollection_AsciiString name = ExprIntrp_Recept.PopName();
182 Handle(Expr_GeneralFunction) thefunc = ExprIntrp_Recept.GetFunction(name);
183 if (thefunc.IsNull()) {
184 ExprIntrp_SyntaxError::Raise();
186 Standard_Integer rank,degree;
187 Handle(Expr_NamedUnknown) thediff;
188 Standard_Integer nbvars = thefunc->NbOfVariables();
190 for (Standard_Integer i=1; i<= ExprIntrp_nbdiff; i++) {
191 rank = ExprIntrp_Recept.PopValue();
192 degree = ExprIntrp_Recept.PopValue();
193 if ((rank > nbvars) || (rank < 1)) {
194 ExprIntrp_SyntaxError::Raise();
196 thediff = thefunc->Variable(rank);
197 thefunc = new Expr_FunctionDerivative(thefunc,thediff,degree);
199 ExprIntrp_Recept.PushFunction(thefunc);
202 extern "C" void ExprIntrp_EndDiffFunction()
204 Handle(Expr_GeneralFunction) thefunc = ExprIntrp_Recept.PopFunction();
205 if (thefunc.IsNull()) {
206 ExprIntrp_SyntaxError::Raise();
208 Standard_Integer nbargs = thefunc->NbOfVariables();
210 Handle(Expr_GeneralExpression) op = ExprIntrp_Recept.Pop();
211 Handle(Expr_UnaryFunction) res =
212 new Expr_UnaryFunction(thefunc,op);
213 ExprIntrp_Recept.Push(res);
215 else if (nbargs == 2) {
216 Handle(Expr_GeneralExpression) arg2 = ExprIntrp_Recept.Pop();
217 Handle(Expr_GeneralExpression) arg1 = ExprIntrp_Recept.Pop();
219 ExprIntrp_SyntaxError::Raise();
221 Handle(Expr_BinaryFunction) res =
222 new Expr_BinaryFunction(thefunc,arg1,arg2);
223 ExprIntrp_Recept.Push(res);
226 Expr_Array1OfGeneralExpression tabarg(1,nbargs);
227 Handle(Expr_GeneralExpression) arg;
228 for (Standard_Integer i = 1; i<= nbargs; i++) {
229 arg = ExprIntrp_Recept.Pop();
231 ExprIntrp_SyntaxError::Raise();
233 tabarg(nbargs-i+1) = arg;
235 Handle(Expr_PolyFunction) res =
236 new Expr_PolyFunction(thefunc,tabarg);
237 ExprIntrp_Recept.Push(res);
241 static Handle(Expr_GeneralExpression) ExprIntrp_StandardFunction(const TCollection_AsciiString& name, const Handle(Expr_GeneralExpression)& op)
243 // return standard functions equivalent corresponding to <name>
244 // with given operand <op> if exists. Returns null value if not.
245 // <name> is not case sensitive
247 Handle(Expr_GeneralExpression) res;
248 if ((name == "abs") || (name == "Abs")) {
249 res = new Expr_Absolute(op);
251 else if ((name == "acos") || (name == "ACos")) {
252 res = new Expr_ArcCosine(op);
254 else if ((name == "asin") || (name == "ASin")) {
255 res = new Expr_ArcSine(op);
257 else if ((name == "atan") || (name == "ATan")) {
258 res = new Expr_ArcTangent(op);
260 else if ((name == "acosh") || (name == "ACosh")) {
261 res = new Expr_ArgCosh(op);
263 else if ((name == "asinh") || (name == "ASinh")) {
264 res = new Expr_ArgSinh(op);
266 else if ((name == "atanh") || (name == "ATanh")) {
267 res = new Expr_ArgTanh(op);
269 else if ((name == "cosh") || (name == "Cosh")) {
270 res = new Expr_Cosh(op);
272 else if ((name == "cos") || (name == "Cos")) {
273 res = new Expr_Cosine(op);
275 else if ((name == "exp") || (name == "Exp")) {
276 res = new Expr_Exponential(op);
278 else if (name == "log") {
279 res = new Expr_LogOf10(op);
281 else if ((name == "Log") || (name == "Ln")) {
282 res = new Expr_LogOfe(op);
284 else if ((name == "sign") || (name == "Sign")) {
285 res = new Expr_Sign(op);
287 else if ((name == "sin") || (name == "Sin")) {
288 res = new Expr_Sine(op);
290 else if ((name == "sinh") || (name == "Sinh")) {
291 res = new Expr_Sinh(op);
293 else if ((name == "sqrt") || (name == "Sqrt")) {
294 res = new Expr_SquareRoot(op);
296 else if ((name == "tan") || (name == "Tan")) {
297 res = new Expr_Tangent(op);
299 else if ((name == "tanh") || (name == "Tanh")) {
300 res = new Expr_Tanh(op);
306 extern "C" void ExprIntrp_EndDerFunction()
308 TCollection_AsciiString name = ExprIntrp_Recept.PopName();
309 Handle(Expr_GeneralExpression) op = ExprIntrp_Recept.Pop();
310 Handle(Expr_GeneralExpression) resstand = ExprIntrp_StandardFunction(name,op);
312 if (!resstand.IsNull()) {
313 Handle(Expr_NamedUnknown) var;
314 Expr_UnknownIterator rit(resstand);
317 ExprIntrp_SyntaxError::Raise();
321 if (var->IsAssigned()) {
328 ExprIntrp_SyntaxError::Raise();
331 Handle(Expr_GeneralExpression) res = resstand->NDerivative(var,ExprIntrp_Recept.PopValue());
332 ExprIntrp_Recept.Push(res);
336 Handle(Expr_NamedFunction) thefunc = ExprIntrp_Recept.GetFunction(name);
337 if (thefunc.IsNull()) {
338 ExprIntrp_SyntaxError::Raise();
340 Standard_Integer nbargs = thefunc->NbOfVariables();
342 ExprIntrp_SyntaxError::Raise();
344 Handle(Expr_NamedUnknown) var = thefunc->Variable(1);
345 Handle(Expr_FunctionDerivative) thefuncder =
346 new Expr_FunctionDerivative(thefunc,var,ExprIntrp_Recept.PopValue());
347 Handle(Expr_UnaryFunction) res =
348 new Expr_UnaryFunction(thefuncder,op);
349 ExprIntrp_Recept.Push(res);
353 extern "C" void ExprIntrp_EndFunction()
355 TCollection_AsciiString name = ExprIntrp_Recept.PopName();
356 Handle(Expr_GeneralExpression) op = ExprIntrp_Recept.Pop();
358 Handle(Expr_GeneralExpression) resstand = ExprIntrp_StandardFunction(name,op);
359 if (!resstand.IsNull()) {
360 ExprIntrp_Recept.Push(resstand->ShallowSimplified());
363 Handle(Expr_NamedFunction) thefunc = ExprIntrp_Recept.GetFunction(name);
364 if (thefunc.IsNull()) {
365 ExprIntrp_SyntaxError::Raise();
367 Standard_Integer nbargs = thefunc->NbOfVariables();
369 Handle(Expr_UnaryFunction) res =
370 new Expr_UnaryFunction(thefunc,op);
371 ExprIntrp_Recept.Push(res);
373 else if (nbargs == 2) {
374 Handle(Expr_GeneralExpression) arg1 = ExprIntrp_Recept.Pop();
376 ExprIntrp_SyntaxError::Raise();
378 Handle(Expr_BinaryFunction) res =
379 new Expr_BinaryFunction(thefunc,arg1,op);
380 ExprIntrp_Recept.Push(res);
383 Expr_Array1OfGeneralExpression tabarg(1,nbargs);
384 Handle(Expr_GeneralExpression) arg;
386 for (Standard_Integer i = 1; i< nbargs; i++) {
387 arg = ExprIntrp_Recept.Pop();
389 ExprIntrp_SyntaxError::Raise();
391 tabarg(nbargs-i) = arg;
393 Handle(Expr_PolyFunction) res =
394 new Expr_PolyFunction(thefunc,tabarg);
395 ExprIntrp_Recept.Push(res);
400 extern "C" void ExprIntrp_NextFuncArg()
405 extern "C" void ExprIntrp_EndFuncArg()
410 extern "C" void ExprIntrp_SumOperator()
412 Handle(Expr_GeneralExpression) op2 = ExprIntrp_Recept.Pop();
413 Handle(Expr_GeneralExpression) op1 = ExprIntrp_Recept.Pop();
414 Handle(Expr_Sum) sres = op1 + op2;
415 Handle(Expr_GeneralExpression) res = sres->ShallowSimplified();
416 ExprIntrp_Recept.Push(res);
419 extern "C" void ExprIntrp_MinusOperator()
421 Handle(Expr_GeneralExpression) op2 = ExprIntrp_Recept.Pop();
422 Handle(Expr_GeneralExpression) op1 = ExprIntrp_Recept.Pop();
423 Handle(Expr_Difference) res = op1 - op2;
424 ExprIntrp_Recept.Push(res->ShallowSimplified());
427 extern "C" void ExprIntrp_DivideOperator()
429 Handle(Expr_GeneralExpression) op2 = ExprIntrp_Recept.Pop();
430 Handle(Expr_GeneralExpression) op1 = ExprIntrp_Recept.Pop();
431 Handle(Expr_Division) res = op1 / op2;
432 ExprIntrp_Recept.Push(res->ShallowSimplified());
435 extern "C" void ExprIntrp_ExpOperator()
437 Handle(Expr_GeneralExpression) op2 = ExprIntrp_Recept.Pop();
438 Handle(Expr_GeneralExpression) op1 = ExprIntrp_Recept.Pop();
439 Handle(Expr_Exponentiate) res = new Expr_Exponentiate(op1,op2);
440 ExprIntrp_Recept.Push(res->ShallowSimplified());
443 extern "C" void ExprIntrp_ProductOperator()
445 Handle(Expr_GeneralExpression) op2 = ExprIntrp_Recept.Pop();
446 Handle(Expr_GeneralExpression) op1 = ExprIntrp_Recept.Pop();
447 Handle(Expr_Product) res = op1 * op2;
448 ExprIntrp_Recept.Push(res->ShallowSimplified());
451 extern "C" void ExprIntrp_UnaryMinusOperator()
453 Handle(Expr_GeneralExpression) op = ExprIntrp_Recept.Pop();
454 Handle(Expr_UnaryMinus) res = new Expr_UnaryMinus(op);
455 ExprIntrp_Recept.Push(res->ShallowSimplified());
458 extern "C" void ExprIntrp_VariableIdentifier()
461 ExprIntrp_GetResult(name);
462 TCollection_AsciiString thename(name);
463 Handle(Expr_NamedExpression) nameexp = ExprIntrp_Recept.GetNamed(thename);
464 if (nameexp.IsNull()) {
465 nameexp = new Expr_NamedUnknown(thename);
466 ExprIntrp_Recept.Use(nameexp);
468 ExprIntrp_Recept.Push(nameexp);
471 extern "C" void ExprIntrp_NumValue()
475 nbcar = ExprIntrp_GetResult(num);
476 Standard_Real value = Atof(num);
477 Handle(Expr_NumericValue) nval = new Expr_NumericValue(value);
478 ExprIntrp_Recept.Push(nval);
481 extern "C" void ExprIntrp_AssignVariable()
483 ExprIntrp_GetResult(ExprIntrp_assname);
486 extern "C" void ExprIntrp_Deassign()
489 ExprIntrp_GetResult(name);
490 TCollection_AsciiString thename(name);
491 Handle(Expr_NamedExpression) nameexp = ExprIntrp_Recept.GetNamed(thename);
492 if (nameexp.IsNull()) {
493 ExprIntrp_SyntaxError::Raise();
495 if (!nameexp->IsKind(STANDARD_TYPE(Expr_NamedUnknown))) {
496 ExprIntrp_SyntaxError::Raise();
498 Handle(Expr_NamedUnknown) var = Handle(Expr_NamedUnknown)::DownCast(nameexp);
502 extern "C" void ExprIntrp_DefineFunction()
504 ExprIntrp_funcdefname = ExprIntrp_Recept.PopName();
505 ExprIntrp_Recept.PushValue(ExprIntrp_nbargs);
508 extern "C" void ExprIntrp_close()
510 ExprIntrp_stop_string();
513 extern "C" void ExprIntrperror(char* msg)
516 ExprIntrp_SyntaxError::Raise(msg);
520 extern "C" void ExprIntrp_EndOfEqual()
522 Handle(Expr_GeneralExpression) memb2 = ExprIntrp_Recept.Pop();
523 Handle(Expr_GeneralExpression) memb1 = ExprIntrp_Recept.Pop();
524 Handle(Expr_Equal) res = new Expr_Equal(memb1,memb2);
525 ExprIntrp_Recept.PushRelation(res);
528 extern "C" void ExprIntrp_EndOfRelation()
530 Handle(Expr_SystemRelation) sys;
531 Handle(Expr_GeneralRelation) currel;
532 Handle(Expr_GeneralRelation) oldrel;
533 while (!ExprIntrp_Recept.IsRelStackEmpty()) {
534 currel = ExprIntrp_Recept.PopRelation();
538 else if (!oldrel.IsNull()) {
539 sys = new Expr_SystemRelation(oldrel);
547 ExprIntrp_Recept.PushRelation(currel);
550 ExprIntrp_Recept.PushRelation(sys);
554 extern "C" void ExprIntrp_EndOfAssign()
556 Handle(Expr_NamedExpression) namexp = ExprIntrp_Recept.GetNamed(ExprIntrp_assname);
557 Handle(Expr_NamedUnknown) namu;
558 if (namexp.IsNull()) {
559 namu = new Expr_NamedUnknown(ExprIntrp_assname);
560 ExprIntrp_Recept.Use(namu);
563 if (!namexp->IsKind(STANDARD_TYPE(Expr_NamedUnknown))) {
564 ExprIntrp_SyntaxError::Raise();
566 namu = Handle(Expr_NamedUnknown)::DownCast(namexp);
568 namu->Assign(ExprIntrp_Recept.Pop());
571 extern "C" void ExprIntrp_EndOfFuncDef()
573 Handle(Expr_GeneralExpression) theexp = ExprIntrp_Recept.Pop();
574 Standard_Integer nbargs = ExprIntrp_Recept.PopValue();
575 Expr_Array1OfNamedUnknown vars(1,nbargs);
576 Expr_Array1OfNamedUnknown internvars(1,nbargs);
578 for (i=nbargs; i > 0; i--) {
579 vars(i) = Handle(Expr_NamedUnknown)::DownCast(ExprIntrp_Recept.Pop());
580 internvars(i) = Handle(Expr_NamedUnknown)::DownCast(vars(i)->Copy());
582 theexp = Expr::CopyShare(theexp); // ATTENTION, PROTECTION BUG STACK
583 for (i=1; i<= nbargs; i++) {
584 if (theexp->Contains(vars(i))) {
585 theexp->Replace(vars(i),internvars(i));
588 if (theexp == vars(i)) {
589 theexp = internvars(i);
593 Handle(Expr_NamedFunction) thefunc =
594 new Expr_NamedFunction(ExprIntrp_funcdefname,
597 ExprIntrp_Recept.Use(thefunc);
600 extern "C" void ExprIntrp_ConstantIdentifier()
603 ExprIntrp_GetResult(name);
604 TCollection_AsciiString thename(name);
605 ExprIntrp_Recept.PushName(thename);
608 extern "C" void ExprIntrp_ConstantDefinition()
610 TCollection_AsciiString name = ExprIntrp_Recept.PopName();
613 nbcar = ExprIntrp_GetResult(num);
614 Standard_Real val = Atof(num);
615 Handle(Expr_NamedConstant) theconst = new Expr_NamedConstant(name,val);
616 ExprIntrp_Recept.Use(theconst);
617 ExprIntrp_Recept.Push(theconst);
621 extern "C" void ExprIntrp_Sumator()
623 Handle(Expr_NumericValue) number = Handle(Expr_NumericValue)::DownCast(ExprIntrp_Recept.Pop());
624 Standard_Integer nb = (Standard_Integer) number->GetValue();
625 Handle(Expr_GeneralExpression) inc = ExprIntrp_Recept.Pop();
626 Handle(Expr_GeneralExpression) first = ExprIntrp_Recept.Pop();
627 Handle(Expr_NamedUnknown) var = Handle(Expr_NamedUnknown)::DownCast(ExprIntrp_Recept.Pop());
628 Handle(Expr_GeneralExpression) theexp = ExprIntrp_Recept.Pop();
629 Standard_Boolean thesame = (theexp == var);
630 Handle(Expr_GeneralExpression) cur = Expr::CopyShare(first);
631 Handle(Expr_GeneralExpression) res;
632 Handle(Expr_GeneralExpression) member;
633 Expr_SequenceOfGeneralExpression seq;
634 for (Standard_Integer i=1; i<= nb; i++) {
639 member = Expr::CopyShare(theexp);
640 member->Replace(var,cur);
643 cur = (cur + inc)->ShallowSimplified();
645 res = new Expr_Sum(seq);
646 ExprIntrp_Recept.Push(res->ShallowSimplified());
649 extern "C" void ExprIntrp_Productor()
651 Handle(Expr_NumericValue) number = Handle(Expr_NumericValue)::DownCast(ExprIntrp_Recept.Pop());
652 Standard_Integer nb = (Standard_Integer) number->GetValue();
653 Handle(Expr_GeneralExpression) inc = ExprIntrp_Recept.Pop();
654 Handle(Expr_GeneralExpression) first = ExprIntrp_Recept.Pop();
655 Handle(Expr_NamedUnknown) var = Handle(Expr_NamedUnknown)::DownCast(ExprIntrp_Recept.Pop());
656 Handle(Expr_GeneralExpression) theexp = ExprIntrp_Recept.Pop();
657 Standard_Boolean thesame = (theexp == var);
658 Handle(Expr_GeneralExpression) cur = Expr::CopyShare(first);
659 Handle(Expr_GeneralExpression) res;
660 Handle(Expr_GeneralExpression) member;
661 Expr_SequenceOfGeneralExpression seq;
662 for (Standard_Integer i=1; i<= nb; i++) {
667 member = Expr::CopyShare(theexp);
668 member->Replace(var,cur);
671 cur = (cur + inc)->ShallowSimplified();
673 res = new Expr_Product(seq);
674 ExprIntrp_Recept.Push(res->ShallowSimplified());