[deal.II] How to use sacado package

2016-10-21 Thread Jaekwang Kim
I'd like to use sacado package for automatic differentiation. However, I have installed deal.ii on my computer with trilinos configuration set off. I looked back deal.ii installation manual and got to know that I have to install trilnos package first. However, when I try to install

Re: [deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-21 Thread Hamed Babaei
Dear Daniel, Please look at the pdf file attached. Thanks On Friday, October 21, 2016 at 4:51:22 PM UTC-5, Daniel Arndt wrote: > > Hamed, > > At least I can't open your file in a way to clearly see what your > bilinearform is, but could just guess from the implementation. Can you > either

Re: [deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-21 Thread Daniel Arndt
Hamed, At least I can't open your file in a way to clearly see what your bilinearform is, but could just guess from the implementation. Can you either write it here (using LaTeX or similar) or upload the file as pdf again? Best, Daniel -- The deal.II project is located at

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

2016-10-21 Thread Oded Yaakobi
Hi Jean-Paul, I wanted to verify that the issue that I encountered is not related to the fact that I am using a version of dealii which is not the most updated one, so I installed the current version of dealii from git. I run two versions of my code with

Re: [deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-21 Thread Hamed Babaei
Dear Wolfgang, In attach you shall find the bilinear form and the assemble part of the code as well. Thanks, On Friday, October 21, 2016 at 11:14:52 AM UTC-5, Wolfgang Bangerth wrote: > > > > On Thursday, October 20, 2016 at 9:48:55 PM UTC-5, Wolfgang Bangerth > wrote: > > > > What is

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

2016-10-21 Thread Jean-Paul Pelteret
Hi Oded, I wasn't able to read in the .cub file (clearly its not backwards compatible). I did manage to import the .inp file and nothing looks amiss. Could you try to set the apply_all_indicators_to_manifolds

Re: [deal.II] question on contact algorithm for step-41

2016-10-21 Thread Wolfgang Bangerth
Hoang, I have some doubts over the contact algorithm used in step-41. Maybe someone can help me to clarify: + In step-41, instead of solving for delta_U and delta_Lambda, we solve for U_k+1 by using a finite difference equation. Is it some form of explicit contact algorithm? I'm not familiar

Re: [deal.II] Re: How to use Trilinos instead of Petsc in step-40

2016-10-21 Thread Hamed Babaei
Dear Wolfgang, On Thursday, October 20, 2016 at 9:48:55 PM UTC-5, Wolfgang Bangerth wrote: > > What is the norm of the first right hand side? You now set the tolerance > in > the first iteration to a fixed value of 1e-16, but how does this compare to > the previous value of

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

2016-10-21 Thread Jean-Paul Pelteret
Hi Oded, Two things: 1. The documentation originally stated that "Files generated using Cubit 11.x, 12.x and 13.x are valid", but I see that its now been updated to include 14.x and 15.x. I personally don't know whether this works with newer versions of Cubit, but hopefully it has indeed been

[deal.II] Re: How to apply boundary values for a particular point on the boundary instead of the whole boundary surface

2016-10-21 Thread Jean-Paul Pelteret
Dear Hamed, The problem here is that VectorTools::interpolate_boundary_values always applied boundary values to surfaces in 3d, or lines in 2d. You cannot directly apply point constraints using this function. What is happening here is that you've set a component mask which says that for ALL of