matteo@montebianco:~/svn/dune/2.2rc/dune-fv/src$ ./test
GridParameterBlock: Parameter 'refinementedge' not specified, defaulting to 'ARBITRARY'.
check_mesh:         checking mesh 'DUNE AlbertaGrid'
                    checking done; no error detected
Cell at x= 0.333333 0.166667
-Face at local coordinate 0.5 0.5
  face.geometry().center()= 0.25 0.25
  face.geometryInOutside().center()= 0.5 0.5
  face.geometryInInside().center()= 0.5 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.166667 0.333333
   using face.geometryInOutside trick = 0.166667 0.333333
-Face at local coordinate 0 0.5
  face.geometry().center()= 0.25 0
  face.geometryInOutside().center()= 0 0.5
  face.geometryInInside().center()= 0 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.166667 0.833333    <==== wrong...
   using face.geometryInOutside trick = 0.166667 -0.166667   <==== GOOD!
 (is a periodic boundary!)
-Face at local coordinate 0.5 0
  face.geometry().center()= 0.5 0.25
  face.geometryInOutside().center()= 0.5 0
  face.geometryInInside().center()= 0.5 0
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.666667 0.333333
   using face.geometryInOutside trick = 0.666667 0.333333

matteo@montebianco:~/svn/dune/2.2rc/dune-fv/src$ ./test
GridParameterBlock: Parameter 'refinementedge' not specified, defaulting to 'ARBITRARY'.
macro_test_3d:      WARNING in ./../3d/macro_3d.c, line 368 
                    not implemented for 3d yet: no check is performed
check_mesh:         checking mesh 'DUNE AlbertaGrid'           
                    checking done; no error detected
Cell at x= 0.375 0.25 0.125                 
-Face at local coordinate 0.333333 0.333333 0.333333  
  face.geometry().center()= 0.5 0.333333 0.166666
  face.geometryInOutside().center()= 0 0.333333 0.333333       
  face.geometryInInside().center()= 0.333333 0.333333 0.333333   
-- global coordinate of nbhr:                            
   as reported by outside.geometry()  = 0.625 0.375 0.25  
   using face.geometryInOutside trick = 0.625 0.375 0.25 
-Face at local coordinate 0 0.333333 0.333333       
  face.geometry().center()= 0.333333 0.166667 0       
  face.geometryInOutside().center()= 0.333333 0.333333 0.333333   
  face.geometryInInside().center()= 0 0.333333 0.333333
-- global coordinate of nbhr:                  
   as reported by outside.geometry()  = 0.25 0.125 0.875    <=== wrong...  
   using face.geometryInOutside trick = 0.25 0.125 -0.125   <=== GOOD!
 (is a periodic boundary!)             
-Face at local coordinate 0.333333 0 0.333333   
  face.geometry().center()= 0.333333 0.166667 0.166667 
  face.geometryInOutside().center()= 0.333333 0 0.333333
  face.geometryInInside().center()= 0.333333 0 0.333333
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.375 0.125 0.25
   using face.geometryInOutside trick = 0.375 0.125 0.25
-Face at local coordinate 0.333333 0.333333 0
  face.geometry().center()= 0.333333 0.333333 0.166667
  face.geometryInOutside().center()= 0.333333 0.333333 0
  face.geometryInInside().center()= 0.333333 0.333333 0
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.25 0.375 0.125
   using face.geometryInOutside trick = 0.25 0.375 0.125


matteo@montebianco:~/svn/dune/2.2rc/dune-fv/src$ ./test
GridParameterBlock: Parameter 'tolerance' not specified, defaulting to `1e-08'!

Created serial ALUCubeGrid<2,2> from macro grid file 'grids/periodicunitcube2.dgf'. 

Cell at x= 0.25 0.25
-Face at local coordinate 1 0.5
  face.geometry().center()= 0.5 0.25
  face.geometryInOutside().center()= 0 0.5
  face.geometryInInside().center()= 1 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.75 0.25
   using face.geometryInOutside trick = 0.75 0.25
-Face at local coordinate 0.5 1
  face.geometry().center()= 0.25 0.5
  face.geometryInOutside().center()= 0.5 0
  face.geometryInInside().center()= 0.5 1
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.25 0.75
   using face.geometryInOutside trick = 0.25 0.75
-Face at local coordinate 0 0.5
  face.geometry().center()= 0 0.25
  face.geometryInOutside().center()= 0 0.5       <==== should be [1,.5]
  face.geometryInInside().center()= 0 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.75 0.25 <=== wrong...
   using face.geometryInOutside trick = 0.25 0.25 <=== wrong
 (is a periodic boundary!)
-Face at local coordinate 0.5 0
  face.geometry().center()= 0.25 0
  face.geometryInOutside().center()= 0.5 0       <==== should be [.5,1]
  face.geometryInInside().center()= 0.5 0
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.25 0.75 <=== wrong...
   using face.geometryInOutside trick = 0.25 0.25 <=== wrong
 (is a periodic boundary!)


matteo@montebianco:~/svn/dune/2.2rc/dune-fv/src$ ./test 
Created serial ALUCubeGrid<3,3> from macro grid file 'grids/periodicunitcube3.dgf'. 

test: /home/matteo/svn/dune/2.2rc/dune-grid/dune/grid/alugrid/3d/faceutility_imp.cc:215: void Dune::ALU3dGridFaceInfo<type, Comm>::updateFaceInfo(const typename Dune::ALU3dImplTraits<elType, Comm>::GEOFaceType&, int, int) [with Dune::ALU3dGridElementType type = (Dune::ALU3dGridElementType)7u, Comm = Dune::No_Comm]: Assertion `bndType_ == periodicBoundary || bndType_ == domainBoundary ? bndId_ > 0 : bndId_ == 0' failed.
Abortito

