[deal.II] Re: a small issue with step44

True, the condition at the bottom should probably read AssertThrow > > > (newton_iteration < parameters.max_iterations_NR, > ExcMessage > > ("No convergence in nonlinear solver!")); Thanks, I'll make a note to fix this. On Thursday, November 17, 2016 at 1:22:54 AM UTC+1, Anup Basak

Re: [deal.II] issues of static executable in 8.4.1 (and latest one on github)

Yi-Chung, I try to compile static libraries by enabling static executable, but it failed. Can anyone help me about this? Here are some issues (and possible way to work around) 8.4.1 always get some error (looks like template instantiation problem) here are examples of errors --

Re: [deal.II] Compilation error

On 11/16/2016 07:40 PM, Lev Karatun wrote: [ 63%] Building CXX object source/base/CMakeFiles/obj_base.release.dir/mpi.cc.o /home/r/russ/lkaratun/aspect/dealii/source/base/mpi.cc: In constructor 'dealii::Utilities::MPI::MPI_InitFinalize::MPI_InitFinalize(int&, char**&,

[deal.II] Compilation error

Hi everyone, I was trying to update dealII to the latest version but got the following error: [ 63%] Building CXX object > source/base/CMakeFiles/obj_base.release.dir/mpi.cc.o > /home/r/russ/lkaratun/aspect/dealii/source/base/mpi.cc: In constructor >

[deal.II] Move vertices in P::D::Triangulation

Hello all, I am using parallel::distributed::triangulation, P4est, and Petsc to solve a 3D solid mechanics problem using deal.ii I am moving my mesh after every time step using the below written code. I have 2 questions regarding this. template void PlasticityContactProblem:: move_mesh

[deal.II] a small issue with step44

Dear developers, I see a small issue in function 'solve_nonlinear_timestep( BlockVector _delta)' of step44. If the 'for' loop corresponding to the Newton iteration is considered in its present form, the exception in the end

Re: [deal.II] Re: mesh quality during deformations using MappingQEulerian

On 11/16/2016 03:30 PM, thomas stephens wrote: but I'm not sure it will redistribute mesh points in tangential directions. I think there's a second step to this that moves vertices around. The important realization is that there are infinitely many ways of parametrizing the same surface.

Re: [deal.II] Re: Periodic boundary condition in Meshworker.

Thanks a lot J-P and Daniel, Let me please see how it works. Sudarshan On Wed, Nov 16, 2016 at 7:16 PM, Daniel Arndt < d.ar...@math.uni-goettingen.de> wrote: > Sudarshan, > > J-P is right. You can use constraints within the MeshWorker framework and > hence also periodic boundary

[deal.II] Re: mesh quality during deformations using MappingQEulerian

On Wednesday, November 16, 2016 at 4:40:05 PM UTC-5, Jean-Paul Pelteret wrote: > > Hi Tom, > > An alternative to implemented a complex mesh update strategy is to define > an initial mesh that, when mapped via your Euler vector, produces a higher > quality mesh. For example, you could remove >

[deal.II] issues of static executable in 8.4.1 (and latest one on github)

Hi I try to compile static libraries by enabling static executable, but it failed. Can anyone help me about this? Here are some issues (and possible way to work around) 8.4.1 always get some error (looks like template instantiation problem) here are examples of errors --

[deal.II] Re: mesh quality during deformations using MappingQEulerian

> > With Mesquite I believe one would essentially redefine a triangulation in > whatever format it wants (I think you provide nodal coordinate and > connectivities), tell it which nodes remain fixed (i.e. provide some > constraints) and then simply specify some algorithm with which it must >

Re: [deal.II] Re: mesh quality during deformations using MappingQEulerian

> > > In the case that the mesh lives in a higher dimensional space, you also > have to enforce -- either as part of the problem formulation, or as a > postprocess step for the output of Mesquite -- that the new nodes still > need to lie on the geometry as described before. In other words,

[deal.II] Re: Periodic Boundary Condition

Hamed, As you can see my assembly in the first post, my system_matrix and > system_rhs are not made directly from cell_system_matrix or > cell_system_rhs, but they are the result of some manipulations on > mass_matrix, laplase_matrix and nl_matrix (the same approach as step_25) > that have

[deal.II] Re: Periodic boundary condition in Meshworker.

Sudarshan, J-P is right. You can use constraints within the MeshWorker framework and hence also periodic boundary conditions. In more detail, you have to provide the assembler with the ConstraintMatrix. In addition to step-45 for periodic boundary conditions, you should also have a look to

