Ejemplo n.º 1
0
 void DecrSortEigenStuff(void)
 {
   Tridiagonal(); //diagonalize the matrix.
   QLAlgorithm(); //
   DecreasingSort();
   GuaranteeRotation();
 }
/********************************************************
* num_sim_body *
* *
* Purpose: Simulate numerical solution for the interior region
* *
* Parameters *
* i, and j. Index trackers
* *
* Returns *
* None (void) *
********************************************************/
PROGRAMDATA num_sim_body(PROGRAMDATA pd)
{
   int y;
   pd = triCopy(pd);
   Tridiagonal(pd.Ny, pd.tri_c, pd.tri_a, pd.tri_b, pd.tri_y);
   for(y=0;y<pd.NyInter;y++){
      pd.pp2[y].Temp = pd.tri_y[y];
   }
   pd.pp2[pd.NyInter].Temp=1;
	return pd;
}