matteo@montebianco:~/svn/dune/2.2rc/dune-fv/src$ ./test
Hello World! This is dune-fv.
This is a sequential program.
GridParameterBlock: Parameter 'tolerance' not specified, defaulting to `1e-08'!

Created serial ALUSimplexGrid<2,2> from macro grid file 'grids/periodicunitcube2.dgf'. 

Cell at x= 0.333333 0.166667
-Face at local coordinate 0.5 0.5
  face.geometry().center()= 0.25 0.25
  face.geometryInOutside().center()= 0.5 0.5
  face.geometryInInside().center()= 0.5 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.166667 0.333333
   using face.geometryInOutside trick = 0.166667 0.333333
-Face at local coordinate 0 0.5
  face.geometry().center()= 0.25 0
  face.geometryInOutside().center()= 0 0.5
  face.geometryInInside().center()= 0 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.166667 0.833333   <=== wrong...
   using face.geometryInOutside trick = 0.166667 -0.166667  <=== GOOD!
 (is a periodic boundary!)
-Face at local coordinate 0.5 0
  face.geometry().center()= 0.5 0.25
  face.geometryInOutside().center()= 0.5 0
  face.geometryInInside().center()= 0.5 0
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.666667 0.333333
   using face.geometryInOutside trick = 0.666667 0.333333

matteo@montebianco:~/svn/dune/2.2rc/dune-fv/src$ ./test

Created serial ALUSimplexGrid<3,3> from macro grid file 'grids/periodicunitcube3.dgf'. 

test: /home/matteo/svn/dune/2.2rc/dune-grid/dune/grid/alugrid/3d/faceutility_imp.cc:215: void Dune::ALU3dGridFaceInfo<type, Comm>::updateFaceInfo(const typename Dune::ALU3dImplTraits<elType, Comm>::GEOFaceType&, int, int) [with Dune::ALU3dGridElementType type = (Dune::ALU3dGridElementType)4u, Comm = Dune::No_Comm]: Assertion `bndType_ == periodicBoundary || bndType_ == domainBoundary ? bndId_ > 0 : bndId_ == 0' failed.
Abortito

YaspGridParameterBlock: Parameter 'overlap' not specified, defaulting to '0'.
Cell at x= 0.25 0.25
-Face at local coordinate 0 0.5
-- is on domain boundary!                  <==== OUCH! WE ASKED FOR PERIODIC!
-Face at local coordinate 1 0.5
  face.geometry().center()= 0.5 0.25
  face.geometryInOutside().center()= 0 0.5
  face.geometryInInside().center()= 1 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.75 0.25
   using face.geometryInOutside trick = 0.75 0.25
-Face at local coordinate 0.5 0
-- is on domain boundary!                  <==== OUCH! WE ASKED FOR PERIODIC!
-Face at local coordinate 0.5 1
  face.geometry().center()= 0.25 0.5
  face.geometryInOutside().center()= 0.5 0
  face.geometryInInside().center()= 0.5 1
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.25 0.75
   using face.geometryInOutside trick = 0.25 0.75

YaspGridParameterBlock: Parameter 'overlap' not specified, defaulting to '0'.
Cell at x= 0.25 0.25 0.25
-Face at local coordinate 0 0.5 0.5
-- is on domain boundary!                  <==== OUCH! WE ASKED FOR PERIODIC!
-Face at local coordinate 1 0.5 0.5
  face.geometry().center()= 0.5 0.25 0.25
  face.geometryInOutside().center()= 0 0.5 0.5
  face.geometryInInside().center()= 1 0.5 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.75 0.25 0.25
   using face.geometryInOutside trick = 0.75 0.25 0.25
-Face at local coordinate 0.5 0 0.5
-- is on domain boundary!                  <==== OUCH! WE ASKED FOR PERIODIC!
-Face at local coordinate 0.5 1 0.5
  face.geometry().center()= 0.25 0.5 0.25
  face.geometryInOutside().center()= 0.5 0 0.5
  face.geometryInInside().center()= 0.5 1 0.5
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.25 0.75 0.25
   using face.geometryInOutside trick = 0.25 0.75 0.25
-Face at local coordinate 0.5 0.5 0
-- is on domain boundary!                  <==== OUCH! WE ASKED FOR PERIODIC!
-Face at local coordinate 0.5 0.5 1
  face.geometry().center()= 0.25 0.25 0.5
  face.geometryInOutside().center()= 0.5 0.5 0
  face.geometryInInside().center()= 0.5 0.5 1
-- global coordinate of nbhr:
   as reported by outside.geometry()  = 0.25 0.25 0.75
   using face.geometryInOutside trick = 0.25 0.25 0.75

With GRIDTYPE=SGRID (both 2d and 3d), same as GRIDTYPE=YASPGRID
With FORCEYASP (constructs the periodic grid without dgfparser), same
as before