Re: [deal.II] Nested reverse and forward-mode Sacado memory leak

2020-08-11 Thread Doug Shi-Dong
Thanks for the answer! I am currently using the master branch of Trilinos, to make sure that a newer version did not fix this. I had some second-guessing with how I had configured Sacado itself since there seems to be some options that changes how the reverse mode reuses the memory through

[deal.II] Nested reverse and forward-mode Sacado memory leak

2020-08-10 Thread Doug Shi-Dong
Hello everyone, I found that Sacado's nested reverse and forward mode AD leads to a memory leak when use multiple times. I have opened an issue on their end: https://github.com/trilinos/Trilinos/issues/7741 and the steps to reproduce are pretty simple. However, I've found that the community

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Doug Shi-Dong
Yes I can confirm, I am using QProjector with a dummy face Quadrature right now it works as intended. Doug -- 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

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Doug Shi-Dong
es Point as input. >> >> The dimensions of points in cells is different from that on faces. That’s >> one of the reasons why I posted this question. >> >> >> If I misunderstand you, please kindly let me know. >> >> >> >> Regards, >&

[deal.II] Re: How to transform points on faces from unit cells to real cells?

2020-06-21 Thread Doug Shi-Dong
sing the member function > "transform_unit_to_real_cell()" to get all support points in a unit cell? > In the step, sort out the face points I need? Am I correct? > > > Regards, > > Lex > > On Sunday, June 21, 2020 at 5:00:35 PM UTC-7, Doug Shi-Dong wrote: >> >> Hello L

[deal.II] Re: How to transform points on faces from unit cells to real cells?

2020-06-21 Thread Doug Shi-Dong
Hello Lex, You were close. https://www.dealii.org/current/doxygen/deal.II/classMapping.html#ae5df63553eb8ed170c3b90524853dd48 The project_real_point_to_unit_point_on_face() is useful is you take a "volume" point within the cell. Since you are mapping from unit to real, you would always know

[deal.II] Re: Strategy to snap the boundary of a triangulation to a manifold

2020-06-09 Thread Doug Shi-Dong
Prof. Blais, Bonjour fellow Montreal-er (McGill), land of the bagel and poutine. I have also given some thoughts to robustly generate our meshes with GMSH. I think it's great to see that we've concluded to similar steps as you described: 1. Provide GMSH some CAD/NURBS parametrization from

[deal.II] Re: LinearOperator MPI Parallel Vector

2020-04-22 Thread Doug Shi-Dong
More so, I just found that the AffineConstraints function have not been instantiated for a mix of TrilinosWrappers::SparseMatrix and LA::dist::Vector, hence it likely has not been tried out/tested. Seems like it's just not a thing to use LA::dist::Vector other than for matrix-free? -- The

[deal.II] LinearOperator MPI Parallel Vector

2020-04-22 Thread Doug Shi-Dong
Hello, I have been using the LinearAlgebra::distributed::Vector class for MPI parallelization since the way it works is more familiar to what I had worked with and seemed more flexible. However, for parallelization, I have to either use a Trilinos or PETSc matrix since the native deal.II

Re: [deal.II] Constrained linear operator on distributed data structures

2019-12-31 Thread Doug Shi-Dong
i Doug, > > If I recall correctly, we added that comment since this is a part of > the library that is not very well tested. The code is probably fine. > > Thanks, > David > > On Tue, Dec 24, 2019 at 2:33 PM Doug Shi-Dong > wrote: > > > > Hello everyon

[deal.II] Re: Converting LinearAlgebra::distributed::Vector to TrilinosWrappers::MPI::Vector

2019-12-24 Thread Doug Shi-Dong
, the resulting Trilinos vector was empty and would of course not have the same parallel distribution as the matrix it was supposed to be multiplied with. Doug On Wednesday, December 4, 2019 at 6:25:49 PM UTC-5, Doug Shi-Dong wrote: > > Hello, > > I am trying to solve a block

[deal.II] Constrained linear operator on distributed data structures

2019-12-24 Thread Doug Shi-Dong
Hello everyone, I see that some functions within the constrained_linear_operator.h have a note saying Currently, this function may not work correctly for distributed data structures. I am currently using it in parallel and seems to be working fine for a very simple linear elasticity problem

[deal.II] Converting LinearAlgebra::distributed::Vector to TrilinosWrappers::MPI::Vector

2019-12-04 Thread Doug Shi-Dong
Hello, I am trying to solve a block system using the LinearOperators and Trilinos matrices and vectors. Currently, all my matrices uses TrilinosWrappers::SparseMatrix, but my vectors use the deal.II LinearAlgebra::distributed::Vector because of how easily it handles ghosted vectors. When

Re: [deal.II] Instantiation of vector_tools_integrate_difference

2019-11-06 Thread Doug Shi-Dong
Not sure why it builds fine locally, but that fixed it! Thank you! Doug On Wednesday, November 6, 2019 at 11:06:47 AM UTC-5, Wolfgang Bangerth wrote: > > On 11/6/19 8:49 AM, Doug Shi-Dong wrote: > > > > Any hints about what could be causing this? > > I thou

[deal.II] Re: Docker MPI

2019-10-23 Thread Doug Shi-Dong
stalled. > > Best, > > Bruno > > On Monday, October 21, 2019 at 9:53:36 PM UTC-4, Doug Shi-Dong wrote: >> >> Hello everyone, >> >> I would like to use Travis CI to test pull requests, and I am currently >> following the steps described in &g

[deal.II] Docker MPI

2019-10-21 Thread Doug Shi-Dong
Hello everyone, I would like to use Travis CI to test pull requests, and I am currently following the steps described in https://github.com/dealii/dealii/wiki/Docker-Images Except I am using the following Docker image: dealii/dealii:v9.1.1-gcc-mpi-fulldepsspack-debugrelease It seems that

[deal.II] Re: Query regarding DoFTools::dof_indices_with_subdomain_association()

2019-10-10 Thread Doug Shi-Dong
Hello Vachan, Sounds like you're implementing nodal DG, hence why you only need values and gradients at the *quadrature* points from the neighbor. Something you might want to consider rather than communicating them at solution points, in case you ever decide to overintegrate. You could still

[deal.II] Re: Integrate the solution over a subdomain

2019-10-08 Thread Doug Shi-Dong
On Tuesday, October 8, 2019 at 8:43:23 AM UTC-4, Yang Liu wrote: > > > Dear DEAL.II users, > > The DEAL.II software provides the function > VectorTools::compute_mean_value() to integrate the solution over the whole > domain. > > But now I would like to integrate the solution over a subdomain.

[deal.II] Re: High order mesh from Gmsh

2019-09-20 Thread Doug Shi-Dong
Hello everyone, Sorry to revive/hijack this thread, but it seems to have the exact information/people that might be able to help. My current goal is to perform some shape optimization with some h(p)-refinement. I am currently using Prof. Heltai's method described in this issue

Re: [deal.II] Instantiation MappingFEField with hp::DoFHandler

2019-09-18 Thread Doug Shi-Dong
On Wednesday, September 18, 2019 at 11:28:45 PM UTC-4, Wolfgang Bangerth wrote: > > > Doug, > > It does not greatly surprise me that it doesn't compile -- what isn't > tested > generally doesn't work. > Makes sense. We'd of course be very happy to accept any patches to make things work.