#include "trussheader.h"

void summary(double x,double y,double f) //simple summary of a vertex
{
      cout<<"\nX= "<<x<<" Y= "<<y<<" Force= "<<f<<endl;
 //     cout<<endl;
}

void difference(double x1,double y1,double x2,double y2,int W,int R)//difference between 2 points
{
      double changex=x2-x1,changey=y2-y1;
      cout<<"Difference in X= "<<changex<<endl;
      cout<<"Difference in Y= "<<changey<<endl;
      cout<<"A vecter from vertix "<<W<<" to vertix "<<R<<" is "<<x2-x1<<"i + "<<y2-y1<<"j\n";
      cout<<"Or a magnitude of "<<sqrt (pow (double (changex), double (2))+pow(double (changey),double (2)))<<" "; //need math.h
      cout<<"at the angle of= "<<(180/3.14159265359)* atan (changey/changex)<<" degrees\n";//must convert to degrees
}

void briefsum(int t)  //simple summary for total# members
{
     cout<<"You have "<<t<<" vertices, thus you have "<<(2*t-3)<<" member sections\n";
}

void test(int x)    //outputs error if user inputs a 1 or less vertex truss
{
if (x<=1)
{cout<<"Error: Must have 2 or more vertices to construct\n";
cout<<" proper truss problem. End Program.\n";
system("PAUSE");
exit(1);}
}

double angle(double x1,double y1,double x2,double y2)//difference between 2 points
{
      double angles;
      double changex=x2-x1;double changey=y2-y1;
      angles=atan(changey/changex);//in radians
      angles=fabs(angles);
      return angles;
}

double unitvector(double x1,double y1,double x2, double y2)//computes unit vector--unused
{
    double unitv;
    unitv=1/sqrt (pow (double (x2-x1), double (2))+pow(double (y2-y1),double (2)));
    return unitv;
}

void purevector(double x1,double y1,double x2,double y2,int W,int R)  //summary of a vector from 1 to 2
{
     cout<<"A vecter from vertix "<<W<<" to vertix "<<R<<" is "<<x2-x1<<"i + "<<y2-y1<<"j\n";
}

double vvector(double y1,double y2,int O,double y3,int T) //chooses the best vertex to use when beginning the calculations of load in that member
{
     double One, Two;
     One=y2-y1;
     Two=y3-y1;
     if (One==0)
        return T;
     else if (Two==0)
        return O;
}

void warning()       //simple note
{        cout<<"*************************************************\n";
         cout<<"*Note: Tension and upwards are taken as positive*\n";
         cout<<"*    : compression and downwards are negative   *\n";
         cout<<"*************************************************\n\n";

}

double truevalue(double Y1, double Y2,double F1,double M1)
{
 int X;
 if (Y2>Y1)
 X=5;
 else
 X=-5;
 if ((X>0) && (F1>0))
   return -M1;
 else if ((X>0) && (F1<0))
   return M1;
 else if ((X<0) && (F1>0))
   return M1;
 else if ((X<0) && (F1<0))
   return -M1;
 else
   exit(1);
}

double tensiontest2(double R1, double Xnew,double Xold,double F1,double M1) //signs correct sign to the load in a member
{
 int X;
 int Y;
 if (Xnew>R1)
 X=5;
 else
 X=-5;
 if (Xold>R1)
 Y=5;
 else
 Y=-5;
 if ((X>0) && (Y>0) && (F1>0)) //+++
   return -M1;
 else if ((X>0) && (Y>0) && (F1<0)) //++-
   return M1;
 else if ((X>0) &&(Y<0)&& (F1<0)) //+--
   return -M1;
 else if ((X<0) && (Y>0) && (F1>0)) //-++
   return M1;
 else if ((X<0) &&(Y>0)&& (F1<0)) //-+-
   return -M1;
 else if ((X<0)&&(Y<0)&&(F1<0))  //---
   return M1;
 else if ((X<0)&&(Y<0)&&(F1>0))   //--+
   return -M1;
 else if ((X>0)&&(Y<0)&&(F1>0))  //+-+
    return M1;
 else
    exit(1);

}

void value1(double X)       //based on results from tensiontest2, this will state if its in compression or tension
{
if (X>0)
cout<<"tension.\n";
else if (X<0)
cout<<"compression.\n";
else
cout<<"error!-value1 function\n";


}

double bestvector(double y1,double y2,int O,double y3,int T)   //chooses best vector after initial vector is chosen
{
     double One, Two;
     One=y2-y1;
     Two=y3-y1;
     if (One==0)
        return O;
     else if (Two==0)
        return T;
}