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

Tue Sep 6 16:56:52 CEST 2011

Hi Benedikt,

the block looks perfectly ok to me. Could you please file a decent bug report
including the DGF file?

Best,

Martin

On 09/06/2011 03:45 PM, Benedikt Oswald wrote:
> 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 <http://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
> <http://www.psi.ch>, Computational Accelerator Scientist, http://amas.web.psi.ch
> Paul Scherrer  Institute (PSI), CH-5232 Villigen, Suisse,
> benedikt.oswald at psi.ch <mailto: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.
> --------------------------------------------------------------------------------------------------------------------------------------------------------------------
> 
> 
> 
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune

-- 
Martin Nolte <nolte at mathematik.uni-freiburg.de>

Universität Freiburg                                   phone: +49-761-203-5630
Abteilung für angewandte Mathematik                    fax:   +49-761-203-5632
Hermann-Herder-Straße 10
79104 Freiburg, Germany