[deal.II] block_write and block_read in dealii::LinearAlgebraTrilinos::MPI::BlockVector

2017-04-05 Thread Lam DANG
Dear All I have a success of *block_write* and *block_read* in storing *Vector::BlockVector* for post processing *https://www.dealii.org/8.4.0/doxygen/deal.II/classVector.html#aef98e4994a6c0b20bb33b0290f8d607c* Now, I would do similar in *dealii::LinearAlgebraTrilinos::MPI::BlockVector* but

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

2017-04-05 Thread Alex Zimmerman
Thanks for the clarification, Jean-Paul. I'm glad to hear that my syntax is correct. I am indeed using an extractor, e.g. here's a line from my code: const Tensor<2, dim> gradv = fe_values[velocity].gradient(j, quad); I'm having some trouble finding an example of the S[i][j] syntax in either

[deal.II] Re: run deal.ii on stampede

2017-04-05 Thread zhen-lin Wang
Hi Artur, I am trying to install these by myself, it would be great that you can share your installation steps, did you configure deal.ii with petsc? If so did you use the available version or also install it by yourself? Thank you, Zhenlin 在 2017年4月5日星期三 UTC-4下午7:56:52,Artur Safin写道: > >

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

2017-04-05 Thread Jean-Paul Pelteret
> > Hi Alex, >> > > Thanks for the clarification, Jean-Paul. > I'm sorry to say that you shouldn't be too quick in this instance to thank me... > I'm glad to hear that my syntax is correct. > > I am indeed using an extractor, e.g. here's a line from my code: > const Tensor<2, dim> gradv

Re: [deal.II] block_write and block_read in dealii::LinearAlgebraTrilinos::MPI::BlockVector

2017-04-05 Thread Wolfgang Bangerth
On 04/05/2017 02:13 AM, Lam DANG wrote: // Now, I would do similar in *dealii::LinearAlgebraTrilinos::MPI::BlockVector* but there is no *block_write* and *block_read* in it. /http://dealii.org/8.4.1/doxygen/deal.II/classTrilinosWrappers_1_1MPI_1_1BlockVector.html/ *Could you help me give any

[deal.II] run deal.ii on stampede

2017-04-05 Thread zhen-lin Wang
Hi, I am trying to run deal.ii on Stampede after loading following modules, I can compile the code: 1) xalt/0.6 3) git/2.7.0 5) intel/15.0.2 7) mvapich2/2.1 9) dealii/8.4.1 11) slepc/3.5 2) TACC 4) cmake/3.7.1 6) cxx11/4.9.18) trilinos/12.6.4 10) petsc/3.5

[deal.II] Re: run deal.ii on stampede

2017-04-05 Thread zhen-lin Wang
Hi Bruno, To load trilinos/12.8.1 I have to load gcc which is against intel/15.0.2 that is necessary to load dealii. I load the python module, but the error is still there, probably I will try to install trilinos and deal.ii bymyself. Thanks, Zhenlin 在 2017年4月5日星期三 UTC-4下午7:34:35,Bruno

[deal.II] Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

2017-04-05 Thread Alex Zimmerman
I am debugging a code where part of my weak form requires indexing into the Tensor<2,dim> that results from getting the gradient of a Tensor<1,dim> (in this case, a velocity vector). I found the convention for the index notation here:

[deal.II] Re: Arbitrary mapping or Manifold?

2017-04-05 Thread Juan Carlos Araujo Cabarcas
Dear all, I recently found the thread: Something wrong with ChartManifold when chartdim=1: https://groups.google.com/forum/#!topic/dealii/Sfo9xKoeRpw a more straight answer to this thread. I took the liberty to copy David's code and modify it, so that the upper face of a square follows:

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

2017-04-05 Thread Jean-Paul Pelteret
Dear Alex, What you've stated is correct. The tensor component $S_{ij}$ is exactly component S[i][j] of a Tensor<2,dim> S if "S" represents $S_{ij} = \partial v_{i} \partial x_{j}$. As for an example, both steps 18

[deal.II] Re: run deal.ii on stampede

2017-04-05 Thread Bruno Turcksin
Zhenlin, On Wednesday, April 5, 2017 at 6:10:44 PM UTC-4, zhen-lin Wang wrote: > > However, errors happens at linking: > > /home1/apps/intel15/mvapich2_2_1/trilinos/12.8.1/lib/libpytrilinos.so: > undefined reference to `PyType_IsSubtype' > >

[deal.II] Re: run deal.ii on stampede

2017-04-05 Thread Artur Safin
Zhenlin, I remember having trouble getting deal2 to compile against the available versions of trilinos. In the end, I just ended up installing both myself. If you want, I can post my installation steps. Rgds, Artur -- The deal.II project is located at http://www.dealii.org/ For mailing