@@ -175,18 +175,17 @@ void HangingNode ::postInitialize()
175175 }
176176 }
177177
178- // if (this->hasDofID(R_u) || this->hasDofID(R_v) || this->hasDofID(R_w)) {
179178 Dof * dof ;
180179 const int nnodes = e -> giveNumberOfNodes ();
181180
181+
182+ //Deal with rotations.
182183 // Compute shape function gradients once
183184 FloatMatrix dNdX ;
184185 fei -> evaldNdx (dNdX , lcoords , FEIElementGeometryWrapper (e ));
185186
186187 // === THETA_X = 0.5 * (du_z/dy - du_y/dz) ===
187188 if (this -> hasDofID (R_u )) {
188- printf ("R_u is still active on node %d!\n" , this -> giveNumber ());
189-
190189 FloatArray coeffs ;
191190 IntArray masterNodeIDs , masterDofIDs ;
192191
@@ -220,22 +219,13 @@ void HangingNode ::postInitialize()
220219 }
221220
222221
223-
224222 if (SlaveDof * sdof = dynamic_cast < SlaveDof * > (this -> giveDofWithID (R_u ))) {
225-
226- printf ("Slave DOF: node %d, dofID %d\n" , this -> giveNumber (), sdof -> giveDofID ());
227- for (int i = 1 ; i <= coeffs .giveSize (); ++ i ) {
228- printf (" Master %d: node %d, dofID %d, coeff = %.4e\n" ,
229- i , masterNodeIDs .at (i ), masterDofIDs .at (i ), coeffs .at (i ));
230- }
231-
232223 sdof -> initialize (masterNodeIDs , masterDofIDs , coeffs );
233224 }
234225 }
235226
236227 // === THETA_Y = 0.5 * (du_x/dz - du_z/dx) ===
237228 if (this -> hasDofID (R_v )) {
238- printf ("R_v is still active on node %d!\n" , this -> giveNumber ());
239229 FloatArray coeffs ;
240230 IntArray masterNodeIDs , masterDofIDs ;
241231
@@ -269,20 +259,12 @@ void HangingNode ::postInitialize()
269259 }
270260
271261 if (SlaveDof * sdof = dynamic_cast < SlaveDof * > (this -> giveDofWithID (R_v ))) {
272-
273- printf ("Slave DOF: node %d, dofID %d\n" , this -> giveNumber (), sdof -> giveDofID ());
274- for (int i = 1 ; i <= coeffs .giveSize (); ++ i ) {
275- printf (" Master %d: node %d, dofID %d, coeff = %.4e\n" ,
276- i , masterNodeIDs .at (i ), masterDofIDs .at (i ), coeffs .at (i ));
277- }
278-
279262 sdof -> initialize (masterNodeIDs , masterDofIDs , coeffs );
280263 }
281264 }
282265
283266 // === THETA_Z = 0.5 * (du_y/dx - du_x/dy) ===
284267 if (this -> hasDofID (R_w )) {
285- printf ("R_w is still active on node %d!\n" , this -> giveNumber ());
286268 FloatArray coeffs ;
287269 IntArray masterNodeIDs , masterDofIDs ;
288270
@@ -316,14 +298,6 @@ void HangingNode ::postInitialize()
316298 }
317299
318300 if (SlaveDof * sdof = dynamic_cast < SlaveDof * > (this -> giveDofWithID (R_w ))) {
319-
320- printf ("Slave DOF: node %d, dofID %d\n" , this -> giveNumber (), sdof -> giveDofID ());
321- for (int i = 1 ; i <= coeffs .giveSize (); ++ i ) {
322- printf (" Master %d: node %d, dofID %d, coeff = %.4e\n" ,
323- i , masterNodeIDs .at (i ), masterDofIDs .at (i ), coeffs .at (i ));
324- }
325-
326-
327301 sdof -> initialize (masterNodeIDs , masterDofIDs , coeffs );
328302 }
329303 }
0 commit comments