ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-24 21:10:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

simd-mapping of point elements

Browse Source
This commit is contained in:
Joachim Schoeberl 2021-02-05 11:42:45 +01:00
parent 221d3f5a9a
commit 22aee3b3a5
1 changed files with 15 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
libsrc/interface/nginterface_v2.cpp
View File

@ -861,7 +861,14 @@ namespace netgen
SIMD<double> * x, size_t sx, SIMD<double> * x, size_t sx,
SIMD<double> * dxdxi, size_t sdxdxi) const SIMD<double> * dxdxi, size_t sdxdxi) const
{ {
cout << "MultiElementtransformation<0,2> simd not implemented" << endl; //cout << "MultiElementtransformation<0,2> simd not implemented" << endl;
PointIndex pi = mesh->pointelements[elnr].pnum;
Point<3> xg = mesh->Point(pi);
if (x)
for (int j = 0; j < npts; j++)
for (int i = 0; i < 2; i++)
x[j*sx+i] = xg(i);
} }
template<> DLL_HEADER void Ngx_Mesh :: template<> DLL_HEADER void Ngx_Mesh ::
@ -870,7 +877,13 @@ namespace netgen
SIMD<double> * x, size_t sx, SIMD<double> * x, size_t sx,
SIMD<double> * dxdxi, size_t sdxdxi) const SIMD<double> * dxdxi, size_t sdxdxi) const
{ {
cout << "multi-eltrafo simd called, 0,1,simd" << endl; //cout << "multi-eltrafo simd called, 0,1,simd" << endl;
PointIndex pi = mesh->pointelements[elnr].pnum;
Point<3> xg = mesh->Point(pi);
if (x)
for (int j = 0; j < npts; j++)
for (int i = 0; i < 1; i++)
x[j*sx+i] = xg(i);
} }
template<> DLL_HEADER void Ngx_Mesh :: template<> DLL_HEADER void Ngx_Mesh ::