<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7655.11">
<TITLE>RE: [Dune] AluCube grid in 2d?</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>

<P><FONT SIZE=2>The seem to be using a dgf files. To avoid the closure you can also add<BR>
the following to your dgf file<BR>
<BR>
GridParameter<BR>
closure none<BR>
#<BR>
<BR>
That should do what Carsten suggested but would also work with other GridType<BR>
(would just be ignored for the new ALUGrid for example since that grid manager can<BR>
only do non conforming refinement).<BR>
<BR>
Best<BR>
Andreas<BR>
<BR>
-----Original Message-----<BR>
From: dune-bounces+a.s.dedner=warwick.ac.uk@dune-project.org on behalf of Matteo Semplice<BR>
Sent: Thu 11/24/2011 9:36 PM<BR>
To: dune@dune-project.org<BR>
Subject: Re: [Dune] AluCube grid in 2d?<BR>
<BR>
Sorry to bother ouy again. I got ug-grid and recompiled dune.<BR>
<BR>
Now compiling the example with<BR>
<BR>
make GRIDDIM=2 GRIDTYPE=UGGRID adaptivefinitevolume<BR>
<BR>
produces a solution on a mesh with a mixture of squares and trapezoids:<BR>
a big square is divided in four little squares or in 3 cells of which a<BR>
square and 2 trapezoids. I tried to picture this below, but I am not an<BR>
ASCII-art expert...<BR>
<BR>
|--------|<BR>
|  \     |<BR>
|   |----|<BR>
|---|----|<BR>
<BR>
Is there a way to force it to use only squares?<BR>
<BR>
     Matteo<BR>
<BR>
PS when compiling I get the warning<BR>
<BR>
In file included from<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/io/file/dgfparser/dgfug.hh:145,<BR>
                  from config.h:390,<BR>
                  from adaptivefinitevolume.cc:1:<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/io/file/dgfparser/dgfug.cc:<BR>
