Re: [deal.II] Calculate cell center distance from a boundary

Drs. Bruno and Wolfgang, *Actually what Wolfgang is suggesting is very doable. We have wrappers for the ArborX library https://dealii.org/current/doxygen/deal.II/classArborXWrappers_1_1BVH.html which allows you to fi

[deal.II] Calculate cell center distance from a boundary

Dear all, Is there a way to get the shortest distance from cell center to a given boundary in p::d::Triangulation? What I really want is the wall normal distance. Any other suggestions are also welcome. Thanking in anticipation Vachan -- The deal.II project is located at http://www.dealii.org

Re: [deal.II] Centroids of a distributed triangulation

Thank you for the responses, Drs. Arndt and Bangerth. It greatly helps that a related function is available for a cell. Vachan: If you'd like to implement such a function, it would be great to > have that integrated into namespace GridTools! I am not entirely sure that these would be required,

[deal.II] Centroids of a distributed triangulation

Dear all, Is there a built-in function to calculate centroid of the area spanned by all the owned cells? I couldn't find any related function in GridTools. Using a bounding box will not work in all cases for my problem. Thanking in anticipation Vachan -- The deal.II project is located at http

[deal.II] Small suggestion to improve GridIn::read_unv()

Dear all, I found this thread which highlights one small issue in importing a unv file (generated by Salome): https://groups.google.com/g/dealii/c/iWnLv4v0V1M/m/38jm5Ck8WXYJ. If I understand correctly, dealii expects the file to start with "-1" and then have the section "2411", instead of all

Re: [deal.II] Use a coarse grid solution as initial condition for a finer grid

That brings a big simplification, at least when the aim is to just restart a simulation and not to map the solution onto some other grid. Thanks for all the help! On Tuesday, August 3, 2021 at 9:18:45 PM UTC+5:30 Wolfgang Bangerth wrote: > On 8/3/21 8:38 AM, vachan potluri wrote: > > > > I also

[deal.II] Can Triangulation::add_periodicity() be called multiple times?

Dear all, I have a rather silly question. Is it possible to call Triangulation::add_periodicity() multiple times with different matched pairs to give the same effect as to when the matched pairs contains *all* the periodic faces information? Or will a later call override the earlier settings?

[deal.II] Unexpected data output with cell data vector

Dear all, I am having some confusion regarding the index used in DataOut for a cell vector. I have created a distributed vector for storing cell data using p::d::Triangulation::global_active_cell_index_partitioner()

[deal.II] Use a coarse grid solution as initial condition for a finer grid

Dear all, I want to save a coarse grid solution, load it and set it as the initial condition for a finer grid simulation. The code I have written is for non-adaptive meshes, the refinement is not happening within the simulation. I wish to start the fine grid simulation using a standalone fine

[deal.II] dealii 9.3.0 make install fails at "Generating mpi.inst" with Invalid instantiation list: missing 'for'

Dear all, I am installing a newer version of dealii on our cluster (cray xc50) and make install fails at "Generating mpi.inst" with the statement Invalid instantiation list: missing 'for' make[2]: *** [source/base/CMakeFiles/obj_base_inst.dir/build.make:135: source/base/mpi.inst] Error 1 make[1]

Re: [deal.II] A data structure for distributed storage of some cell "average"

