[deal.II] A small question about subdomain renumbering

Dear all, I have a question about DoFRenumbering::subdomain_wise(dof_handler). In step-17 tutorial for example, this function doesn't use the argument of n_mpi_processes, how does it know how many subdomains it should divide the dofs into? After we renumbered the subdomain and get the locally_own

Re: [deal.II] Re: Set up boundary id in gmsh

The pleasure is mine. You will find that people the deal.II community are very helpfull :)! On Wednesday, 5 August 2020 22:44:11 UTC-4, yuesu jin wrote: > > Dear Bruno, > Thank you very much!!! This solved my problem. This is very helpful! > Best regards, > Yuesu > > On Wed, Aug 5, 2020 at 9:3

Re: [deal.II] Re: Set up boundary id in gmsh

Dear Bruno, Thank you very much!!! This solved my problem. This is very helpful! Best regards, Yuesu On Wed, Aug 5, 2020 at 9:33 PM Bruno Blais wrote: > The best way to proceed is to set the boundary conditions using Physical > Entities in GMSH. > An example is always better. For example, the

[deal.II] Re: Set up boundary id in gmsh

The best way to proceed is to set the boundary conditions using Physical Entities in GMSH. An example is always better. For example, the following GMSH code : lc = 2.0e-1; lf = 2.0e-1; RO=1; RI=0.25; Point(0) = {0, 0, 0, lc}; Point(1) = {RO, 0, 0, lc}; Point(2) = {0, -RO , 0, lc}; Point

[deal.II] Set up boundary id in gmsh

Dear all, I want to set up one Dirchlet boundary condition and one Neumann boundary condition on a 2d mesh, which is generated by gmsh. How can I set up the mesh file in order to tell dealii with the boundary id when I input the .msh file into dealii? Thanks! Best regards, -- Yuesu Jin, Ph.D st

Re: [deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

Hi, Wolfgang, Thank you for your reply. The complete code can be found on GitHub and we also have the instruction for compiling the code on TACC in UTA. This code is in the branch thm_seg_distributed. You may like to download it and build it. https://github.com/cb-geo/thm/tree/thm_seg_distrib

Re: [deal.II] Output of Gauss point stress tensor

Thanks for the guidance. I tried replacing the " source/particles/particle_handler.cc

Re: [deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

On 8/5/20 10:42 AM, 孙翔 wrote: Please see the attached file, which is the head file that used the library. Thank you very much. That's not enough -- you're using a file "interpolation.h" that I don't have, and parts of the rest of that file depend on this. You need to give me something that'

Re: [deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

Hi, Wolfgang, Please see the attached file, which is the head file that used the library. Thank you very much. Best, Xiang On Wednesday, 5 August 2020 09:14:51 UTC-7, Wolfgang Bangerth wrote: > > On 8/5/20 10:13 AM, 孙翔 wrote: > > > > I have included it actually, and also the distributed/grid

Re: [deal.II] Output of Gauss point stress tensor

Dear Prof. Wolfgang, Thanks for the guidance. I tried replacing the " source/particles/particle_handler.cc " file with the one present in "https://github.com/dealii/dea

Re: [deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

On 8/5/20 10:13 AM, 孙翔 wrote: I have included it actually, and also the distributed/grid_refinement.h was also included. Because I want to import the outside mesh to build up the triangulation, I also used Gridin to work on the triangulation, which is like as the following: Then I have no o

Re: [deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

Hi, Wolfgang, I have included it actually, and also the distributed/grid_refinement.h was also included. Because I want to import the outside mesh to build up the triangulation, I also used Gridin to work on the triangulation, which is like as the following: void CoupledTH::make_grid_and_dofs(

Re: [deal.II] Reading a Tensor from parameter file

> Thank you. This really shortens the code. Is it also possible to set the > delimiter to a user desired value, say "|" instead of ",". All functions of the Convert class take a unique pointer to a pattern. You can specify your own pattern, and use that in the conversion. The code is all in de

Re: [deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

On 8/5/20 3:01 AM, 孙翔 wrote: error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’ I suspect that you forgot in #include Best W. -- Wo

Re: [deal.II] Reading a Tensor from parameter file

Dear Luca, Thank you. This really shortens the code. Is it also possible to set the delimiter to a user desired value, say "|" instead of ",". Also, where exactly is the procedure for conversion from string to deali::Tensor or some other type say std::vector defined? I could not find it in the

[deal.II] Re: Getting started and posting guidelines for new users

Dear Kaleem, It would be good to open your own topic with your question so that you may get help. Best Bruno On Wednesday, 5 August 2020 05:30:44 UTC-4, kaleem iqbal wrote: > > The following error shown in step-49 during make run. > Exception on processing: > > -

[deal.II] Doctoral Programme in Mathematical Analysis, Modelling and Applications at SISSA, Trieste, Italy

Dear All, I would like to bring at your attention that the second (and last) deadline to apply for a Phd position in Mathematical Analysis, Modelling and Applications at SISSA, Scuola Internazionale Superiore di Studi Avanzati Trieste, Italy is August 20, 2020 at noon (Rome time). Research ac

Re: [deal.II] Reading a Tensor from parameter file

> Could you please provide an MWE to describe how the > Patterns::Convert::to_value() function would work in this case. In the dealii/tests/parameter_handler/patterns_05.cc there are many examples that show how to use these. > Is it must to use prm.add_parameter() to be able to do so? > I usua

[deal.II] Re: Getting started and posting guidelines for new users

The following error shown in step-49 during make run. Exception on processing: An error occurred in line <1430> of file in function void dealii::GridIn::read_msh(std::istream&) [with int dim = 2; int spacedim = 2; std::istream = std

[deal.II] error: cannot declare field ‘CoupledTH<3>::triangulation’ to be of abstract type ‘dealii::parallel::distributed::Triangulation<3, 3>’

Hi, I want to declare a variable triangulation as the type dealii::parallel::distributed::Triangulation<3, 3>. The code is : in main function, using namespace dealii; Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); CoupledTH<3> coupled_TH_solver(1); coupled_TH_solver.run(

Re: [deal.II] ERROR WHEN RUNNING THE CODE ON MULTIPLE NODES OF HPC

Got it. Thanks a lot. On Monday, 3 August 2020 06:10:40 UTC-7, Wolfgang Bangerth wrote: > > On 8/2/20 11:26 PM, 孙翔 wrote: > > > > Yes, it cannot run on a cluster. Both of them run in release mode. I'm > also > > curious about the error. I debug the code by outputting some specific > values.