[deal.II] Problems in handling the userdata of Distributed Triangulation

I add the userdata to each cell using the cell->user_pointer. It seems that the userdata of some cells which are changed during the refining and coarsening procedure will lose. How can make sure that the refined children cells(the coarsened parent cell) inherit the userdata from the parent (t

[deal.II] Re: Problems in handling the userdata of Distributed Triangulation

Dear Jean， Thanks so much for your suggestions. I’m really very happy for such new development. I understand that using the CellDataStorage and TransferableQuadraturePointData classes is better than user_pointer when the mesh requires coarsening and refining. *Yet I do not understan

[deal.II] Re: Problems in handling the userdata of Distributed Triangulation

Dear Denis, Thanks for your replies. I do not have some data attached to ghost cells. But I have a variable that changes with time and its change depends on the data of neighbor cells. For example, const double value1 = 1, value2 =-1; for (cell = tr.begin_active(); cell != tr.end();

Re: [deal.II] Re: Problems in handling the userdata of Distributed Triangulation

. If I figure out this problem I will post it to the forum. If you have made some progress on this problem, please tell me promptly. Thanks so much. Cheers, Jack -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d

Re: [deal.II] Re: Problems in handling the userdata of Distributed Triangulation

> > Dear Prof. Bangerth, > Many thanks for your suggestions and help! > I will try this with your reference of Aspect. I hope I can figure out this > smoothly. > Thanks again! > Regards, > Jack > 在 2016年11月3日星期四 UTC+8下午10:54:44，Wolfgang Bangerth写道： > > On 1

Re: [deal.II] Re: Problems in handling the userdata of Distributed Triangulation

break; } } catch (typename Mapping::ExcTransformationFailed &) {} } } #endif const types::LevelInd found_cell = std::make_pair(cell->level(),cell->index()

Re: [deal.II] Re: Problems in handling the userdata of Distributed Triangulation

在 2016年11月9日星期三 UTC+8上午10:06:07，Jack写道： > > Hi, everyone > > I’m not familiar to ASPECT, could someone tell me how to get the data of > ghost cells? > In the ASPECT, I find that it is implemented through the following > functions to obtain ghost particles > in ~/a

[deal.II] Re: New: step-57

ut I did not get any good idea from the tutorial > examples (step-22, step-32, etc.) for Stokes flow. > > Thanks very much for your help! > > > > Jack > 在 2017年2月11日星期六 UTC+8上午12:47:34，Timo Heister写道： > > I am happy to announce that we added a new tutorial program:

[deal.II] Re: New: step-57

e tried to extend them to solve Navier-Stokes equations by adding the convection term into the left hand matrix but always failed. Thanks very much! Regards, Jack 在 2017年3月26日星期日 UTC+8下午5:45:48，Daniel Arndt写道： > > Jack, > > Is it possible to extend your code to time depend

[deal.II] Re: New: step-57

e tried to extend them to solve Navier-Stokes equations by adding the convection term into the left hand matrix but always failed. Thanks very much! Regards, Jack [1] https://en.wikipedia.org/wiki/Stokes_flow 在 2017年3月26日星期日 UTC+8下午5:45:48，Daniel Arndt写道： > > Jack, > >

Re: [deal.II] Re: New: step-57

Hi Timo, Thanks for your help. I have presented the formulations for time-dependent NS equations in the way you have done in step-57. Is it reasonable or you have a better suggestion? Thanks very much for your advices and help.

[deal.II] Re: Crack propagation

Hi S. A. Mohseni, I'm also interested on the cracking mechanics, where is such an example for crack propagation ? Thanks ! 在 2017年3月30日星期四 UTC+8下午9:13:16，Seyed Ali Mohseni写道： > > Dear Timo Heister and Thomas Wick, > > I am trying to run your phase-field crack propagation example, namely the

[deal.II] Re: Crack propagation

Hi Daniel, Thank you very much! Regards, Jack 在 2017年5月14日星期日 UTC+8下午5:07:26，Daniel Arndt写道： > > Jack, > > have a look at https://github.com/tjhei/cracks > > Best, > Daniel > > Am Freitag, 12. Mai 2017 07:57:06 UTC+2 schrieb Jack: >> >> Hi S. A. Mo

