Re: [deal.II] Re: Thermoelastic Problem

2019-08-05 Thread Muhammad Mashhood
Dear Prof. Wolfgang, Thank you for your concern and understanding. Yes you are right. I tried something similar instead of copying the triangulation. I used *const SmartPointer >* to transfer triangulation between both thermal and the solid mechanics classes

[deal.II] Reading mesh from vtk file

2019-08-05 Thread Andreas Rupp
Dear all, I am trying to read a grid from a vtk file created by Paraview. Unfortunately, this does not work. Moreover, reading a vtk file that has been created using deal.ii also does not work for me. It seems to me that deal.ii reads the file but does not connect the faces to form a mesh.

Re: [deal.II] How to efficiently add / remove cells (maybe with hp::DoFHandler)

2019-08-05 Thread Mathias Anselmann
Dear Professor Bangerth, thanks again for your reply, I will start answering your questions and give you more insight into what I try to do at the moment: That's exactly the question you have to answer first. What do you *want* to > happen on these cells? If you know what you want to happen,

[deal.II] DG Matrix-free: getting face iterator form

2019-08-05 Thread Michał Wichrowski
Dear all, I'm trying to use matrix-free with discontinuous Galerkin elements for variable coefficient problem. The face integral require access to value of coefficient. For the cell I've done the following thing: //Vectors of values, stored in object: AlignedVector< VectorizedArray > mu;

Re: [deal.II] Reading mesh from vtk file

2019-08-05 Thread Andreas Rupp
Dear Luca, Thank you very much for your quick reply. However, I have to admit that I do not fully understand it. The vtk is created using GridOut grid_writer; std::ofstream output_file("./DealOutput.vtk"); grid_writer.write_vtk(triangulationSelf, output_file); Thus, I assume that

Re: [deal.II] Reading mesh from vtk file

2019-08-05 Thread Luca Heltai
Dear Andreas, You are trying to read a vtk file that was generated by DataOut. Try to write the file with GridOut::write_vtk. DataOut produces files that are split cell wise (to allow for some flexibility in discontinuos fields. Luca > Il giorno 5 ago 2019, alle ore 03:05, Andreas Rupp

[deal.II] deal.II Newsletter #89

2019-08-05 Thread Rene Gassmoeller
Hello everyone! This is deal.II newsletter #89. 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: #8456: Further cleanups of the SVG GridOut code. (proposed by bangerth)

Re: [deal.II] deal.II Newsletter #89

2019-08-05 Thread Wolfgang Bangerth
Ah, we're back on track with a substantial number of pull requests. It will be interesting to see how much we get done during the deal.II workshop this week -- I'm hoping for a very long list! W. On 8/5/19 12:03 PM, Rene Gassmoeller wrote: > Hello everyone! > > This is deal.II newsletter

[deal.II] Postdoc Position (deal.II related), International School for Advanced Studies, Trieste, Italy

2019-08-05 Thread luca.heltai
Dear all, I’d like to advertise an opening for a Postdoctoral position (18 months) in my group, to work with and on deal.II. The deadline is pretty tight… please send me a private message if you are interested. Luca. --- Title of the position: Adaptive methods for solving multi-physics

Re: [deal.II] constraining dofs across distributed mesh

2019-08-05 Thread Reza Rastak
Thank you Danial for the explanation. I used the IndexSet::add_index() method and its seems to be working so far. Reza On Wednesday, July 10, 2019 at 8:42:51 PM UTC-7, Daniel Arndt wrote: > > > [...] >> It works in the serial case. However, it doesn't work when I use multiple >> MPI processes.

[deal.II] Re: Dealii Installtion fails because it uses the wrong MPI Version of the intel compiler. It uses debug_mt instead of release_mt.

2019-08-05 Thread 张嘉宁
Hi, recently, I have the same problem. When I complie the PETsc, it link the debug_mt/linmpi.so. So I think the reason is that Petsc is complied in Debugging, not Release. So I add the paras when configure the Petsc,--with-debugging=0 and give the opt-paras, COPTFLAG="-O2" CXXOPTFLAG="-O2"

Re: [deal.II] Re: Serendipity elements

2019-08-05 Thread Jonathan Russ
Professor Bangerth - Okay thank you for your advice! I'll take a look at the FE_Poly class as you suggest. Thanks again, Jonathan On Sun, Aug 4, 2019 at 9:14 PM Wolfgang Bangerth wrote: > On 8/4/19 5:47 PM, Jonathan Russ wrote: > > > > Thank you for your reply. The shape functions of the