Re: [deal.II] Re: dealii compilation on Cray machine

Hi Daniel, Ok I do that, now I get the original compilation error for which I started this email thread. *Please verify that both the operating system and the processor support Intel(R) AVX512DQ, AVX512F, AVX512CD, AVX512BW and AVX512VL instructions.* *source/numerics/CMakeFiles/obj_numerics_ins

Re: [deal.II] Re: dealii compilation on Cray machine

erify that it still works :-D > > Best, > Matthias > > > > On Fri, Apr 5, 2019, at 14:45 CDT, Bruno Turcksin < > bruno.turck...@gmail.com> wrote: > > > Le ven. 5 avr. 2019 à 14:48, Phani Motamarri a écrit > : > >> > >> No login node is in

Re: [deal.II] Re: dealii compilation on Cray machine

pr 6, 2019 at 3:45 PM Martin Kronbichler < kronbichler.mar...@gmail.com> wrote: > > On 05.04.19 21:45, Bruno Turcksin wrote: > > Le ven. 5 avr. 2019 à 14:48, Phani Motamarri a écrit > : > >> No login node is intel-broadwell and compute node is intel-skylake. > Prob

Re: [deal.II] Re: dealii compilation on Cray machine

iday, April 5, 2019 at 1:21:45 PM UTC-4, Phani Motamarri wrote: >> >> Hi Daniel, >> >> Many thanks for the reply. As you suggested, I tried a minimal >> configuration just using >> *cmake -DDEAL_II_ALLOW_AUTODETECTION=OFF >> -DDEAL_II_ALLOW_

Re: [deal.II] Re: dealii compilation on Cray machine

Hi Daniel, Many thanks for the reply. As you suggested, I tried a minimal configuration just using *cmake -DDEAL_II_ALLOW_AUTODETECTION=OFF -DDEAL_II_ALLOW_PLATFORM_INTROSPECTION=OFF -DCMAKE_C_COMPILER=cc -DCMAKE_CXX_COMPILER=CC -DCMAKE_Fortran_COMPILER=ftn -DMPI_C_COMPILER=cc -DMPI_CXX_COMPILER=C

Re: [deal.II] Re: dealii compilation on Cray machine

Thank you Daniel for the reply. It was my bad of pointing to the wrong scalapack path. Now the configure step passes. When I start doing "make" , I get the following error after finishing 35% of the build process: *Scanning dependencies of target obj_sundials_release* *[ 34%] Building CXX object

Re: [deal.II] numbering of degrees of freedom

Thank you Prof Wolfgang for the reply. Yes I meant DoF indices. Since get the level 1 mesh globally refined from level 0 (hyper_cube), if I make both of these levels active at the same time, does dealii default numbering scheme gets the dofs on the lowest level numbered first followed by the number

[deal.II] numbering of degrees of freedom

Hi, Regarding the distribution of indices for degrees of freedom it says in the manual " Since the triangulation is traversed starting with the cells of the coarsest active level and going to more refined levels, the lowest numbers for dofs are given to the largest cells as well as their bounding

[deal.II] Debug error in Kelly Error Estimate function using parallel distributed vectors

Hi, I have just written a plan vanilla interface to call the Kelly Error Estimate function passing a parallel distributed vector as input vector type. But strangely enough when I run on MPI on multiple procs in debug mode, I get the following error: ---

Re: [deal.II] Sparse ILU

i On Mon, Apr 2, 2018 at 2:28 PM, Wolfgang Bangerth wrote: > On 04/02/2018 11:48 AM, Phani Motamarri wrote: > >> We have an ongoing research project with Caltech where we are trying to >> explore ways of approximating the subspace spanned by k lowest eigen values >> of a

Re: [deal.II] Sparse ILU

We have an ongoing research project with Caltech where we are trying to explore ways of approximating the subspace spanned by k lowest eigen values of a given positive definite matrix(FEM discretized matrix) in a sparse format using the incomplete Cholesky factor with reduced computational complexi

[deal.II] Sparse ILU

Dear All, If I have a SparseMatrix and compute SparseILU restricting the sparsity to the given matrix using the functions/options given in https://www.dealii.org/developer/doxygen/deal.II/classSparseILU.html, I was wondering if there is anyway to access the Sparse matrix obtained after the SparseI

Re: [deal.II] p4est critical number of elements in each Proc

PI bug? Could you share your solution? I have a similar > deadlock. > > Thank you! > > Best, > > Qing > > On Tuesday, December 12, 2017 at 9:29:36 AM UTC-7, Phani Motamarri wrote: >> >> Thanks for your reply Dr Wolfgang. It was my bad, I had a MPI bug in my >>

[deal.II] Computing Overlap Matrix

