Re: [deal.II] Alternative to local_range in LinearAlgebra::distributed::Vector

> Thanks. If I dont do any reordering then, continuous nature holds, right ? > Does this hold even after a grid refinement ? Yes. If you call `dof_handler.distribute_dofs()`, then the index range is contiguous. > I can create this info myself > > // First global index on this partition >  

Re: [deal.II] Alternative to local_range in LinearAlgebra::distributed::Vector

> On 02-Feb-2019, at 8:41 PM, Wolfgang Bangerth wrote: > > No :-) The reason why it was scheduled for removal is that it only works in > case the locally owned range is really just a single interval. But that's not > always the case. For example, you may have reordered degrees of freedom in

[deal.II] Re: Role of CUDAWrappers::MatrixFree::Data and CUDAWrappers::SharedData in CUDAWrappers::MatrixFree:cell_loop

Steve, On Wednesday, February 6, 2019 at 2:49:26 PM UTC-5, Stephen DeWitt wrote: > The one area that I don't really understand for > matrix_free_matrix_vector_01.cu is how information flows into the > HelmholtzOperator >

[deal.II] Role of CUDAWrappers::MatrixFree::Data and CUDAWrappers::SharedData in CUDAWrappers::MatrixFree:cell_loop

Hello all, Following up from a previous post about where the GPU capabilities are, I've been going through the CUDA matrix-free matrix-vector product tests as Bruno suggested. For the most part

[deal.II] Re: command line compilation error on a linux machine using spack installation of deal.ii dev

> > Dear all, > I have found the solution thanking to Alberto Sartori who suggested that the problem could be a repeated definition between my main.cc and something else somewhere in deal.ii. In my CMakeLists.txt I defined DIMENSION, defining IMENSION the compilation works fine. So we ju

[deal.II] Accuracy of Dirichlet condition for p in step-20

Hi all, I am trying to solve a system of equations that do this: Stokes to solve for v_r and p_r for one fluid (viscous rock), I use these solutions on the RHS of a Darcy type equation solved like step-20 for the pressure p_f in the fluid in the domain. Using the 3 solutions, I update another c

Re: [deal.II] Automatic Differentiation of the Finite Element Jacobian

Doug, in addition to what J-P already wrote: > We had chosen to use deal.ii since most of it looked templated in terms of > accepting the float type. However, I noticed that the FEValuesBase returns a > 'double', instead of a generic 'value_type' for the Jacobian term. That depends on the con

Re: [deal.II] Alternative to local_range in LinearAlgebra::distributed::Vector

> On 2/2/19 12:32 AM, Praveen C wrote: >> Thank you for this tip. For my need, I only require the starting global >> index >> on each partition since I am using local_element to do some operations like >> multiplying by a local time step which is different for each cell. The >> local_range fu