:) ! Thanks On Friday, June 11, 2021 at 11:40:37 AM UTC+5:30 vachanpo...@gmail.com wrote: > Hello, > > I am having an issue in using DataOut for such vector in a parallel > process. I am attaching a MWE which captures my problem. I am encountering > a segmentation fault (signal 11

Re: [deal.II] A data structure for distributed storage of some cell "average"

If I want to add such a vector to DataOut, will the regular DataOut::add_data_vector() work? Or is something else required to be done? On Saturday, May 29, 2021 at 8:44:01 AM UTC+5:30 vachanpo...@gmail.com wrote: > Ok, right. Will have a new installation then! > > On Sat, 29 May, 20

Re: [deal.II] Ordering of polynomials in FE_DGQLegendre<3>

Please have a look at pr12331 . On Thursday, May 27, 2021 at 6:56:58 PM UTC+5:30 Wolfgang Bangerth wrote: > On 5/25/21 11:26 PM, vachan potluri wrote: > > And the output shows all errors to be 0! If this is ok, I will create a > pr > > with a patch t

Re: [deal.II] A data structure for distributed storage of some cell "average"

other alternative (that I can use on 9.1.1)? Otherwise I will have to update the installation. On Thursday, May 27, 2021 at 7:01:12 PM UTC+5:30 Wolfgang Bangerth wrote: > On 5/26/21 10:20 PM, vachanpo...@gmail.com wrote: > > > > I have a quantity which has a single value in a ce

[deal.II] Re: Outputting DataOut Object into Personal Object Rather Than Outputing to File

Chen, I could not understand your motivation for doing this. If you want to "save" the data to be read again, then you can have a look at SolutionTransfer class. https://www.dealii.org/current/doxygen/deal.II/classSolutionTransfer.html Vachan On Wednesday, May 26, 2021 at 5:28:54 PM UTC+5:30

[deal.II] A data structure for distributed storage of some cell "average"

Dear All, I have a quantity which has a single value in a cell. Let me call this an "average". I want to be able to access the average of neighbouring cells from every cell. What data structure can I use? Since don't have mesh refinement, I was thinking of mapping this average to the cell ind

Re: [deal.II] Ordering of polynomials in FE_DGQLegendre<3>

be the same for FE_DGQLegendre. Apologies for my oversight. On Tuesday, May 25, 2021 at 10:37:24 AM UTC+5:30 vachanpo...@gmail.com wrote: > Wolfgang, > > Thanks a lot for your reply. What I actually need is a change of basis > from Lagrange polynomials (nodal) to Legendre polynomials

Re: [deal.II] Ordering of polynomials in FE_DGQLegendre<3>

.II/polynomial_8cc_source.html#l00877> (). I will also try your suggestion to figure out the ordering for a 3d element once I get some time. Thanks again! On Monday, May 24, 2021 at 9:02:29 PM UTC+5:30 Wolfgang Bangerth wrote: > On 5/24/21 5:36 AM, vachanpo...@gmail.com wrote: > > &g

[deal.II] Ordering of polynomials in FE_DGQLegendre<3>

Dear All, I wasn't able to figure out the ordering of polynomials stored in FE_DGQLegendre<3>. If I am correct, this class is constructed by FE_Poly's constructor using TensorProductPolynomials of Polynomials::Legendre. If this is indeed so, then I wasn't able to find the ordering mentioned in

Re: [deal.II] Matching dofs across cell at a common face for DG finite element

s not aware that this is such a difficult aspect to get right. The alternative you suggested will work for me. Thanks again for the reply! Vachan On Tuesday, May 11, 2021 at 7:29:41 PM UTC+5:30 Wolfgang Bangerth wrote: > On 5/10/21 4:53 AM, vachanpo...@gmail.com wrote: > > > > My

[deal.II] Matching dofs across cell at a common face for DG finite element

Dear all, I have a simple question. I am solving a DG problem locally and have maps established between cell-local and face-local dof numbering in a cell. Let me explain what I mean. For a 2nd order 3D DG element, I have the following maps constructed according to documentation of FE_DGQ

Re: [deal.II] make install fails with ld signal 9

David Wells > -- > *From:* dea...@googlegroups.com on behalf of > vachanpo...@gmail.com > *Sent:* Tuesday, May 4, 2021 2:45 AM > *To:* deal.II User Group > *Subject:* [deal.II] make install fails with ld signal 9 > > Dear all, > >

Re: [deal.II] Gmsh unexpected behaviour - axes directions swapped

I would probably always recommend using anyway) works for you. > > Best, > Daniel > > > > Am Fr., 7. Mai 2021 um 02:58 Uhr schrieb vachanpo...@gmail.com < > vachanpo...@gmail.com>: > >> I think I have found an explanation. >> >> Since the mesh generated

Re: [deal.II] Gmsh unexpected behaviour - axes directions swapped

erence axes directions, the output was unexpected to me. Any comments on whatever I concluded so far would be greatly helpful to me! On Thursday, May 6, 2021 at 10:41:21 AM UTC+5:30 vachanpo...@gmail.com wrote: > I have just read the gmsh manual and verified that the physical tags to > sur

Re: [deal.II] Gmsh unexpected behaviour - axes directions swapped

<< face_1->boundary_id() << "\n"; } if(face_2->at_boundary()){ std::cout << "Found face2 of local id " << 2*periodic_dir+1 << " at boundary with bid " << face_2->boundary_id() <<

[deal.II] Re: make install fails with ld signal 9

Bruno, As you suggested, I was able to successfully install the release version. Thanks! On Tuesday, May 4, 2021 at 5:54:57 PM UTC+5:30 bruno.t...@gmail.com wrote: > Vachan, > > > On Tuesday, May 4, 2021 at 2:45:30 AM UTC-4 vachanpo...@gmail.com wrote: > >> Dear all, >

[deal.II] Gmsh unexpected behaviour - axes directions swapped

Dear all, I have encountered a strange problem while dealing with physical boundary ids. I have a simple cube mesh (annexure 1) which I am exporting to msh4 format for reading through GridIn. The mesh is constructed by extruding a square 2D mesh in x-y plane to 3D. I was facing problems with u

[deal.II] make install fails with ld signal 9

Dear all, I am trying to do a new installation of dealii. `make install` exits while compiling libdeal_II.g.so with ld signal 9. I had a look at this question and noted that this might be due to RAM shortage, but my device has >2

Re: [deal.II] Unable to collect periodic faces

> triangulation after reading it in and modify the boundary indicators > manually if necessary. > > Best, > Daniel > > Am Di., 20. Apr. 2021 um 03:37 Uhr schrieb vachanpo...@gmail.com < > vachanpo...@gmail.com>: > >> Dear all, >> >> I have a si

[deal.II] Unable to collect periodic faces

Dear all, I have a simple cube gmsh file (annexed) which I want to read, and assign a periodic BC to it. I export the mesh to 'msh2' format and then read it using GridIn. I then call collect_periodic_faces()

Re: [deal.II] Setting manifold for a truncated cone mesh

This has worked! Many thanks! On Wednesday, April 7, 2021 at 9:11:04 PM UTC+5:30 peterrum wrote: > I am not sure. But it is worth a try! > > Peter > > On Wednesday, 7 April 2021 at 13:30:59 UTC+2 vachanpo...@gmail.com wrote: > >> Peter, >> >> Thanks for your

Re: [deal.II] Setting manifold for a truncated cone mesh

gt; > Hope this helps! > > Peter > > On Wednesday, 7 April 2021 at 12:46:04 UTC+2 vachanpo...@gmail.com wrote: > >> Sorry for continuing on this thread, my current question is not fully >> relevant, but related to this. >> >> I have a coarse mesh for a holl

Re: [deal.II] Setting manifold for a truncated cone mesh

Sorry for continuing on this thread, my current question is not fully relevant, but related to this. I have a coarse mesh for a hollow cylinder I have read this file and applied a cylindrical manifold by first calling Triangulation::set_all_manifold_ids() and then Triangulation::set_manifold()

Re: [deal.II] Setting manifold for a truncated cone mesh

_manifold_ids(0); > } > triangulation.set_manifold(0, CylindricalManifold<3>()); > > set_all_manifold_ids() makes sure that also the inner cells use the > manifold when refining. You should be able to copy the code above and > modify it according to your needs. > > Best, > Da