[deal.II] It seems that parallel::distributed::SolutionTransfer does not support hp::DoFHandler, how can I transfer solution after mesh changed

? Many thanks! Best, Jack -- 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"

[deal.II] Re: It seems that parallel::distributed::SolutionTransfer does not support hp::DoFHandler, how can I transfer solution after mesh changed

Hi Jean-Paul, Thanks for your help. My program failed to compile, so I didn't run it and get such an assertion. I see that the un-supporting of hp::DoFHandler for parallel distributed calculation in hp namespace. Thanks! Regards, Jack 在 2017年5月15日星期一 UTC+8下午10:27:22，Jean-Paul Pelte

[deal.II] Problems in solving two PDE systems using MPI+Thread parallel computing

you very much. Best, Jack -- 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.

[deal.II] Re: Problems in solving two PDE systems using MPI+Thread parallel computing

d such problem will be avoided and I do not need to change my code for initialization of MPI and TBB myself? Thanks for your kind help! Regards, Jack 在 2017年5月19日星期五 UTC+8下午8:25:55，Bruno Turcksin写道： > > Jack, > > are your solvers using MPI? This looks similar to this problem &g

Re: [deal.II] Re: Problems in solving two PDE systems using MPI+Thread parallel computing

rsion, it does not has the libs of libmpi_usempi.so which is ignored since v2.0. Thank you for your help Regards, Jack 在 2017年5月21日星期日 UTC+8上午1:33:51，Bruno Turcksin写道： > > Jack, > > 2017-05-19 22:06 GMT-04:00 Jack >: > > I should update my OpenMPI to the latest version a

Re: [deal.II] Re: Problems in solving two PDE systems using MPI+Thread parallel computing

the way is it possible to copy Vector and Matrix of TrilinosWrappers to another ones those are initialized with another MPI communicator? Thanks, All best， Jack 在 2017年5月22日星期一 UTC+8下午9:50:32，Wolfgang Bangerth写道： > > > Jack, > unrelated to the question about versions of M

Re: [deal.II] Re: Problems in solving two PDE systems using MPI+Thread parallel computing

computing time. I’m not going to try this idea at present. When I’m free, I will test its feasibility and post here the result. Thank you very much! All my best wishes, Jack 在 2017年5月24日星期三 UTC+8上午1:59:17，Wolfgang Bangerth写道： > > > Jack, > > > “The way to do this i

Re: [deal.II] Re: Problems in solving two PDE systems using MPI+Thread parallel computing

！ 在 2017年5月24日星期三 UTC+8下午10:38:57，Wolfgang Bangerth写道： > > On 05/23/2017 07:24 PM, Jack wrote: > > Originally, I suspected that when solving two linear systems > simultaneously by > > two threads would reduce time. But now it seems that this idea increases > the > &g

[deal.II] 5120 Segmentation fault (core dumped)

I am compiling dealII and get the error: /bin/sh: line 1: 5120 Segmentation fault (core dumped) ../../bin/expand_instantiations /home/jurombo/Downloads/build/share/deal.II/template-arguments < /home/jurombo/Downloads/dealii/source/opencascade/boundary_lib.inst.in > /home/jurombo/Download

Re: [deal.II] Docker and Eclipse

Managed to install from the latest git repo and it went perfectly well: MathW15Lnx:/home/jurombo/Downloads/dealii/build # make test Running quicktests... Test project /home/jurombo/Downloads/dealii/build/tests/quick_tests Start 1: step.debug Start 2: step.release Start 3: a

[deal.II] MPI and Docker

The current docker image is not finding mpi.h *[ 50%] Building CXX object CMakeFiles/cdr.dir/cdr.cc.o In file included from /usr/include/deal.II/base/conditional_ostream.h:19, from /home/dealii/dealii/cdr/cdr.cc:1: /usr/include/deal.II/base/config.h:490:12: fatal error

Re: [deal.II] step-35 in modified

