Dear all,
I want to create a C0-continuous field from interpolated input data. The input 
field is obtained from a FVM code and was filtered with a gaussian in space and 
interpolated to the quadrature points. I store the field in a  ContField2D
m_bf = Array<OneD, MultiRegions::ExpListSharedPtr>(m_spacedim + 2);
m_bf[i] = 
MemoryManager<MultiRegions::ContField2D>::AllocateSharedPtr(m_session, 
m_graph);
copy the interpolated data into its m_phys array
Vmath::Vcopy(nq, filteredField[i], 1, m_bf[i]->UpdatePhys(), 1);
and then try to smooth the field. Until now, i just did:
m_bf[i]->IProductWRTBase(m_bf[i]->GetPhys(), m_bf[i]->UpdateCoeffs());
m_bf[i]->MultiplyByElmtInvMass(m_bf[i]->GetCoeffs(), m_bf[i]->UpdateCoeffs());
m_bf[i]->BwdTrans(m_bf[i]->GetCoeffs(), m_bf[i]->UpdatePhys());
which yields the field depicted in DG.png. However, this does not ensure C0-
continuity, so i tried to call ExpList::SmoothField afterward:
m_bf[i]->SmoothField(m_bf[i]->UpdatePhys());
However, this produces a distorted field as shown in smoothed.png.
I am using a 6th order modified expansion basis:
<E COMPOSITE="C[10]" NUMMODES="6" FIELDS="p,u,v,w" TYPE="MODIFIED"/>
and the mesh is depicted in smoothed_mesh.png.
What is the cause of this problem and what can i do do fix it?
Cheers,
Kilian
-- 
Kilian Lackhove, M.Sc.
Fachgebiet für Energie- und Kraftwerkstechnik
L1|08 114
Technische Universität Darmstadt
Jovanka-Bontschits-Straße 2
D-64287 Darmstadt
Germany
Tel.: +49 6151 16 - 28915
Fax: +49 6151 16 - 6555
e-mail: lackhove(a)ekt.tu-darmstadt.de