[deal.II] Nonhomogeneous Dirichlet Boundary conditions using a Dirichlet lift

2017-02-08 Thread Giulia Deolmi
Dear deal.ii users, is there someone who implemented Nonhomogeneous Dirichlet Boundary conditions using a Dirichlet lift? Thanks a lot in advance, Kind regards, Giulia -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/

[deal.II] ParameterHandler: MultipleSelection of couples

2017-02-08 Thread Franco Milicchio
Dear all, I am using the ParameterHandler class to select multiple parameters and lists, but I don't know if I could have something like the following. What I'd like is to specify multiple parameters that, sort of nested, have a string and an integer. As a concrete example, I want to specify st

Re: [deal.II] Nonhomogeneous Dirichlet Boundary conditions using a Dirichlet lift

2017-02-08 Thread Praveen C
Hello Giulia The usual way of applying Dirichlet bc in deal.II essentially does a lifting approach. If u = g on boundary then the lifting is u_{g,h}(x) = sum_(i on boundary) g(x_i) \phi_i(x) Did you want to use a different lifting ? Best praveen On Wed, Feb 8, 2017 at 8:38 PM, Giulia Deolmi

Re: [deal.II] ParameterHandler: MultipleSelection of couples

2017-02-08 Thread Praveen C
Hi In a similar case I have used this function http://www.dealii.org/developer/doxygen/deal.II/namespaceUtilities.html#a8d799bb35ac16d206818c88e82afbfae You can use this twice to get the result you want, combined with http://www.dealii.org/developer/doxygen/deal.II/namespaceUtilities.html#af501

Re: [deal.II] ParameterHandler: MultipleSelection of couples

2017-02-08 Thread Franco Milicchio
Thanks, Praveen. I could use those, and in my tests that I've done now they work great, but this means I won't be able to detect errors automatically. This deal.II feature is the best thing, dropping it would be not so cool. If not possible, I will just do this by hand. Thank you! Franco

Re: [deal.II] Fully distributed triangulation (level 0)

2017-02-08 Thread Timo Heister
> Are you willing to share that code, Timo? Yes, we will be creating a PR for that soon. > I suspect that if implemented right, it should not be terribly difficult to > do refinement of the mesh, but because you can't repartition the coarse > mesh, it will quickly become unbalanced if processors

[deal.II] Assemble Righthand Side for vector-valued problem

2017-02-08 Thread Jaekwang Kim
Hi all, I was testing my code, with method of manufactured solution. I have a source term, which is vector value, defined as .. template class RightHandSide : public Function { public: RightHandSide () : Function(dim+1) {} virtual void vector_value (const Point &p,

Re: [deal.II] Assemble Righthand Side for vector-valued problem

2017-02-08 Thread Wolfgang Bangerth
On 02/08/2017 12:35 PM, Jaekwang Kim wrote: */_local_rhs(i) += fe_values.shape_value(i,q) *_/* */_rhs_values[q] *_/* */_fe_values.JxW(q);_/* rhs_values[q] is a Vector. You need to say which component of it you wa

Re: [deal.II] Assemble Righthand Side for vector-valued problem

2017-02-08 Thread Jaekwang Kim
Thanks, I understand what you meant!! I fixed it. Add up this.. can I ask your intuition for my another problem? I am using deal-ii to analyze Non-newtonian fluids, which varies viscosity. For example, my viscosity eta is function of (u) field. I am solving this with iteration method.. (Picard

Re: [deal.II] Assemble Righthand Side for vector-valued problem

2017-02-08 Thread Wolfgang Bangerth
On 02/08/2017 01:15 PM, Jaekwang Kim wrote: I am using deal-ii to analyze Non-newtonian fluids, which varies viscosity. For example, my viscosity eta is function of (u) field. I am solving this with iteration method.. (Picard) After coding, I am testing my code with method of manufacture. For s