Re: [deal.II] Re: Coarse Cell IDs Distributed Triangulation

2021-11-11 Thread Jonathan Russ
Amazing. Thank you very much! Jonathan On Thursday, November 11, 2021 at 4:33:47 PM UTC-5 Wolfgang Bangerth wrote: > On 11/11/21 1:08 PM, Jonathan Russ wrote: > > Does the same cell in the bottom left corner of the domain have exactly > > the same CellId as the CellId it had wh

[deal.II] Re: Coarse Cell IDs Distributed Triangulation

2021-11-11 Thread Jonathan Russ
rg/current/doxygen/deal.II/classCellAccessor.html#aee7d5ec9cac8af95a136554501ea4f0a> > . > > Hope this helps, > Marc > > On Thursday, November 11, 2021 at 6:30:30 AM UTC-7 Jonathan Russ wrote: > >> Hello - >> >> I really would like to have a unique

[deal.II] Coarse Cell IDs Distributed Triangulation

2021-11-11 Thread Jonathan Russ
Hello - I really would like to have a unique ID for each coarse cell regardless of the number of MPI processes. I have a distributed (not fullydistributed) triangulation. Testing the coarse_cell->index

Re: [deal.II] Extract constant modes slow

2020-01-24 Thread Jonathan Russ
for taking your time. Jonathan On Friday, January 24, 2020 at 12:49:34 PM UTC-5, Wolfgang Bangerth wrote: > > On 1/23/20 5:00 PM, Jonathan Russ wrote: > > > > I noticed that the Trilinos ML preconditioner accessed through the > > TrilinosWrappers benefits heavily from supp

[deal.II] Extract constant modes slow

2020-01-23 Thread Jonathan Russ
Hello - I noticed that the Trilinos ML preconditioner accessed through the TrilinosWrappers benefits heavily from supplying the constant modes or the null-space of the laplace operator. Using the extract_constant_modes functionality I notice that this can take quite a bit of time for a vector

Re: [deal.II] Iterative Linear Solver Tolerance

2019-09-14 Thread Jonathan Russ
Professor Bangerth - Thank you very much for your very helpful reply! I was thinking that the addition of large vector of numbers at or below machine precision would accumulate a lot of numerical roundoff error but I see what you mean. The Eisenstat Walker algorithm looks like exactly what I

[deal.II] Iterative Linear Solver Tolerance

2019-09-13 Thread Jonathan Russ
Hello - I have a simple question I was hoping to get advice on related to setting the linear solver tolerance for an iterative solver. In the tutorial examples (and in most of the codes I have seen and written) the solver tolerance is some fraction of the initial right hand side 2-norm (e.g.

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

Re: [deal.II] Re: Serendipity elements

2019-08-04 Thread Jonathan Russ
uires a significant amount of additional effort? Thank you again, Jonathan On Sunday, August 4, 2019 at 3:58:07 PM UTC-6, Wolfgang Bangerth wrote: > > On 8/3/19 1:17 PM, Jonathan Russ wrote: > > > > How difficult would it be to add serendipity elements to dealii? They >

[deal.II] Re: Serendipity elements

2019-08-03 Thread Jonathan Russ
Hello - How difficult would it be to add serendipity elements to dealii? They are very useful in solid mechanics applications and can greatly reduce the computational cost. Is it very difficult to add an element that is the same is FE_Q but without the few extra basis functions when quadratic

[deal.II] Re: Output global variable such as "time"

2019-04-20 Thread Jonathan Russ
Perfect! Thank you so much! On Saturday, April 20, 2019 at 9:31:51 PM UTC-4, Daniel Arndt wrote: > > Jonathan, > > have a look at DataOutBase::write_pvd_record() ( > https://www.dealii.org/developer/doxygen/deal.II/namespaceDataOutBase.html#a6f1c052ba49fd44cd8e3f35ba871aebd > ). > > Best, >

[deal.II] Output global variable such as "time"

2019-04-20 Thread Jonathan Russ
Hello - Is it possible to output a "global variable" such as time to, say, the paraview output files so that I can easily see at what solution time the current output data corresponds to? I only see the ability to output vectors of data that correspond to either the number of nodes or elements

[deal.II] Re: Saving distributed vector and FullMatrix to binary files

2018-10-23 Thread Jonathan Russ
Professor Bangerth - Thank you for your reply. I may try this but I also plan to look into HDF5 to see if that might be effective. Thank you again! Jonathan On Friday, October 19, 2018 at 1:20:37 PM UTC-4, Jonathan Russ wrote: > > Good afternoon - > > I am implementing an algorith

[deal.II] Saving distributed vector and FullMatrix to binary files

2018-10-19 Thread Jonathan Russ
Good afternoon - I am implementing an algorithm in which I need to save some data and the distributed solution vector at each time increment for later use. Since this is potentially a lot of data that I cannot store in available memory I would like to be able to write this data out to binary

Re: [deal.II] FE w/ Gauss Point DoFs

2018-03-07 Thread Jonathan Russ
gt; mass matrix can be made diagonal by choosing basis functions that are not > either zero or one at the vertices of the cell, but instead are zero or one > at a given set of quadrature points. > > > J-P > > On 07 Mar 2018, at 18:09, Jonathan Russ <jr3...@columbia.edu >

Re: [deal.II] Re: FE w/ Gauss Point DoFs

2018-03-07 Thread Jonathan Russ
, 2018 at 11:04:33 AM UTC-5, Bruno Turcksin wrote: > > Jonathan, > > 2018-03-07 10:43 GMT-05:00 Jonathan Russ <jr3...@columbia.edu > >: > >> >> Thank you for your fast reply. I'm sorry I'm not doing a great job of >> explaining this. Here is an example from el

[deal.II] Re: FE w/ Gauss Point DoFs

2018-03-07 Thread Jonathan Russ
of freedom. On Wednesday, March 7, 2018 at 10:46:12 AM UTC-5, Jonathan Russ wrote: > > Some people refer to this as an element with internal nodes. > > On Wednesday, March 7, 2018 at 10:43:13 AM UTC-5, Jonathan Russ wrote: >> >> Hi Bruno - >> >> Thank yo

[deal.II] Re: FE w/ Gauss Point DoFs

2018-03-07 Thread Jonathan Russ
Some people refer to this as an element with internal nodes. On Wednesday, March 7, 2018 at 10:43:13 AM UTC-5, Jonathan Russ wrote: > > Hi Bruno - > > Thank you for your fast reply. I'm sorry I'm not doing a great job of > explaining this. Here is an example from elasticity that i

[deal.II] FE w/ Gauss Point DoFs

2018-03-07 Thread Jonathan Russ
Hello - I apologize since this question is likely obvious to most, but a quick search didn't reveal the answer to me. Is there a FiniteElement in deal.II that allows the addition of degrees of freedom at gauss points so that the DoFHandler object also manages them? For example, I would like

[deal.II] Re: Output incorrect along parallel partition

2016-06-07 Thread Jonathan Russ
Hi Artur - Unfortunately I don't have any hanging node constraints in the model. I do have constraints but there aren't any along the mesh partition. I still can't figure out why my output looks the way it does. The solve function looks like this: unsigned int EigenvalueProblem::solve () {