Commit 1ae1adc6 authored by Gabriele Rocco's avatar Gabriele Rocco
Browse files

3DHomo1D Floquet analysis; removed printing lines


git-svn-id: https://gforge.sci.utah.edu/svn/nektar/trunk@3798 305cdda6-5ce1-45b3-a98d-dfc68c8b3305
parent c3b7dad5
......@@ -2329,7 +2329,6 @@ namespace Nektar
const LibUtilities::PointsKey pkey2(npoints[cnt+b],LibUtilities::eGaussLobattoLegendre);
pkey = pkey2;
}
cout << "basis " << basis[b] << endl;
LibUtilities::BasisKey bkey(basis[b],nmodes[cnt+b],pkey);
bkeyvec.push_back(bkey);
}
......@@ -2482,7 +2481,6 @@ namespace Nektar
std::vector< std::vector<LibUtilities::PointsType> > &pointstype)
{
cout << "Inside set Expansion" << endl;
int i,j,k,g,h,cnt,id;
GeometrySharedPtr geom;
......
......@@ -610,23 +610,15 @@ namespace Nektar
}
FILE *pFile3;
pFile3= fopen("Base_U.txt", "w");
for(int s=0; s< m_base[0]->GetNpoints(); ++s)
{
fprintf(pFile3, "base[%i]= %10.20lf\t \n",s,m_base[0]->GetPhys()[s]);
}
fclose(pFile3);
}
//@ToDo: Post Processing for ModifiedBase
//std::string outname ="BaseFlow.bse";
//WriteFldBase(outname);
if(m_session->DefinesParameter("N_slices"))
{
m_nConvectiveFields = m_base.num_elements()-1;
for(int i=0; i<m_nConvectiveFields;++i)
{
......@@ -931,7 +923,9 @@ namespace Nektar
for (int i = 2; i < m_slices; i += 2)
{
phase = (i>>1) * BetaT;
Vmath::Svtvp(npoints, cos(phase),&inarray[i*npoints],1,&outarray[0],1,&outarray[0],1);
//-cos(phase) changed
Vmath::Svtvp(npoints, -cos(phase),&inarray[i*npoints],1,&outarray[0],1,&outarray[0],1);
Vmath::Svtvp(npoints, -sin(phase), &inarray[(i+1)*npoints], 1, &outarray[0], 1,&outarray[0],1);
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment