On 02/06/2014 04:13 PM, John Peterson wrote:
On Thu, Feb 6, 2014 at 2:07 PM, David Knezevic
<dkneze...@seas.harvard.edu <mailto:dkneze...@seas.harvard.edu>> wrote:
On 02/06/2014 03:27 PM, John Peterson wrote:
On Thu, Feb 6, 2014 at 11:28 AM, David Knezevic
<dkneze...@seas.harvard.edu <mailto:dkneze...@seas.harvard.edu>>
wrote:
Ah, good catch! If I change use_coupled_neighbor_dofs to do
"return libMesh::on_command_line
("--implicit_neighbor_dofs")", then it works fine, e.g. the
bandwidth in my test code with FIRST MONOMIAL becomes 4.
Excellent.
So I'd like to be able to enforce
--implicit_neighbor_dofs=false, even if all variables are
DISCONTINUOUS.
Making it take a following argument instead of just being a
boolean flag sounds good...
Any preferences on how to change the API to allow that?
Perhaps a second boolean argument to
use_coupled_neighbor_dofs called "use_commandline_value"?
The following is untested, but seems to compile. The API stays
the same but we additionally look for a following argument to
--implicit_neighbor_dofs which, if false, immediately returns
false from DofMap::use_coupled_neighbor_dofs without checking if
all the variables are discontinuous.
diff --git a/src/base/dof_map.C b/src/base/dof_map.C
index 7c50688..74b1a42 100644
--- a/src/base/dof_map.C
+++ b/src/base/dof_map.C
@@ -1471,6 +1471,24 @@ bool
DofMap::use_coupled_neighbor_dofs(const MeshBase& mesh) const
bool implicit_neighbor_dofs =
libMesh::on_command_line ("--implicit_neighbor_dofs");
+ // If the user specifies --implicit_neighbor_dofs=false, then
+ // presumably he knows what he is doing and we won't try to
+ // automatically turn it on even when all the variables are
+ // discontinuous.
+ if (implicit_neighbor_dofs)
+ {
+ // No flag provided defaults to 'true'
+ bool flag = true;
+ flag = libMesh::command_line_value
("--implicit_neighbor_dofs", flag);
+
+ if (!flag)
+ {
+ // The user said --implicit_neighbor_dofs 0, so he
knows
+ // what he is doing and really doesn't want it.
+ return false;
+ }
+ }
+
This looks good to me. Do you want to push it to the repo?
Sure, I can push it once I've had a chance to test it a bit...
OK, sounds good.
------------------------------------------------------------------------------
Managing the Performance of Cloud-Based Applications
Take advantage of what the Cloud has to offer - Avoid Common Pitfalls.
Read the Whitepaper.
http://pubads.g.doubleclick.net/gampad/clk?id=121051231&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