Re: [deal.II] Setting manifold for a truncated cone mesh

at should work. There is nothing special about using GridGenerator to > create geometries that can not be achieved by reading in a custom mesh and > setting the boundary ids, manifold ids and manifold. > > Best, > Daniel > > Am Di., 6. Apr. 2021 um 07:42 Uhr schrieb vach

[deal.II] Setting manifold for a truncated cone mesh

Dear all, I have a rather silly question; I am not confident about my understanding of manifolds. Suppose I want to replicate what GridGenerator::truncated_cone() does, but I have a

Re: [deal.II] Doubt in collect_periodic_faces() documentation

> > On 30. Mar 2021, at 14:22, Wolfgang Bangerth > wrote: > > On 3/30/21 11:20 AM, vachanpo...@gmail.com wrote: > > I would like to. However, I do not know how to write patches. Nor have I > any experience with using github collaboratively; I have only used it for >

Re: [deal.II] Doubt in collect_periodic_faces() documentation

Prof. Bangerth I would like to. However, I do not know how to write patches. Nor have I any experience with using github collaboratively; I have only used it for my personal repos. I would love to contribute if the experts can provide a little guidance. I have also noticed a small mistake in t

Re: [deal.II] Doubt in collect_periodic_faces() documentation

> > Best, > Daniel > > Am Fr., 26. März 2021 um 02:00 Uhr schrieb vachanpo...@gmail.com < > vachanpo...@gmail.com>: > >> Hi all, >> >> The documentation of collect_periodic_faces() >> <https://www.dealii.org/current/doxygen/deal.II/names

[deal.II] Doubt in collect_periodic_faces() documentation

Hi all, The documentation of collect_periodic_faces() reads Instead of defining a 'first' and 'second' boundary with the help of two boundary_ids this function defines a 'left' boundary a

[deal.II] Periodicity argument and domain in FunctionManifold

Dear All, I want to perform a 3D sector (possibly even a single strip of cells) simulation of flow over objects like blunt cone. As far as I am aware, dealii currently cannot directly read higher order meshes (I had a look at pr10163 and this discu

Re: [deal.II] Unable to visualise high order polynomial variation in ParaView

te: > >> Did you call >> >> data_out.build_patches(mapping, >> subdivisions,DataOut::curved_inner_cells); >> >> with subdivisions > 1? >> >> L. >> >> > On 20 Jan 2021, at 12:06, vachanpo...@gmail.com >> wrote: >> > >>

[deal.II] Unable to visualise high order polynomial variation in ParaView

Dear all, I have a working code that outputs pvtu/vtu files for visualisation. I have read the wiki page for writing and visualising high order output. Currently, I only have straight-edge quadrilaterals so I the o

Re: [deal.II] Build error | resizing a vector of FunctionParsers

rser instances) so that no > copies or moves are required. > > Best, > David > -- > *From:* dea...@googlegroups.com on behalf of > vachanpo...@gmail.com > *Sent:* Wednesday, January 13, 2021 1:23 AM > *To:* deal.II User Group > *Subject:*

Re: [deal.II] deal.ii Installation : expand_instantiations Error

Hi Aaditya and Johannes, My reply might be too late to be of help. About an year ago, I ran into a host of problems for installing dealii on Cray X50. I have posted all the details to get a working installation here https://groups.google.com/g/dealii/c/MCYyPrZNyjg/m/bVS1woSXBgAJ. The workaroun

[deal.II] Build error | resizing a vector of FunctionParsers

Dear all, I require to use a vector of an array of FunctionParsers in my code std::vector< std::array, 4> > fps_vec; At a later point when I know its size, I call fps_vec.resize(size); This line produces a long build error message (attached) starting with invalid application of ‘sizeof’ to i