[deal.II] Re: Hanging node constraints and periodic constraints together causing an issue

2018-01-16 Thread Sambit Das 
Thank you, Dr. Arndt.

Best,
Sambit

On Tuesday, January 16, 2018 at 11:16:08 AM UTC-6, Daniel Arndt wrote:
>
> Sambit,
>
> I created an issue at https://github.com/dealii/dealii/issues/5741 with a 
> modified version of your example.
>
> Best,
> Daniel
>

-- 
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 message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] step-42 clarification

2018-01-16 Thread Wolfgang Bangerth 

On 01/16/2018 12:59 AM, Alberto Salvadori wrote:

Dear community,

I am studying with pleasure step-42 and I got a bit confused in the 
method PlasticityContactProblem::solve_newton. At the very end of it, we 
find these instructions:


old_active_set = active_set;
previous_residual_norm = residual_norm;
if (Utilities::MPI::sum 
((active_set 
== old_active_set) ? 0 : 1,

mpi_communicator) == 0)
{
pcout << " Active set did not change!" << std::endl;
if (residual_norm < 1e-10)
break;
}


Now, I am a bit confused by the if statement and I am clearly missing 
something. The first instruction makes old_active_set equal to 
active_set. Accordingly, wouldn't active_set == old_active_set  always 
be true? What am I missing?


You're right -- the code would suggest that the if-condition is always 
true and that what matter is only the inner if.


I'm not sure I have the current email address of Joerg Frohne who wrote 
step-42. He would be the one who could possibly provide a better answer.


Best
 W.

--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 message because you are subscribed to the Google Groups "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] Nedelec Elements and non-tangential Dirichlet data

2018-01-16 Thread Wolfgang Bangerth 

On 01/16/2018 01:41 AM, Pascal Kraft wrote:
I am currently using a FeSystem composed of two 3D Fields (real and 
imaginary E-field) and I want to impose Dirichlet conditions of the sort 
E(x,y,z) = E_{in}(x,y,z) on the input interface (in an xy-plane). 
Earlier the z-component had been 0 so I did not run into real problems 
and could use the project_boundary_values_curl_conforming_l2 function 
(for the x and y components) and my own code for the z-components 
(setting them 0 in the first cell layer). Now however I would like to 
impose different values on the z-component. Since my cells have variable 
thickness in z-direction (i.e. hanging node constraints appear) I would 
like to use some library function for the interpolation.

Are there any feasible solutions to this?
Ideally I would like to keep using 
project_boundary_values_curl_conforming_l2 for the tangential components 
and just add some code for the z-components.


The question doesn't quite have enough information:
* Are you using Nedelec elements for the two fields?
* If the input face is in the xy-plane, and you are using a finite 
element that is only curl-conforming (such as the Nedelec element), 
isn't it correct that you can only prescribe *tangential* components for 
boundary values? In other words, you cannot prescribe z-components anyway.


Best
 W.

--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 message because you are subscribed to the Google Groups "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[deal.II] Re: Hanging node constraints and periodic constraints together causing an issue

2018-01-16 Thread Daniel Arndt 
Sambit,

I created an issue at https://github.com/dealii/dealii/issues/5741 with a 
modified version of your example.

Best,
Daniel

-- 
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 message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[deal.II] deal.II Newsletter #17

2018-01-16 Thread Rene Gassmoeller 
Hello everyone!

This is deal.II newsletter #17.
It automatically reports recently merged features and discussions about the 
deal.II finite element library.


## Below you find a list of recently proposed or merged features:

#5738: Fix a grammar error. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/5738

#5737: Fix warnings: do not use Assert in CUDA kernel (proposed by Rombur) 
https://github.com/dealii/dealii/pull/5737

#5736: Change .cuh extension to h (proposed by Rombur) 
https://github.com/dealii/dealii/pull/5736

#5735: Add cuSPARSE support in the build system (proposed by Rombur) 
https://github.com/dealii/dealii/pull/5735

#5734: Extend scalapack more (proposed by BenBrands) 
https://github.com/dealii/dealii/pull/5734

#5733: Minor change in FE Enriched element (proposed by niveshd; merged) 
https://github.com/dealii/dealii/pull/5733

#5732: Follow-up of #5655: add mmult and Tmmult to petsc sparse matrix and 
petsc parallel sparse matrix (proposed by Jie-Cheng) 
https://github.com/dealii/dealii/pull/5732

#5731: Increase precision in FE_RT_Bubbles tests (proposed by masterleinad) 
https://github.com/dealii/dealii/pull/5731

#5730: Make some variables const (references) (proposed by masterleinad; 
merged) https://github.com/dealii/dealii/pull/5730

#5729: Use nullptr instead of 0 or NULL (proposed by masterleinad; merged) 
https://github.com/dealii/dealii/pull/5729

#5728: Be explicit about virtual functions called in con-/destructors (proposed 
by masterleinad) https://github.com/dealii/dealii/pull/5728

#5726: Use automatic memory management, rather than doing it by hand. (proposed 
by bangerth; merged) https://github.com/dealii/dealii/pull/5726

#5725: Fix a formula that trips up doxygen. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/5725

