#include <aerodata.h>

Inheritance diagram for AeroMemory:

Collaboration diagram for AeroMemory:

Public Member Functions
	AeroMemory (DriveCaller *pt)

virtual	~AeroMemory (void)

void	Predict (int i, doublereal alpha, doublereal &alf1, doublereal &alf2)

void	Update (int i)

void	SetNumPoints (int i)

int	GetNumPoints (void) const

Protected Member Functions
virtual int	StorageSize (void) const =0

Protected Attributes
DriveCaller *	pTime

Private Attributes
doublereal *	a

doublereal *	t

integer	iPoints

int	numUpdates

Detailed Description

Definition at line 59 of file aerodata.h.

Constructor & Destructor Documentation

AeroMemory::AeroMemory ( DriveCaller * pt )

Definition at line 44 of file aerodata.cc.

References NO_OP.

 : a(0), t(0), iPoints(0), numUpdates(0), pTime(pt)
 {
         NO_OP;
 }

AeroMemory::~AeroMemory ( void )

virtual

Definition at line 50 of file aerodata.cc.

References a, iPoints, pTime, SAFEDELETE, and SAFEDELETEARR.

 {
         if (iPoints > 0) {
                 if (pTime) {
                         SAFEDELETE(pTime);
                 }
 
                 if (a) {
                         SAFEDELETEARR(a);
                 }
         }
 }

Member Function Documentation

int AeroMemory::GetNumPoints ( void ) const

Definition at line 162 of file aerodata.cc.

References iPoints.

Referenced by TheodorsenAeroData::AssRes().

 {
         return iPoints;
 }

void AeroMemory::Predict	(	int	i,
		doublereal	alpha,
		doublereal &	alf1,
		doublereal &	alf2
	)

Definition at line 65 of file aerodata.cc.

References a, ASSERT, DriveCaller::dGet(), iPoints, numUpdates, order, polcoe(), pTime, StorageSize(), and t.

Referenced by STAHRAeroData::GetForces(), C81AeroData::GetForces(), C81MultipleAeroData::GetForces(), and C81InterpolatedAeroData::GetForces().

 {
         /* FIXME: this should be s, but I don't want a malloc here */
         doublereal coe[3];
         int s = StorageSize();
 #ifdef USE_POLCOE
         /*
          * FIXME: actually this is not the order of the polynomial,
          * but the number of coefficients, e.g. a second-order polynomial
          * has three coefficients :)
          */
         integer         order = s;
 #endif /* USE_POLCOE */
 
         ASSERT(s > 0);
         ASSERT(i < iPoints);
 
         doublereal *aa = a + s*i;
         doublereal *tt = t + s*i;
 
         aa[s-1] = alpha;
         tt[s-1] = pTime->dGet();
 
         if (numUpdates >= s) {
 #ifdef USE_POLCOE
                 __FC_DECL__(polcoe)(tt, aa, &order, coe);
 #else /* !USE_POLCOE */
                 coe[2] = ((aa[2]-aa[0])/(tt[2]-tt[0])
                                 - (aa[1]-aa[0])/(tt[1]-tt[0]))/(tt[2]-tt[1]);
                 coe[1] = (aa[1]-aa[0])/(tt[1]-tt[0]) - (tt[1]+tt[0])*coe[2];
 
 #if 0   /* not used ... */
                 coe[1] = aa[0] - tt[0]*coe[1] - tt[0]*tt[0]*coe[2];
 #endif
 #endif /* !USE_POLCOE */
 
 #if 0
                         std::cerr << "aa[0:2]= " << aa[0] << "," << aa[1] << "," << aa[2] << std::endl
                                 << "tt[0:2]= " << tt[0] << "," << tt[1] << "," << tt[2] << std::endl
                                 << "coe[0:2]=" << coe[0] << "," << coe[1] << "," << coe[2] << std::endl;
 #endif /* 0 */
 
                 alf1 = coe[1]+2.*coe[2]*tt[2];
                 alf2 = 2.*coe[2];
         } else {
                 alf1 = 0.;
                 alf2 = 0.;
         }
 }

Here is the call graph for this function:

void AeroMemory::SetNumPoints ( int i )

Definition at line 141 of file aerodata.cc.

References a, iPoints, SAFENEWARR, StorageSize(), and t.

Referenced by AeroData::AeroData().

 {
         int s = StorageSize();
 
         iPoints = i;
 
         if (s > 0) {
                 SAFENEWARR(a, doublereal, 2*s*iPoints);
                 t = a + s*iPoints;
 
 #ifdef HAVE_MEMSET
                 memset(a, 0, 2*s*iPoints*sizeof(doublereal));
 #else /* !HAVE_MEMSET */
                 for (int i = 0; i < 2*s*iPoints; i++) {
                         a[i] = 0.;
                 }
 #endif /* !HAVE_MEMSET */
         }
 }

Here is the call graph for this function:

virtual int AeroMemory::StorageSize ( void ) const

protectedpure virtual

Implemented in AeroData.

Referenced by Predict(), SetNumPoints(), and Update().

void AeroMemory::Update ( int i )

Definition at line 116 of file aerodata.cc.

References a, ASSERT, iPoints, numUpdates, StorageSize(), and t.

 {
         int s = StorageSize();
         if (s > 0) {
                 /*
                  * shift back angle of attack and time
                  * for future interpolation
                  */
                 doublereal *aa = a + s*i;
                 doublereal *tt = t + s*i;
 
                 ASSERT(i < iPoints);
 
                 for (int j = 1; j < s; j++) {
                         aa[j-1] = aa[j];
                         tt[j-1] = tt[j];
                 }
 
                 if (i == 0) {
                         numUpdates++;
                 }
         }
 }