MBDyn-1.7.3
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups
ginacdrive.h
Go to the documentation of this file.
1 /* $Header: /var/cvs/mbdyn/mbdyn/mbdyn-1.0/mbdyn/base/ginacdrive.h,v 1.7 2017/01/12 14:46:09 masarati Exp $ */
2 /*
3  * MBDyn (C) is a multibody analysis code.
4  * http://www.mbdyn.org
5  *
6  * Copyright (C) 1996-2017
7  *
8  * Pierangelo Masarati <masarati@aero.polimi.it>
9  * Paolo Mantegazza <mantegazza@aero.polimi.it>
10  *
11  * Dipartimento di Ingegneria Aerospaziale - Politecnico di Milano
12  * via La Masa, 34 - 20156 Milano, Italy
13  * http://www.aero.polimi.it
14  *
15  * Changing this copyright notice is forbidden.
16  *
17  * This program is free software; you can redistribute it and/or modify
18  * it under the terms of the GNU General Public License as published by
19  * the Free Software Foundation (version 2 of the License).
20  *
21  *
22  * This program is distributed in the hope that it will be useful,
23  * but WITHOUT ANY WARRANTY; without even the implied warranty of
24  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
25  * GNU General Public License for more details.
26  *
27  * You should have received a copy of the GNU General Public License
28  * along with this program; if not, write to the Free Software
29  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
30  */
31 
32 /* Classe di drivers che usano gradi di liberta' nodali
33  *
34  * Il valore del drive viene fatto dipendere da un grado di liberta' nodale
35  * in modo del tutto generale e trasparente. Ovvero in ingresso dati si
36  * associa il drive al grado di liberta' j-esimo del nodo i-esimo, con
37  * verifiche di consistenza. In esecuzione, la funzione propria dGet() del
38  * drive restituisce il valore del grado di liberta' associato
39  * (o della derivata, se il grado di liberta' e' differenziale).
40  */
41 
42 #ifndef GINACDRIVE_H
43 #define GINACDRIVE_H
44 
45 #include "drive.h"
46 #include <ginac/ginac.h>
47 
49 {
50 private:
51  // parameter symbols
52  GiNaC::symbol * gVar;
53 
54  // expression
55  GiNaC::ex gExpr;
56 
57  // derivative
58  GiNaC::ex gExprDVar;
59 
60 public:
61  GiNaCDriveCaller(const DriveHandler* pDH,
62  const std::string& var, const std::string& expression);
63  virtual ~GiNaCDriveCaller(void);
64 
65  /* Copia */
66  virtual DriveCaller* pCopy(void) const;
67 
68  virtual std::ostream& Restart(std::ostream& out) const;
69 
70  inline doublereal dGet(const doublereal& dVar) const;
71 
72  /* this is about drives that are differentiable */
73  virtual bool bIsDifferentiable(void) const;
74  virtual doublereal dGetP(const doublereal& dVar) const;
75 };
76 
77 inline doublereal
79 {
80  GiNaC::lst l;
81 
82  l.append(*gVar == dVar);
83 
84  GiNaC::ex f_expr = gExpr.subs(l);
85 
86  return GiNaC::ex_to<GiNaC::numeric>(f_expr).to_double();
87 }
88 
89 struct GiNaCDCR : public DriveCallerRead {
90  DriveCaller *
91  Read(const DataManager* pDM, MBDynParser& HP, bool bDeferred);
92 };
93 
94 #endif // GINACDRIVE_H
95 
DriveCaller * Read(const DataManager *pDM, MBDynParser &HP, bool bDeferred)
Definition: ginacdrive.cc:121
virtual ~GiNaCDriveCaller(void)
Definition: ginacdrive.cc:71
GiNaCDriveCaller(const DriveHandler *pDH, const std::string &var, const std::string &expression)
Definition: ginacdrive.cc:39
GiNaC::ex gExpr
Definition: ginacdrive.h:55
virtual std::ostream & Restart(std::ostream &out) const
Definition: ginacdrive.cc:95
virtual DriveCaller * pCopy(void) const
Definition: ginacdrive.cc:78
GiNaC::symbol * gVar
Definition: ginacdrive.h:52
GiNaC::ex gExprDVar
Definition: ginacdrive.h:58
virtual doublereal dGet(void) const
Definition: drive.h:489
virtual doublereal dGetP(void) const
Definition: drive.h:501
virtual bool bIsDifferentiable(void) const
Definition: ginacdrive.cc:103
double doublereal
Definition: colamd.c:52