ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-13 17:48:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

little polish

Browse Source
This commit is contained in:
Joachim Schöberl 2019-08-21 18:38:13 +02:00
parent ca57759604
commit 5ada92bcdb
2 changed files with 7 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
libsrc/occ/occmeshsurf.cpp
View File

@ -210,13 +210,9 @@ namespace netgen
A_ = A_ * R_; A_ = A_ * R_;
Ainv_ = Trans(R_) * Ainv_; Ainv_ = Trans(R_) * Ainv_;
for (int i = 0; i < 2; i++) Amat = A_;
for (int j = 0; j < 2; j++) Amatinv = Ainv_;
{
Amat(i,j) = A_(i,j);
Amatinv(i,j) = Ainv_(i,j);
}
// temp = Amatinv * (psp2-psp1); // temp = Amatinv * (psp2-psp1);

7
libsrc/visualization/vssolution.cpp
View File

@ -2725,10 +2725,11 @@ namespace netgen
{ {
// NgProfiler::RegionTimer reg2 (timer2); // NgProfiler::RegionTimer reg2 (timer2);
int nse = mesh->GetNSE(); // int nse = mesh->GetNSE();
for (int i = 0; i < nse; i++) // for (int i = 0; i < nse; i++)
for (SurfaceElementIndex i : mesh->SurfaceElements().Range())
{ {
ELEMENT_TYPE type = mesh->SurfaceElement(i+1).GetType(); ELEMENT_TYPE type = (*mesh)[i].GetType();
double val; double val;
bool considerElem = (type == QUAD) bool considerElem = (type == QUAD)
? GetSurfValue (sol, i, -1, 0.5, 0.5, comp, val) ? GetSurfValue (sol, i, -1, 0.5, 0.5, comp, val)