Commit 92300889 authored by Robin Basso's avatar Robin Basso

new

parent ba0c50cf
......@@ -172,7 +172,7 @@ void DriverModifiedArnoldi::v_Execute(ostream &out)
m_comm->AllReduce(alpha[0], Nektar::LibUtilities::ReduceSum);
alpha[0] = std::sqrt(alpha[0]);
Vmath::Smul(ntot, 1.0/alpha[0], Kseq[0], 1, Kseq[0], 1);
cout<<"alpha[0] = "<<alpha[0];
// Scale angle if needed
if(m_BlowingSuction)
{
......@@ -195,7 +195,7 @@ cout<<"alpha[0] = "<<alpha[0];
m_comm->AllReduce(alpha[i], Nektar::LibUtilities::ReduceSum);
alpha[i] = std::sqrt(alpha[i]);
Vmath::Smul(ntot, 1.0/alpha[i], Kseq[i], 1, Kseq[i], 1);
cout<<"alpha[i] = "<<alpha[i];
// Scale angle if needed
if(m_BlowingSuction)
{
......
......@@ -909,7 +909,7 @@ namespace Nektar
m_bsbcParams->m_angleVel[0] = (1/m_bsbcParams->m_alpha) * sigma * exp(sigma * time);
m_bsbcParams->m_angle = (1/m_bsbcParams->m_alpha) * exp(sigma * time);
////////////////////////////////////////////////////////
cout<<"alpha = "<<m_bsbcParams->m_alpha<<endl;
if( m_bsbcParams->m_doOutput )
{
//m_bsbcParams->m_filterForces->Update(m_fields, time);
......
Markdown is supported
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