From: david <david>
Date: Thu, 24 Jan 2002 15:51:26 +0000 (+0000)
Subject: Latest update from JSBSim.
X-Git-Url: https://git.mxchange.org/?a=commitdiff_plain;h=45807b31268b410c4cb6e4af20c5759b69e73cc9;p=flightgear.git

Latest update from JSBSim.
---

diff --git a/src/FDM/JSBSim/FGInitialCondition.cpp b/src/FDM/JSBSim/FGInitialCondition.cpp
index 29dc1f10f..a76ef2a92 100644
--- a/src/FDM/JSBSim/FGInitialCondition.cpp
+++ b/src/FDM/JSBSim/FGInitialCondition.cpp
@@ -501,6 +501,9 @@ bool FGInitialCondition::getMachFromVcas(double *Mach,double vcas) {
 bool FGInitialCondition::getAlpha(void) {
   bool result=false;
   double guess=theta-gamma;
+  
+  if(vt < 0.01) return 0;
+  
   xlo=xhi=0;
   xmin=fdmex->GetAircraft()->GetAlphaCLMin();
   xmax=fdmex->GetAircraft()->GetAlphaCLMax();
@@ -521,6 +524,9 @@ bool FGInitialCondition::getAlpha(void) {
 bool FGInitialCondition::getTheta(void) {
   bool result=false;
   double guess=alpha+gamma;
+  
+  if(vt < 0.01) return 0;
+  
   xlo=xhi=0;
   xmin=-89;xmax=89;
   sfunc=&FGInitialCondition::GammaEqOfTheta;
@@ -555,7 +561,6 @@ double FGInitialCondition::GammaEqOfTheta(double Theta) {
 double FGInitialCondition::GammaEqOfAlpha(double Alpha) {
   double a,b,c;
   double sAlpha,cAlpha;
-
   sAlpha=sin(Alpha); cAlpha=cos(Alpha);
   a=wdown + vt*cAlpha*cbeta + uw;
   b=vt*sphi*sbeta + vw*sphi;
@@ -666,6 +671,7 @@ bool FGInitialCondition::solve(double *y,double x)
   i=0;
   while ((fabs(d) > eps) && (i < 100)) {
     d=(x3-x1)/d0;
+    cout << "f3-f1= " << f3-f1 << endl;
     x2 = x1-d*d0*f1/(f3-f1);
     
     f2=(this->*sfunc)(x2)-x;