Dear All, I have a STL array of parallel distributed vectors *std::vector X;* In the above VectorType denotes parallel::distributed::Vector. The size of X (numVectorsX) can be close to 5,000 to 10,000 for very large problems. All these vectors in X form a subspace and I would like to evalu

Re: [deal.II] p4est critical number of elements in each Proc

Thanks for your reply Dr Wolfgang. It was my bad, I had a MPI bug in my code which was causing this problem. Now everything works fine! On Mon, Dec 11, 2017 at 12:32 AM, Wolfgang Bangerth wrote: > On 12/10/2017 09:39 PM, Phani Motamarri wrote: > >> Thank you Dr WolfGang for your rep

Re: [deal.II] p4est critical number of elements in each Proc

PM, Phani Motamarri wrote: > >> >> I was just wondering in a MPI run, is there a critical number of >> finite-elements to be present in the base mesh on each processor when doing >> an adaptive refinement of the base mesh using p4est in deal ii. I had a >> base un

[deal.II] p4est critical number of elements in each Proc

Dear All, I was just wondering in a MPI run, is there a critical number of finite-elements to be present in the base mesh on each processor when doing an adaptive refinement of the base mesh using p4est in deal ii. I had a base uniform mesh of 64 elements and when I try to do adaptive refinement u

Re: [deal.II] Re: Handling Dirichlet Boundary conditions with Hanging Nodes

Thanks Daniel for pointing me to this! On Sat, Dec 9, 2017 at 9:05 AM, Daniel Arndt < daniel.ar...@iwr.uni-heidelberg.de> wrote: > Phani, > > >> I have a general query. Is there anything I should be particularly >> careful when I have to handle non-homogenous dirichlet boundary on a mesh >> with

Re: [deal.II] Re: Seg fault using Extracting hanging node dofs in MPI

Oh thank you very much for the pull request. For now, I was extracting DoFs corresponding to hanging nodes by using is_constrained(dof).in parallel. On Sat, Dec 9, 2017 at 6:26 AM, Daniel Arndt < daniel.ar...@iwr.uni-heidelberg.de> wrote: > Phani, > > apparently no one used this function for a

[deal.II] Seg fault using Extracting hanging node dofs in MPI

Dear All, I am trying to create a vector of type bool to identify which dofs correspond to hanging nodes. I use the function DoFTools::extract_hanging_node_dofs(dofHandler,selectedDofsHanging). As mentioned in the manual, I set the size of the vector selectedDofsHanging to dofHandler.n_dofs().

[deal.II] Handling Dirichlet Boundary conditions with Hanging Nodes

Hi, I have a general query. Is there anything I should be particularly careful when I have to handle non-homogenous dirichlet boundary on a mesh with hanging nodes on the boundary? I am working in the matrix-free data structure format. I am trying to impose non-homogenous boundary conditions on th

Re: [deal.II] Creating periodic Boundary Conditions after doing a uniform subdivision

Sure, I will definitely try to make a patch to improve the error message. So if I understand, I have my own git local branch of deal ii develop branch and I modify the code and commit it in my local branch. After that how do I get it to central system for your review? Thanks Phani On Fri, Nov 10,

Re: [deal.II] Creating periodic Boundary Conditions after doing a uniform subdivision

_iterator>>> &) [with dim = 3, spacedim = 3] The violated condition was: this->n_levels() == 1 Additional information: The triangulation is refined! -- Thanks Phani On Fri, Nov 10, 2017 at 9:07 AM, Wolfgang Bangerth wrote: > O

Re: [deal.II] Creating periodic Boundary Conditions after doing a uniform subdivision

icity > condition is retained (and propagated consistently to the refined levels) > when performing subsequent refinement. > > I hope that this helps. > > Kind regards, > Jean-Paul > > > On 10 Nov 2017, at 09:31, Phani Motamarri wrote: > > Hello all, > >

[deal.II] Creating periodic Boundary Conditions after doing a uniform subdivision

Hello all, I have an existing triangulation and I have enforced periodic boundary conditions and created the relevant constraintMatrix using the following five steps just as in the manual: 1. Create a mesh 2. Gather the periodic faces using GridTools::collect_periodic_faces()

[deal.II] Re: Compilation error during the use of CellDataStorage datastructure

Many thanks for all your inputs. Creating a struct encapsulating a double solved all my problems. Further fixing the assert condition as *(std::is_base_of::value || **std::is_same::value)* also solved my problem. Thanks and Regards Phani On Friday, October 27, 2017 at 8:31:53 PM UTC-4, Phani

[deal.II] Compilation error during the use of CellDataStorage datastructure

Hi, I am trying to create a quadrature Point data in the form of CellDataStorage. My quadrature point data are of of type "double" and hence I do the following CellDataStorage::active_cell_iterator,double> rhoQuadData; typename DoFHandler<3>::active_cell_iterator cell_start = dofHandler.