In static member function 'static Dune::UGGrid<dim>*<BR>
Dune::MacroGrid::Impl<Dune::UGGrid<dim> >::generate(Dune::MacroGrid&,<BR>
const char*, ompi_communicator_t*) [with int dim = 2]':<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/io/file/dgfparser/dgfparser.hh:141:  <BR>
instantiated from 'GridType* Dune::MacroGrid::createGrid() [with<BR>
GridType = Dune::UGGrid<2>]'<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/io/file/dgfparser/dgfparser.hh:210:  <BR>
instantiated from 'Dune::DGFGridFactory<GridType>::DGFGridFactory(const<BR>
std::string&, ompi_communicator_t*) [with G = Dune::UGGrid<2>]'<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/io/file/dgfparser/dgfparser.hh:354:  <BR>
instantiated from 'Dune::GridPtr<GridType>::GridPtr(const std::string&,<BR>
ompi_communicator_t*) [with GridType = Dune::UGGrid<2>]'<BR>
adaptivefinitevolume.cc:107:   instantiated from here<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/io/file/dgfparser/dgfug.cc:28:<BR>
warning: 'Dune::UGGrid<dim>::UGGrid(unsigned int) [with int dim = 2]' is<BR>
deprecated (declared at<BR>
/home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/uggrid.hh:427)<BR>
<BR>
On 24/11/2011 17:27, Oliver Sander wrote:<BR>
> Use UGGrid if you want 2d cube grids.<BR>
> --<BR>
> Oliver<BR>
><BR>
> Am 24.11.2011 17:21, schrieb Andreas Dedner:<BR>
>   <BR>
>> The error is correct at least with alu 1.23. There's a alu 1.5 on the<BR>
>> way which includes cubes in 2d and surface grids 2d/3d - should be<BR>
>> released within a week or two.<BR>
>><BR>
>> Andreas<BR>
>><BR>
>> On 24/11/11 16:13, Matteo Semplice wrote:<BR>
>>     <BR>
>>> Hi.<BR>
>>> Does the error reported below point to some mistake in my library<BR>
>>> setup/compilation or is it the intended behaviour (i.e. dune cannot<BR>
>>> handle ALUCube grid in 2d?)<BR>
>>><BR>
>>> So, I downloaded ALUGrid-1.23 , configured it with<BR>
>>> ./configure --prefix=/home/matteo/<BR>
>>> --includedir=/home/matteo/include/alugrid<BR>
>>> and installed with "make all" and "make install". No errors were<BR>
>>> reported.<BR>
>>><BR>
>>> Then I downloaded dune-common, grid, istl and grid-howto, and<BR>
>>> configured, compiled, installed with dunecontrol.opts<BR>
>>><BR>
>>> PKG_CONFIG_PATH="/home/matteo/lib/pkgconfig"<BR>
>>> CONFIGURE_FLAGS="--enable-parallel<BR>
>>> --with-alugrid=/home/matteo/ALUGrid-1.23 --prefix=$HOME"<BR>
>>> (There is no newline in the definition of CONFIGURE_FLAGS in the actual<BR>
>>> file)<BR>
>>><BR>
>>> Both the previously compiled alugrid libraries and the alberta libraries<BR>
>>> already in my system were found. Again, I noticed no errors.<BR>
>>><BR>
>>> Then, I moved to the dune-grid-2.1.0 directory and tried to recompile<BR>
>>> adaptivefinitevolume with GRIDTYPE set to ALUGRID_SIMPLEX and<BR>
>>> ALUGRID_CONFORM, as suggested in the grid-howto. They compile and run<BR>
>>> correctly.<BR>
>>><BR>
>>> Next I tried ALUGRID_CUBE and got the following error.<BR>
>>> I noticed that the error mentions dune/grid/alugrid/2d/alugrid.hh,<BR>
>>> dune/grid/alugrid.hh and lugrid/3d/capabilities.hh... is this ok?<BR>
>>><BR>
>>> Thanks for any tip...<BR>
>>><BR>
>>> Matteo<BR>
>>><BR>
>>> ==========<BR>
>>> matteo@montebianco:~/dune2.1/dune-grid-howto-2.1.0$ make GRIDDIM=2<BR>
>>> GRIDTYPE=ALUGRID_CUBE adaptivefinitevolume<BR>
>>> g++ -std=c++0x -DHAVE_CONFIG_H -I.<BR>
>>> -I/home/matteo/dune2.1/dune-common-2.1.0<BR>
>>> -I/home/matteo/dune2.1/dune-grid-2.1.0<BR>
>>> -I/home/matteo/dune2.1/dune-istl-2.1.0 -I. -I/usr/include<BR>
>>> -I/usr/lib/openmpi/include -I/usr/lib/openmpi/include/openmpi -pthread<BR>
>>> -DMPIPP_H -DENABLE_MPI=1 -I/usr/include<BR>
>>> -I/home/matteo/dune2.1/dune-common-2.1.0<BR>
>>> -I/home/matteo/dune2.1/dune-grid-2.1.0<BR>
>>> -I/home/matteo/dune2.1/dune-istl-2.1.0 -I/usr/lib/openmpi/include<BR>
>>> -I/usr/lib/openmpi/include/openmpi -pthread -DMPIPP_H -DENABLE_MPI=1<BR>
>>> -I/home/matteo/dune2.1/dune-common-2.1.0 -DGRIDDIM=2 -DWORLDDIM=2<BR>
>>> -DALUGRID_CUBE -I/usr/include/alberta -DALBERTA_DIM=2 -DENABLE_ALBERTA<BR>
>>> -I/usr/lib/openmpi/include -I/usr/lib/openmpi/include/openmpi -pthread<BR>
>>> -DMPIPP_H -DENABLE_MPI=1 -I/home/matteo/include/alugrid<BR>
>>> -I/home/matteo/include/alugrid/serial<BR>
>>> -I/home/matteo/include/alugrid/duneinterface -DENABLE_ALUGRID<BR>
>>> -I/home/matteo/include/alugrid/parallel -I/usr/include -DENABLE_GMP=1<BR>
>>> -I/home/matteo/dune2.1/dune-grid-2.1.0 -I/usr/include/superlu<BR>
>>> -DENABLE_SUPERLU -I/home/matteo/dune2.1/dune-istl-2.1.0 -g -O2 -MT<BR>
>>> adaptivefinitevolume-adaptivefinitevolume.o -MD -MP -MF<BR>
>>> .deps/adaptivefinitevolume-adaptivefinitevolume.Tpo -c -o<BR>
>>> adaptivefinitevolume-adaptivefinitevolume.o `test -f<BR>
>>> 'adaptivefinitevolume.cc' || echo './'`adaptivefinitevolume.cc<BR>
>>> In file included from<BR>
>>> /home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/alugrid/2d/alugrid.hh:8,<BR>
>>> from /home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/alugrid.hh:11,<BR>
>>> from config.h:457,<BR>
>>> from adaptivefinitevolume.cc:1:<BR>
>>> /home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/alugrid/2d/grid.hh: In<BR>
>>> instantiation of 'Dune::ALU2dGrid<2, 2, (ALUGridSpace::ElementType)1u>':<BR>
>>> /home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/alugrid/3d/capabilities.hh:21:<BR>
>>><BR>
>>> instantiated from 'Dune::ALUCubeGrid<2, 2>'<BR>
>>> adaptivefinitevolume.cc:104: instantiated from here<BR>
>>> /home/matteo/dune2.1/dune-grid-2.1.0/dune/grid/alugrid/2d/grid.hh:192:<BR>
>>> error: static assertion failed: "ALU2dGrid only implemented for<BR>
>>> triangles."<BR>
>>> make: *** [adaptivefinitevolume-adaptivefinitevolume.o] Error 1<BR>
>>> ===============<BR>
>>><BR>
>>> _______________________________________________<BR>
>>> Dune mailing list<BR>
>>> Dune@dune-project.org<BR>
>>> <A HREF="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</A><BR>
>>>       <BR>
>><BR>
>> _______________________________________________<BR>
>> Dune mailing list<BR>
>> Dune@dune-project.org<BR>
>> <A HREF="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</A><BR>
>>     <BR>
> _______________________________________________<BR>
> Dune mailing list<BR>
> Dune@dune-project.org<BR>
> <A HREF="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</A><BR>
>   <BR>
<BR>
<BR>
--<BR>
Matteo Semplice                         Dip. di Fisica e Matematica<BR>
Phone: 031-2386132                      Università dell'Insubria<BR>
Fax:   031-2386209                      Via Valleggio, 11<BR>
                                        22100 Como<BR>
<BR>
<BR>
_______________________________________________<BR>
Dune mailing list<BR>
Dune@dune-project.org<BR>
<A HREF="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</A><BR>
<BR>
</FONT>
</P>

</BODY>
</HTML>