[deal.II] deal.II cmake stuck at "Performing Test DEAL_II_HAVE_USABLE_FLAGS_DEBUG"

Dear all I am trying to install on a new machine running opensuse tumbleweed. I installed packages using spack $ spack find *==>* 37 installed packages. -- linux-opensuse20161212-x86_64 / gcc@6 arpack-ng@3.4.0 m4@1.4.17 pkg-config@0.29.1

Re: [deal.II] build_patches(degree) for higher order dg elements in parallel

Justin, I've attached a small test case that simply initializes a parallel vector and then attempts to write it to file via DataOutFaces. The finite element used only lives on the faces of the cells so DataOutFaces in appropriate. This test case provided me with some insight into my problem.

Re: [deal.II] build_patches(degree) for higher order dg elements in parallel

Wolfgang, I am using deal.ii 8.5.0-pre so the patch is in my version. I forgot to change my value of reputations for my triangulation to take the vector arguments. The updated code produces the error on the coarse mesh, but works as desired on the finer mesh. Since I am still new to parallel

[deal.II] make test

I have compiled dealii with the following; cmake -DDEAL_II_WITH_MPI=ON \ -DDEAL_II_WITH_TRILINOS=ON\ -DTRILINOS_DIR=~/trilinos/ \ -DCMAKE_INSTALL_PREFIX=~/dealII_1 \ ~/dealii-8.4.1 I installed trilinos12.4.2. I could compile and install dealii with the above command.

[deal.II] obtaining solution values at any point using Quadrature< dim >::Quadrature(const Point< dim > & point)

Dear all, I was trying to write a roots finding method along cell faces in 2D, so I need to know the solution values at any point of cell faces. Is this the right way to use Quadrature < dim >:: Quadrature

Re: [deal.II] obtaining solution values at any point using Quadrature< dim >::Quadrature(const Point< dim > & point)

On 12/15/2016 06:11 PM, Jiaqi ZHANG wrote: std::cout<<"testing quadrature points:\n"; const Point v0 = cell->vertex(0); Quadrature q_point(v0); Quadrature objects are defined on the reference cell [0,1]^d. So you will need to pass a point in reference coordinates to the quadrature object, but

Re: [deal.II] obtaining solution values at any point using Quadrature< dim >::Quadrature(const Point< dim > & point)

Got it. Thanks so much, Dr. Bangerth. Jiaqi 2016-12-16 0:06 GMT-05:00 Wolfgang Bangerth : > On 12/15/2016 06:11 PM, Jiaqi ZHANG wrote: > >> std::cout<<"testing quadrature points:\n"; >> const Point v0 = cell->vertex(0); >> Quadrature q_point(v0); >> > > Quadrature

Re: [deal.II] build_patches(degree) for higher order dg elements in parallel

Justin, I am using deal.ii 8.5.0-pre so the patch is in my version. I forgot to change my value of reputations for my triangulation to take the vector arguments. The updated code produces the error on the coarse mesh, but works as desired on the finer mesh. Since I am still new to parallel

[deal.II] Re: deal.II cmake stuck at "Performing Test DEAL_II_HAVE_USABLE_FLAGS_DEBUG"

Additionally I added Fortran compiler and fixed "suitesparse" to "suite-sparse", but I am still stuck at the same place as before. Should I use "Filesystem View" to compile my own dealii with spack packages ? If I install dealii@develop using spack, how do I update dealii in future ? Thanks

[deal.II] Re: deal.II cmake stuck at "Performing Test DEAL_II_HAVE_USABLE_FLAGS_DEBUG"

Dear Praveen, Unfortunately I don't have any suggestions about CMake hanging. I've just reinstalled deal.II through Spack on both my Apple and Linux machines and I didn't run into the problem you're experiencing. Denis has some experience with installing deal.II though spack on a number of

[deal.II] Re: Tensor output via ConditionalOStream pcout

Dear S.A. Mohseni, There is no need to duplicate code. ;-) FullMatrix also has a print function[1] and its implementation[2] is more or less identical to what you are doing. Best, Daniel [1] https://www.dealii.org/8.4.1/doxygen/deal.II/classFullMatrix.html#a6427bd9b23b03c590b9881c7dc551552

[deal.II] Re: Tensor output via ConditionalOStream pcout

I see. Then I will output the local matrices and check the overall structure accordingly. I wrote a simple function to output the cell matrix. Maybe someone as noob as me, may find it useful :) template void SolidMechanics::print_matrix (FullMatrix& matrix, char* name) { std::cout << std::endl

[deal.II] Re: Tensor output via ConditionalOStream pcout

I see. Then I will output the local matrices and check the overall structure accordingly. I wrote a simple function to output the cell matrix. Maybe someone as noob as me, may find it useful :) template void SolidMechanics::print_matrix (FullMatrix& matrix, const unsigned int dofs_per_cell ,

[deal.II] Re: Tensor output via ConditionalOStream pcout

Dear Seyed, Thanks for sharing your print function. For future reference, the FullMatrix class has two print

Re: [deal.II] Re: ifem

Dear Joaquín, >From your previous posts I assume that all of PeFT, Pe[qs] and F[qs] represent rank 2 tensors. If this is correct then no, the contraction indices that you've specified are not correct. They're out of the range [0,1] that are valid for rank-2 tensors, and this would presumably