Extruding a model with friction = frictionweertman leads to an inconsistent model.
%...
md.friction = frictionweertman
md.friction.C = c * ones(md.mesh.numberofvertices, 1);
md.friction.m = m * ones(md.mesh.numberofelements, 1);
%...
md = md.extrude(md, 7, 1);
%...
md = solve(md, 'tr');
leads to consistency error, length of friction.C/m doesn't match number of vertices/elements.
I believe this patch fixes it:
Index: src/m/classes/frictionweertman.m
===================================================================
--- src/m/classes/frictionweertman.m (revision 27645)
+++ src/m/classes/frictionweertman.m (working copy)
@@ -23,8 +23,8 @@
end
end % }}}
function self = extrude(self,md) % {{{
- md.friction.C=project3d(md,'vector',md.friction.C,'type','node','layer',1);
- md.friction.m=project3d(md,'vector',md.friction.m,'type','element');
+ self.C=project3d(md,'vector',self.C,'type','node','layer',1);
+ self.m=project3d(md,'vector',self.m,'type','element');
end % }}}
function self = setdefaultparameters(self) % {{{