[deal.II] Re: Periodic boundary condition in Meshworker.

Sudarshan, Does MeshWorker take a constraints matrix? I suspect that it does . If so then yes, you should be able to implement periodic boundary conditions. I

[deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Tom, An alternative to implemented a complex mesh update strategy is to define an initial mesh that, when mapped via your Euler vector, produces a higher quality mesh. For example, you could remove

[deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Tom, Yes, you're essentially correct. With Mesquite I believe one would essentially redefine a triangulation in whatever format it wants (I think you provide nodal coordinate and connectivities), tell it which nodes remain fixed (i.e. provide some constraints) and then simply specify some

[deal.II] Periodic boundary condition in Meshworker.

Is there any way to implement periodic boundary condition with MeshWorker,? Any help is appreciated. Thanks a lot, Best -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this

[deal.II] The "make run" command on the new version of deal.ii (8.4.1) complains about "libmpi.cxx.1.dylib" version.

I am trying to run the deal.II version 8.4.1 on a mac with OS X Yosemite version 10.10.5: cmake works fine as well make, but when I "make run" the code, it complains that the version of "libmpi.cxx.1.dylib" is 2.0.0 and the required version of "libmpi.cxx.1.dylib" must be 3.0.0 or above. I

[deal.II] Re: Periodic Boundary Condition

Dear Daniel, As you can see my assembly in the first post, my system_matrix and system_rhs are not made directly from cell_system_matrix or cell_system_rhs, but they are the result of some manipulations on mass_matrix, laplase_matrix and nl_matrix (the same approach as step_25) that have

[deal.II] Re: mesh quality during deformations using MappingQEulerian

J-P, I should not have given Mesquite such short shrift in my reply yesterday. With regard to the MappingQEulerian approach, I would imagine Mesquite would take the mapped points as the mesh points (basically, what my figures show as the mesh), then smooth those out a bit, and return new

Re: [deal.II] Re: question on doxygen

>> where did you need to change this? In deal.II somewhere? Or just in your own >> program? We already do this, see _mathjax_relpath in doc/doxygen/CMakeLists.txt -- Timo Heister http://www.math.clemson.edu/~heister/ -- The deal.II project is located at http://www.dealii.org/ For mailing

Re: [deal.II] Re: question on doxygen

Just on my own program. On Wednesday, November 16, 2016 at 9:44:59 AM UTC-5, Wolfgang Bangerth wrote: > > On 11/16/2016 07:41 AM, Michael Harmon wrote: > > > > I had to change the http from: > > > > MATHJAX_REPATH = http://cdn.mathjax.org/mathjax/latest/ > > > > to https: > > > >

Re: [deal.II] Re: question on doxygen

On 11/16/2016 07:41 AM, Michael Harmon wrote: I had to change the http from: MATHJAX_REPATH = http://cdn.mathjax.org/mathjax/latest/ to https: MATHJAX_REPATH = http*s*://cdn.mathjax.org/mathjax/latest/ Michael, where did you need to change this? In deal.II somewhere? Or just in your own

[deal.II] Re: hp HDG in deal.II

Samuel, I am attempting to implement a hp hybridized dg-method for solving > ellipctic problems in deal.II. I am relying on steps 27 and 51 of the > deal.II tutorial programs. However, when I try to group some FE_FaceQ > elements in an hp::FeCollection and distribute the degrees of freedom, I

[deal.II] Re: question on doxygen

Thanks Timo! I had to change the http from: MATHJAX_REPATH = http://cdn.mathjax.org/mathjax/latest/ to https: MATHJAX_REPATH = http*s*://cdn.mathjax.org/mathjax/latest/ On Tuesday, November 15, 2016 at 12:03:55 PM UTC-5, Michael Harmon wrote: > > Hi, > > I used deal for my phd work, loved

Re: [deal.II] Cluster with Infiniband network

On 11/16/2016 04:31 AM, Ashkan Dorostkar wrote: Could you elaborate on why I am using twice more memory on each node if I use more than one core per node? If you use two cores per node, then you are running two programs within the same amount of memory, each of which will allocate their own

[deal.II] hp HDG in deal.II

Hello all I am attempting to implement a hp hybridized dg-method for solving ellipctic problems in deal.II. I am relying on steps 27 and 51 of the deal.II tutorial programs. However, when I try to group some FE_FaceQ elements in an hp::FeCollection and distribute the degrees of freedom, I get