Skip to content
Snippets Groups Projects
Commit e3e9b96d authored by Porrmann, Maik's avatar Porrmann, Maik
Browse files

printBasis colors, add Todo comment

parent 29af37c0
No related branches found
No related tags found
No related merge requests found
......@@ -197,6 +197,8 @@ namespace Dune
transformator_.apply(out, x, element_->geometry());
}
// TODO sfinae protected hessian method
//! \brief Polynomial order of the shape functions
auto order() const
{
......
......@@ -90,8 +90,10 @@ namespace Dune
<< labelPos[0] + 0.1 * dirMatrix[i][0] << ","
<< labelPos[1] + 0.1 * dirMatrix[i][1] << " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
}
......@@ -103,15 +105,19 @@ namespace Dune
<< pos[0] + 0.1 * dirMatrix[i][0] << "," << pos[1] + 0.1 * dirMatrix[i][1]
<< " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
plotfile << "set arrow from " << pos[0] + 0.1 * dirMatrix[i][0] << ","
<< pos[1] + 0.1 * dirMatrix[i][1] << " to " << pos[0] + 0.2 * dirMatrix[i][0]
<< "," << pos[1] + 0.2 * dirMatrix[i][1] << " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
}
......@@ -123,16 +129,20 @@ namespace Dune
<< pos[0] + 0.1 * dirMatrix[0][0] << "," << pos[1] + 0.1 * dirMatrix[0][1]
<< " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
plotfile << "set arrow from " << pos[0] + 0.1 * dirMatrix[0][0] << ","
<< pos[1] + 0.1 * dirMatrix[0][1] << " to "
<< pos[0] + 0.1 * dirMatrix[0][0] + 0.1 * dirMatrix[1][0] << ","
<< pos[1] + 0.1 * dirMatrix[0][1] + 0.1 * dirMatrix[1][1] << " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
}
......@@ -145,15 +155,19 @@ namespace Dune
<< pos[0] + 0.1 * dirMatrix[i][0] << "," << pos[1] + 0.1 * dirMatrix[i][1]
<< " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
plotfile << "set arrow from " << pos[0] + 0.1 * dirMatrix[i][0] << ","
<< pos[1] + 0.1 * dirMatrix[i][1] << " to " << pos[0] + 0.2 * dirMatrix[i][0]
<< "," << pos[1] + 0.2 * dirMatrix[i][1] << " lt rgb "
<< (lv.tree().finiteElement().localCoefficients().isDirichlet(localIndex)
? "\"green\""
: "\"red\"")
? "\"green\" lw 4"
: (lv.tree().finiteElement().localCoefficients().isClamped(localIndex))
? "\"blue\" lw 3"
: "\"red\" lw 2")
<< std::endl;
}
}
......@@ -275,7 +289,8 @@ namespace Dune
plotfile << "set object circle front at " << pos[0] << "," << pos[1]
<< " size first 0.08 fc rgb \"blue\" lw 4" << std::endl;
}
writeElementInfo(plotfile, localView, pos, key, i); // writes directions
if (entity.hasBoundaryIntersections())
writeElementInfo(plotfile, localView, pos, key, i); // writes directions
}
else if (key.codim() == dim - 1) // suppose this is normal derivative dof //TODO this is not
// valid for the vast majority of finite elements
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment