Manav Bhatia

Thanks, Roy. I am using h-adaptivity with quad/hex elements. So, I am guessing that I will need to use `enable-nodeconstraints` (?). I will give the PR a shot. If things...

By moving, do you mean change the nodal coordinate, or move it from rank i to rank j on the communicator? My workflow is: * initialize mesh/eq-sys/system * Add extra...

I see. I am not moving (nodal coordinates) any mesh points, but I am doing h-adaptivity. So if I understand correctly, I can turn off node-constraint in configuration? Out of...

Thanks Roy! Would using `semilocal_elem_begin/semilocal_elem_end` instead of `local_elem_begin/local_elem_end` include ghosted elements on each processor?

changing the iterator to `semilocal_elements_begin` gives the following errors: ``` libc++abi.dylib: Pure virtual function called! [InfiHorizon:49775] *** Process received signal *** libc++abi.dylib: Pure virtual function called! [InfiHorizon:49774] Signal: Abort trap:...

Looks like the issue was because I had written the loop like this: ``` MeshBase::element_iterator it = mesh.semilocal_elements_begin(), end = mesh.semilocal_elements_end(); for ( ; it != end; it++) mesh.delete_elem(*it); ```...

I have tried both smallest_magnitude and largest_magnitude with no success. I have played around with number of modes from 2 to 20.

chol(m0) and chol(m1) returns without errors on octave. However, m0 is symmetric, while m1 is not.

Here is a related thread on octave concerning the same issue: http://octave.1599824.n4.nabble.com/Issues-with-eigs-td4689314.html#a4689320 . While different people are able to converge to somewhat different eigenvalues on different releases of octave, I...

v 3.6.2 (installed from macports on Mac 10.13.6).