Still working on modified step-35, but now starting the modification by Muhammad Mohebujjaman. https://github.com/Mohebujjaman/Projection-Method-for-NSE I want to resolve the issue of the iterators defined here: *template voidNavierStokesProjection::assemble_one_cell_of_gra

Re: [deal.II] step-35 in modified

Does this replace all the four instances of *std::get *or just the last two? On Friday, 4 March 2022 at 22:39:14 UTC+2 Wolfgang Bangerth wrote: > On 3/4/22 12:12, jack urombo wrote: > > / > > std::get<0> (SI.iterators)->get_dof_indices > (data.vel_local_dof

Re: [deal.II] step-35 in modified

5 March 2022 at 01:02:42 UTC+2 Wolfgang Bangerth wrote: > On 3/4/22 14:37, jack urombo wrote: > > ** > > > > Does this replace all the four instances of *std::get *or just the > > last two? > > All of the

Re: [deal.II] step-35 in modified

I need suggestions on how to properly define or indications of where I getting wrong. The code is attached. On Saturday, 5 March 2022 at 05:42:26 UTC+2 jack urombo wrote: > Yeah this: > > > > *scratch.fe_val_vel.reinit(std::get<0> (*SI)); > scratch.fe_val_pres.reini

[deal.II] Re: MPI and Docker

2021 at 9:54:42 AM UTC+1 jack urombo wrote: > >> The current docker image is not finding mpi.h >> >> >> >> >> >> >> >> >> >> >> >> *[ 50%] Building CXX object CMakeFiles/cdr.dir/cdr.cc.o In file included >> fro

[deal.II] 3D Periodic BCs

I am trying to set periodic boundary conditions for a 3D problem like *velocity->set_periodic_boundary_condition(0, 1, 2); velocity->set_periodic_boundary_condition(2, 3, 1); pressure->set_periodic_boundary_condition(0, 1, 2); pressure->set_periodic_boundary_condition(2, 3, 1);* The bou

Re: [deal.II] 3D Periodic BCs

e already been closed Aborting! * On Friday, 8 July 2022 at 00:17:17 UTC+2 d.arnd...@gmail.com wrote: > Jack, > > step-45 https://www.dealii.org/current/doxygen/deal.II/step_45.h

[deal.II] Dimension Mismatch

i am using the implimenation by Jaekwang Kim to model polymer flows. The major change that I want to do is to output convergence tables like in Step-7 and Step-22, which means I have to have an ExactSolution. The code is encountering an error because the dimension of the solutions do not match

[deal.II] Re: Dimension Mismatch

or&, const dealii::Quadrature&, const Nor mType&, const dealii::Function*, double)’ The issues seem to be emanating from the form of the *exact_solution *and how the components are selected. Ant suggestions on resolving this? --

[deal.II] Re: Dimension Mismatch

The code is attached this message. On Wednesday, 8 March 2023 at 05:38:43 UTC+2 jack urombo wrote: > Took a brief hiatus. and managed to pass a huddle following Marc's > suggestion. > > Now I have set to construct convergence tables for the polymer flow > problem I am

Re: [deal.II] Re: Dimension Mismatch

n() #8 ./oldroyd_bn: main make[3]: *** [CMakeFiles/run.dir/build.make:71: CMakeFiles/run] Aborted (core dumped) make[2]: *** [CMakeFiles/Makefile2:116: CMakeFiles/run.dir/all] Error 2 make[1]: *** [CMakeFiles/Makefile2:123: CMakeFiles/ru

[deal.II] Example not Compiling in deallii on Docker

I have installed dealii on docker and when itry t run the examples i get the error below. What am i doing wrong? /usr/local/opt/spack/linux-ubuntu16.04-x86_64/gcc-5.4.0/dealii-9.1.1-kbq6c5p67nir5zwpx5lbevwutndfivxz/share/deal.II/examples/step-6/CMakeFiles/CMakeTmp/testCCompiler.c try_c

[deal.II] Compiling release 9.5.0

