44 namespace MultiRegions
47 ExpList3DHomogeneous2D::ExpList3DHomogeneous2D():
59 const bool dealiasing):
72 const bool dealiasing,
89 for(j = 0; j < nel; ++j)
97 for(n = 1; n < (ny*nz); ++n)
100 for(j = 0; j < nel; ++j)
102 (*m_exp).push_back((*
m_exp)[j]);
123 if(DeclareLinesSetCoeffPhys)
128 for(
int n = 0; n <
m_lines.num_elements(); ++n)
141 const std::vector<unsigned int> &eIDs,
142 const bool DeclareLinesSetCoeffPhys):
147 if(DeclareLinesSetCoeffPhys)
150 std::vector<unsigned int> eIDsLine;
151 int nel = eIDs.size()/
m_lines.num_elements();
153 for (
int i = 0; i < nel; ++i)
155 eIDsLine.push_back(eIDs[i]);
164 for(
int n = 0; n <
m_lines.num_elements(); ++n)
183 int ncoeffs_per_line =
m_lines[0]->GetNcoeffs();
184 int npoints_per_line =
m_lines[0]->GetTotPoints();
186 int nyzlines =
m_lines.num_elements();
195 int nel =
m_lines[0]->GetExpSize();
201 for(cnt = n = 0; n < nyzlines; ++n)
204 m_lines[n]->SetPhysArray(tmparray =
m_phys + npoints_per_line*n);
206 for(i = 0; i < nel; ++i)
209 m_phys_offset[cnt] =
m_lines[n]->GetPhys_Offset(i) + n*npoints_per_line;
210 m_offset_elmt_id[cnt++] =
m_lines[n]->GetOffset_Elmt_Id(i) + n*nel;
241 (*m_exp)[eid]->GetCoords(x);
244 for(m = 0; m < nzlines; ++m)
246 for(j = 0; j < nylines; ++j)
248 for(n = 0; n < npoints; ++n)
250 Vmath::Fill(1,x[n],tmp_xc = xc0 + n +(j*npoints) + (m*npoints*nylines), 1);
251 Vmath::Fill(1,y[j],tmp_xc = xc1 + n +(j*npoints) + (m*npoints*nylines), 1);
252 Vmath::Fill(1,z[m],tmp_xc = xc2 + n +(j*npoints) + (m*npoints*nylines), 1);
280 int npoints =
m_lines[0]->GetTotPoints();
301 for(m = 0; m < nzlines; ++m)
303 for(j = 0; j < nylines; ++j)
305 for(n = 0; n < npoints; ++n)
307 Vmath::Fill(1,x[n],tmp_xc = xc0 + n +(j*npoints) + (m*npoints*nylines), 1);
308 Vmath::Fill(1,y[j],tmp_xc = xc1 + n +(j*npoints) + (m*npoints*nylines), 1);
309 Vmath::Fill(1,z[m],tmp_xc = xc2 + n +(j*npoints) + (m*npoints*nylines), 1);
325 int nquad0 = (*m_exp)[expansion]->GetNumPoints(0);
332 coords[1] = coords[0] + nquad0*nquad1*nquad2;
333 coords[2] = coords[1] + nquad0*nquad1*nquad2;
335 GetCoords(expansion,coords[0],coords[1],coords[2]);
337 outfile <<
"Zone, I=" << nquad0 <<
", J=" << nquad1 <<
",K="
338 << nquad2 <<
", F=Block" << std::endl;
340 for(j = 0; j < 3; ++j)
342 for(i = 0; i < nquad0*nquad1*nquad2; ++i)
344 outfile << coords[j][i] <<
" ";
346 outfile << std::endl;
354 int nquad0 = (*m_exp)[expansion]->GetNumPoints(0);
357 int ntot = nquad0*nquad1*nquad2;
358 int ntotminus = (nquad0-1)*(nquad1-1)*(nquad2-1);
364 GetCoords(expansion,coords[0],coords[1],coords[2]);
366 outfile <<
" <Piece NumberOfPoints=\""
367 << ntot <<
"\" NumberOfCells=\""
368 << ntotminus <<
"\">" << endl;
369 outfile <<
" <Points>" << endl;
370 outfile <<
" <DataArray type=\"Float32\" "
371 <<
"NumberOfComponents=\"3\" format=\"ascii\">" << endl;
373 for (i = 0; i < ntot; ++i)
375 for (j = 0; j < 3; ++j)
377 outfile << coords[j][i] <<
" ";
382 outfile <<
" </DataArray>" << endl;
383 outfile <<
" </Points>" << endl;
384 outfile <<
" <Cells>" << endl;
385 outfile <<
" <DataArray type=\"Int32\" "
386 <<
"Name=\"connectivity\" format=\"ascii\">" << endl;
387 for (i = 0; i < nquad0-1; ++i)
389 for (j = 0; j < nquad1-1; ++j)
391 for (k = 0; k < nquad2-1; ++k)
393 outfile << k*nquad0*nquad1 + j*nquad0 + i <<
" "
394 << k*nquad0*nquad1 + j*nquad0 + i + 1 <<
" "
395 << k*nquad0*nquad1 + (j+1)*nquad0 + i + 1 <<
" "
396 << k*nquad0*nquad1 + (j+1)*nquad0 + i <<
" "
397 << (k+1)*nquad0*nquad1 + j*nquad0 + i <<
" "
398 << (k+1)*nquad0*nquad1 + j*nquad0 + i + 1 <<
" "
399 << (k+1)*nquad0*nquad1 + (j+1)*nquad0 + i + 1 <<
" "
400 << (k+1)*nquad0*nquad1 + (j+1)*nquad0 + i <<
" " << endl;
405 outfile <<
" </DataArray>" << endl;
406 outfile <<
" <DataArray type=\"Int32\" "
407 <<
"Name=\"offsets\" format=\"ascii\">" << endl;
408 for (i = 0; i < ntotminus; ++i)
410 outfile << i*8+8 <<
" ";
413 outfile <<
" </DataArray>" << endl;
414 outfile <<
" <DataArray type=\"UInt8\" "
415 <<
"Name=\"types\" format=\"ascii\">" << endl;
416 for (i = 0; i < ntotminus; ++i)
421 outfile <<
" </DataArray>" << endl;
422 outfile <<
" </Cells>" << endl;
423 outfile <<
" <PointData>" << endl;
439 for(
int m = 0; m < nzlines; ++m)
441 for(
int n = 0; n < nylines; ++n)
443 errL2 =
m_lines[n+(m*nylines)]->
L2(inarray + cnt, soln + cnt);
Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expa...
virtual void v_WriteTecplotZone(std::ostream &outfile, int expansion)
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
NekOptimize::GlobalOptParamSharedPtr m_globalOptParam
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
const boost::shared_ptr< LocalRegions::ExpansionVector > GetExp() const
This function returns the vector of elements in the expansion.
Array< OneD, NekDouble > m_phys
The global expansion evaluated at the quadrature points.
NekDouble m_lhom_z
Width of homogeneous direction z.
Array< OneD, NekDouble > m_coeffs
Concatenation of all local expansion coefficients.
int GetExpSize(void)
This function returns the number of elements in the expansion.
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
Array< OneD, ExpListSharedPtr > m_lines
Vector of ExpList, will be filled with ExpList1D.
virtual void v_WriteVtkPieceHeader(std::ostream &outfile, int expansion, int istrip)
Array< OneD, int > m_coeff_offset
Offset of elemental data into the array m_coeffs.
LibUtilities::BasisSharedPtr m_homogeneousBasis_z
Base expansion in z direction.
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
boost::shared_ptr< LocalRegions::ExpansionVector > m_exp
The list of local expansions.
int GetTotPoints(void) const
Returns the total number of quadrature points m_npoints .
boost::shared_ptr< ExpList1D > ExpList1DSharedPtr
Shared pointer to an ExpList1D object.
virtual NekDouble v_L2(const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
void SetCoeffPhys(void)
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Array< OneD, int > m_phys_offset
Offset of elemental data into the array m_phys.
virtual void v_GetCoords(Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1, Array< OneD, NekDouble > &coord_2)
ExpList3DHomogeneous2D()
Default constructor.
int m_ncoeffs
The total number of local degrees of freedom. m_ncoeffs .
Array< OneD, int > m_offset_elmt_id
Array containing the element id m_offset_elmt_id[n] that the n^th consecutive block of data in m_coef...
NekDouble m_lhom_y
Width of homogeneous direction y.
virtual ~ExpList3DHomogeneous2D()
Destructor.
LibUtilities::SessionReaderSharedPtr m_session
Session.
void Sadd(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Add vector y = alpha + x.
LibUtilities::BasisSharedPtr m_homogeneousBasis_y
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
This class is the abstraction of a one-dimensional multi-elemental expansions which is merely a colle...
void GetCoords(Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coord_2=NullNekDouble1DArray)
This function calculates the coordinates of all the elemental quadrature points . ...
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
boost::shared_ptr< MeshGraph > MeshGraphSharedPtr
NekDouble L2(const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
This function calculates the error with respect to soln of the global spectral/hp element approximat...
Abstraction of a one-dimensional multi-elemental expansion which is merely a collection of local expa...
Describes the specification for a Basis.