Re: [deal.II] compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Weixiong Zheng
在 2017年6月6日星期二 UTC-7下午12:07:18,Wolfgang Bangerth写道: > > On 06/06/2017 12:10 PM, Weixiong Zheng wrote: > > It runs in serial. The error occurs when using multiple processors. > > To answer if it really happens when I call compress (), I added a bit of tracking: for (unsigned int k=0; kid() << "

Re: [deal.II] integrating a scalar in a vertical column (e.g. finding hydrostatic pressure at a point)

2017-06-06 Thread Wolfgang Bangerth
Sean, Thanks for the wonderfully clear answer. For my function f, I can basically treat the line integral like a "normal" volume integral, setting the direction using the delta functions. This would be the right approach for a distributed adaptive setting as well, no? Yes. You'd write the i

Re: [deal.II] compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Weixiong Zheng
This is how I setup my sparsity pattern: if (discretization=="DFEM") fe = new FE_DGQ (p_order); else fe = new FE_Q (p_order); dof_handler.distribute_dofs (*fe); locally_owned_dofs = dof_handler.locally_owned_dofs (); DynamicSparsityPattern dsp (locally_owned_dofs); if (discretization==

Re: [deal.II] compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Wolfgang Bangerth
On 06/06/2017 12:10 PM, Weixiong Zheng wrote: It runs in serial. The error occurs when using multiple processors. And it really happens when you call `compress()`? The only thing I can think of is that you didn't set up the sparsity pattern correctly. Can you show us the code that sets up the

Re: [deal.II] compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Weixiong Zheng
It runs in serial. The error occurs when using multiple processors. -- 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

Re: [deal.II] convergence adaptive refinement

2017-06-06 Thread Wolfgang Bangerth
On 06/06/2017 04:43 AM, Lam DANG wrote: Dear All I use /write_ucd/ with /flags.write_faces = true/, but the result of .ucd file is *lack of information of face id 0.* How I can output the information of face id 0? The piece of code: /Triangulation<2> triangulation; GridGenerator::hyp

Re: [deal.II] compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Wolfgang Bangerth
On 06/06/2017 03:13 AM, Weixiong Zheng wrote: Now, I got the error: compressing can go through direction 0 but not 1 with the error at the end of this post. Any ideas how to fix it? Thanks in advance! Weixiong [0]PETSC ERROR: - Error Message -

Re: [deal.II] integrating a scalar in a vertical column (e.g. finding hydrostatic pressure at a point)

2017-06-06 Thread Sean McGovern
Hi Wolfgang, Thanks for the wonderfully clear answer. For my function f, I can basically treat the line integral like a "normal" volume integral, setting the direction using the delta functions. This would be the right approach for a distributed adaptive setting as well, no? It makes no assumpt

Re: [deal.II] Re: How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-06-06 Thread Praveen C
In my applications, I add this to CMakeLists.txt # In release mode, disable some warning messages IF(CMAKE_BUILD_TYPE MATCHES Release) SET( CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wno-tautological-compare" ) SET( CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wno-unused-parameter" ) SET( CMAKE_CXX_FLAGS

Re: [deal.II] Re: How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-06-06 Thread Bruno Turcksin
Juan Carlos On Tuesday, June 6, 2017 at 5:34:09 AM UTC-4, Juan Carlos Araujo Cabarcas wrote: > > The warnings come from my code, which is in developing stage! > In the installation I had before, > [-Wunused-variable],[-Wunused-parameter],[-Wunused-but-set-variable] did > not show up. It's a bu

Re: [deal.II] convergence adaptive refinement

2017-06-06 Thread Lam DANG
Dear All I use *write_ucd* with *flags.write_faces = true*, but the result of .ucd file is *lack of information of face id 0.* How I can output the information of face id 0? The piece of code: *Triangulation<2> triangulation;GridGenerator::hyper_cube (triangulation, 0, 1, true);

Re: [deal.II] Re: How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-06-06 Thread Juan Carlos Araujo Cabarcas
Thanks for answering! The warnings come from my code, which is in developing stage! In the installation I had before, [-Wunused-variable],[-Wunused-parameter],[-Wunused-but-set-variable] did not show up. It's a bummer that changing deal.II version will make such a difference! As I am not profi

[deal.II] Re: compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Weixiong Zheng
在 2017年6月6日星期二 UTC-7上午2:13:28,Weixiong Zheng写道: > > Dear All, > > I am doing a radiation transport problem, which has multiple individual > equations (number determined by the number of directions we input). As the > number of directions is arbitrary, I put the sparsematries are contained in >

[deal.II] compress error for vector of PETScWrappers::MPI::SparseMatrix

2017-06-06 Thread Weixiong Zheng
Dear All, I am doing a radiation transport problem, which has multiple individual equations (number determined by the number of directions we input). As the number of directions is arbitrary, I put the sparsematries are contained in std::vectors as raw pointers (using LA namespace to stand for