The new release has this error on compiling. [ 10%] Building CXX object bundled/tbb-2018_U2/src/CMakeFiles/bundled_tbb_debug.dir/tbb/ar ena.cpp.o In file included from /home/ju.rom.bo/downloads/dealii-9.5.0/bundled/tbb-2018_U2/src/tbb/s cheduler_common.h:43, from /home/ju.rom.

Re: [deal.II] Compiling release 9.5.0

Here is the detailed log. On Sat, 8 Jul 2023 at 17:37, Matthias Maier wrote: > Hi, > > Would you mind to also send us the "detailed.log" file located in your > build directory? > > Best, > Matthias > > > > On Sat, Jul 8, 2023, at 03:16 CDT, jack ur

[deal.II] How to use a numerical boundary condition in deal.II

Hello: I want to know how to use a numerical boundary condition in deal.II. For example, I want to solve a set of diffusion equations: \partial c_1 / \partial t = d_1 \partial c_1^2 / \partial x^2 \partial c_2 / \partial t = d_2 \partial c_2^2 / \partial x^2 on the one dimensional domain [0, L]. A

[deal.II] Received signal SIGSEGV, Segmentation fault

Hello, I am trying to compile DEAL.II on a server. There is no error in the compilation, but when I try to run an example, segmentation fault occurs. I create a environment by conda, in this environment, I installed gcc, gxx， gfortan, make, cmake, lapack, blas, scalapack, petsc, symengine I also

[deal.II] Minimal deal.ii application loads many surprising dynamic libraries on macos

stand the issue well enough to know if there is anything I or the deal.ii developers can do about it! In any case, 0.25s is a very small price to pay for so much functionality and expressiveness. Thank you, -Jack -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum

Re: [deal.II] Minimal deal.ii application loads many surprising dynamic libraries on macos

p in about half the time (also just built against MPI). -Jack On Friday, April 19, 2024 at 2:00:29 PM UTC-7 Wolfgang Bangerth wrote: On 4/18/24 23:29, Jack Coughlin wrote: > > As I say, I don't understand the issue well enough to know if there is > anything I or the deal.ii develo

Re: [deal.II] Matrix-free: two-point volume fluxes for entropy-stable DG?

Martin, I haven't tried implementing it with FEValues yet, but when I do I will loop you in. -Jack On Wednesday, April 24, 2024 at 1:28:07 PM UTC-7 Martin Kronbichler wrote: > Dear Jack, > > Two-point fluxes are something we have not yet included in deal.II. I > would be v

Re: [deal.II] Matrix-free: two-point volume fluxes for entropy-stable DG?

g the quadrature rule and FiniteElement. Best, -Jack On Wednesday, April 24, 2024 at 1:28:07 PM UTC-7 Martin Kronbichler wrote: > Dear Jack, > > Two-point fluxes are something we have not yet included in deal.II. I > would be very interested if we could realize this feature. I mu

[deal.II] adaptive mesh refinement for Navier Stokes

I am working on a adaptive finite elements for flows. For starters I want to base my simulation on Step 35 and introduce the refine_mesh( )  like in Step 26 to it: See attached code snippet. When I make I get the error message: /[ 50%] Building CXX object CMakeFiles/step-35a.dir/step-35a.cc.o

[deal.II] step-35 in modified

Hie. I am working on flow problems with a model based on the NSE.  The project involves using adaptive finite elementts, and first steps are  to: 1) have a mesh adapting strategy 2) compute errors 3) obtain convergence rates I am basing these on the dealii step-35

Re: [deal.II] step-35 in modified

|   class NavierStokesProjection// //  | ^~// //make[3]: *** [CMakeFiles/step-35a.dir/build.make:63: CMakeFiles/step-35a.dir/step-35a.cc.o] Error 1// / On 11/07/2021 21:24, Wolfgang Bangerth wrote: On 7/10/21 4:23 AM, andrew jack urombo wrote: I have attached

[deal.II] Docker and Eclipse

I am running dealii on Docker, and now would like to debug my projects using Eclipse. 1. How do you that so that Eclipse recognise the dealii that is in a container and work in sync? Why I am not running a compiled dealii is because I am getting//this error with log ending: //usr/lib64/gcc