[Dune] Periodic Face Transformation Block in DGF Format, what is wrong ?

[Benedikt Oswald]

Dear Dune, when I use PERIODICFACETRANSFORMATION a  in a otherwise perfectly legal DGF File
in order to denote a grid that is:

periodic in x by 1 unit
periodic in y by 2 units
and not periodic in y direction


PERIODICFACETRANSFORMATION
% set periodic boundaries in x,y,z direction with periodicity px= 1.0 and py=2.0
1 0 0, 0 1 0, 0 0 1 + 1.0 0.0 0.0
1 0 0, 0 1 0, 0 0 1 + 0.0 2.0 0.0




then the DGF parser complains with:

Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.
Dune reported error: DGFException [PeriodicFaceTransformationBlock:periodicfacetrans.cc:32]: Error in block PERIODICFACETRANSFORMATION (line 2): Not enough entries in matrix row 0.


What is wrong with my specification ? Indeed, I derived it from the DGF format description for 2D:

PeriodicFaceTransformation 
Each line describes an affine transformation that shall be used to glue grid boundaries together. The transformation is denoted as matrix + shift. The following 2d example describes a shift by the first unit vector:
  1 0, 0 1 + 1 0

I use ALUGrid for simplices.


thanks for having a look! Greetings, Benedikt

PS: dune versions are:

odysseus:hades oswald$ svn info ../dune-common/
Path: ../dune-common
URL: https://svn.dune-project.org/svn/dune-common/trunk
Repository Root: https://svn.dune-project.org/svn/dune-common
Repository UUID: 419231ca-4411-0410-aba2-fc6cf8c4f609
Revision: 6491
Node Kind: directory
Schedule: normal
Last Changed Author: sander
Last Changed Rev: 6491
Last Changed Date: 2011-08-30 10:56:22 +0200 (Tue, 30 Aug 2011)

odysseus:hades oswald$ svn info ../dune-grid
Path: ../dune-grid
URL: https://svn.dune-project.org/svn/dune-grid/trunk
Repository Root: https://svn.dune-project.org/svn/dune-grid
Repository UUID: f82f27de-4411-0410-b156-930a902cda50
Revision: 7671
Node Kind: directory
Schedule: normal
Last Changed Author: mnolte
Last Changed Rev: 7671
Last Changed Date: 2011-09-03 15:55:55 +0200 (Sat, 03 Sep 2011)

odysseus:hades oswald$ svn info ../dune-grid-howto/
Path: ../dune-grid-howto
URL: https://svn.dune-project.org/svn/dune-grid-howto/trunk
Repository Root: https://svn.dune-project.org/svn/dune-grid-howto
Repository UUID: 971b7fa1-7d12-0410-b1b7-9c127676b1f7
Revision: 360
Node Kind: directory
Schedule: normal
Last Changed Author: sander
Last Changed Rev: 358
Last Changed Date: 2011-05-26 17:50:48 +0200 (Thu, 26 May 2011)

odysseus:hades oswald$ svn info ../dune-istl/
Path: ../dune-istl
URL: https://svn.dune-project.org/svn/dune-istl/trunk
Repository Root: https://svn.dune-project.org/svn/dune-istl
Repository UUID: 8c7282f0-4411-0410-8e4f-cdd3e049afb8
Revision: 1495
Node Kind: directory
Schedule: normal
Last Changed Author: sander
Last Changed Rev: 1495
Last Changed Date: 2011-08-30 11:03:49 +0200 (Tue, 30 Aug 2011)

odysseus:hades oswald$ gcc --v
Using built-in specs.
Target: x86_64-apple-darwin10.4.0
Configured with: ./configure --prefix=/Users/oswald/extlib/gcc/gcc-4.4.5 --with-gmp=/Users/oswald/extlib/gmp-5.0.1 --with-mpfr=/Users/oswald/extlib/mpfr-3.0 --with-mpc=/Users/oswald/extlib/mpc-0.8.2 --with-ppl=/Users/oswald/extlib/ppl-0.11 --with-cloog=/Users/oswald/extlib/cloog-ppl-trunk
Thread model: posix
gcc version 4.4.5 (GCC) 
odysseus:hades oswald$ uname .a
usage: uname [-amnprsv]
odysseus:hades oswald$ uname -a
Darwin odysseus.psi.ch 10.8.0 Darwin Kernel Version 10.8.0: Tue Jun  7 16:33:36 PDT 2011; root:xnu-1504.15.3~1/RELEASE_I386 i386
odysseus:hades oswald$ 



--------------------------------------------------------------------------------------------------------------------------------------------------------------------
Benedikt Oswald, Dr. sc. techn., dipl. El. Ing. ETH, www.psi.ch, Computational Accelerator Scientist, http://amas.web.psi.ch
Paul Scherrer  Institute (PSI), CH-5232 Villigen, Suisse, benedikt.oswald at psi.ch, +41(0)56 310 32 12

"Passion is required for any great work, and for the revolution passion and audacity are required in big doses.",
 Ernesto 'Che' Guevara, Letter to his parents.
--------------------------------------------------------------------------------------------------------------------------------------------------------------------

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20110906/89ef031a/attachment.htm>