#5724: Fix grammar in parameter description (proposed by bartniday; merged) 
https://github.com/dealii/dealii/pull/5724

#5723: Use complete grammar. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/5723

#5722: New Finite Element RT+Bubbles (PR #4540 part II) (proposed by eldarkh; 
merged) https://github.com/dealii/dealii/pull/5722

#5721: Add another test. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/5721

#5720: Fix bug in MappingFEField. (proposed by luca-heltai; merged) 
https://github.com/dealii/dealii/pull/5720

#5719: Update some comments around renumbering DoFs. (proposed by bangerth; 
merged) https://github.com/dealii/dealii/pull/5719

#5718: DoFRenumbering::hierarchical() for parallel::shared::Triangulation 
(proposed by bangerth) https://github.com/dealii/dealii/pull/5718

#5717: Make sure distribute_dofs doesn't copy unnecessarily (proposed by 
masterleinad; merged) https://github.com/dealii/dealii/pull/5717

#5714:  Replace HyperShellBoundary by SphericalManifold in tests (proposed by 
masterleinad) https://github.com/dealii/dealii/pull/5714

#5713: suppress PHG warning saying: PHG_EDGE_SIZE_THRESHOLD is low (proposed by 
vishalkenchan; merged) https://github.com/dealii/dealii/pull/5713

#5712: Make the instructions for the p4est installation script more precise 
(proposed by masterleinad; merged) https://github.com/dealii/dealii/pull/5712

#5711: Simplify implementation for normal_vector and get_normals_at_vertices in 
SphericalManifold (proposed by masterleinad) 
https://github.com/dealii/dealii/pull/5711

#5710: Update documentation by using Manifold instead of Boundary (proposed by 
masterleinad; merged) https://github.com/dealii/dealii/pull/5710

#5709: New polynomial space RT+Bubbles (PR #4540 part I) (proposed by eldarkh; 
merged) https://github.com/dealii/dealii/pull/5709

#5708: fix Zoltan debug output (proposed by davydden; merged) 
https://github.com/dealii/dealii/pull/5708

#5706: Update output of bits/cone_03 (proposed by masterleinad; merged) 
https://github.com/dealii/dealii/pull/5706

#5705: Extended step-26 documentation with regard to positivity preservation. 
(proposed by marcfehling; merged) https://github.com/dealii/dealii/pull/5705

#5704:  add template wrappers for ScaLAPACK functions (proposed by BenBrands; 
merged) https://github.com/dealii/dealii/pull/5704

#5703:  Replace HyperShellBoundary by SphericalManifold in tests (proposed by 
masterleinad; merged) https://github.com/dealii/dealii/pull/5703

#5702: ROL: bug fix applyBinary function in Rol::VectorAdaptor (proposed by 
vishalkenchan; merged) https://github.com/dealii/dealii/pull/5702

#5701: Fix failing tests after #5691 (proposed by masterleinad; merged) 
https://github.com/dealii/dealii/pull/5701

#5700: Make AssertCuda use Assert instead of AssertThrow (proposed by Rombur; 
merged) https://github.com/dealii/dealii/pull/5700

#5699: Do not use pedantic when cuda support is enabled (proposed by Rombur; 
merged) https://github.com/dealii/dealii/pull/5699

#5698: Rebase Opencascade Boundary classes (proposed by masterleinad; merged)

[deal.II] Nedelec Elements and non-tangential Dirichlet data

2018-01-16 Thread Pascal Kraft 
I am currently using a FeSystem composed of two 3D Fields (real and 
imaginary E-field) and I want to impose Dirichlet conditions of the sort 
E(x,y,z) = E_{in}(x,y,z) on the input interface (in an xy-plane). Earlier 
the z-component had been 0 so I did not run into real problems and could 
use the project_boundary_values_curl_conforming_l2 function (for the x and 
y components) and my own code for the z-components (setting them 0 in the 
first cell layer). Now however I would like to impose different values on 
the z-component. Since my cells have variable thickness in z-direction 
(i.e. hanging node constraints appear) I would like to use some library 
function for the interpolation.
Are there any feasible solutions to this? 
Ideally I would like to keep using 
project_boundary_values_curl_conforming_l2 for the tangential components 
and just add some code for the z-components.

-- 
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 message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[deal.II] step-42 clarification

2018-01-16 Thread Alberto Salvadori 
Dear community,

I am studying with pleasure step-42 and I got a bit confused in the method 
PlasticityContactProblem::solve_newton. At the very end of it, we find 
these instructions:

old_active_set = active_set;
previous_residual_norm = residual_norm;
 
if (Utilities::MPI::sum 
((active_set
 
== old_active_set) ? 0 : 1,
mpi_communicator) == 0)
{
pcout << " Active set did not change!" << std::endl;
if (residual_norm < 1e-10)
break;
}


Now, I am a bit confused by the if statement and I am clearly missing 
something. The first instruction makes old_active_set equal to active_set. 
Accordingly, wouldn't active_set == old_active_set  always be true? What am 
I missing?

Thanks for your help

Alberto

-- 

Informativa sulla Privacy: http://www.unibs.it/node/8155

-- 
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 message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.