Anuga-user Digest, Vol 42, Issue 3

10 views
Skip to first unread message

anuga-use...@lists.sourceforge.net

unread,
Feb 14, 2012, 5:10:41 AM2/14/12
to anuga...@lists.sourceforge.net
Send Anuga-user mailing list submissions to
anuga...@lists.sourceforge.net

To subscribe or unsubscribe via the World Wide Web, visit
https://lists.sourceforge.net/lists/listinfo/anuga-user
or, via email, send a message with subject or body 'help' to
anuga-use...@lists.sourceforge.net

You can reach the person managing the list at
anuga-us...@lists.sourceforge.net

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Anuga-user digest..."


Today's Topics:

1. Mesh file (Mathieu Lap?tre)
2. Re: Mesh file (Stephen Roberts)
3. Adding an Interior Hole to geometry (marcor...@libero.it)
4. Re: Adding an Interior Hole to geometry (Stephen Roberts)


----------------------------------------------------------------------

Message: 1
Date: Fri, 10 Feb 2012 20:03:39 +0100
From: Mathieu Lap?tre <lapotre...@hotmail.com>
Subject: [Anuga-user] Mesh file
To: <anuga...@lists.sourceforge.net>
Message-ID: <DUB114-W63E4B5E1D...@phx.gbl>
Content-Type: text/plain; charset="iso-8859-1"


Hello,
Is there a simple way to visualize the mesh from the .msh file instead of the .sww file on Windows?Cheers,
Mathieu Lapotre
-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 2
Date: Mon, 13 Feb 2012 09:32:42 +1100
From: Stephen Roberts <stephen...@anu.edu.au>
Subject: Re: [Anuga-user] Mesh file
To: anuga...@lists.sourceforge.net
Message-ID: <4F383E0A...@anu.edu.au>
Content-Type: text/plain; charset="iso-8859-1"

Hi Mathieu,

There is a primitive mesh viewer available in the anuga distribution.

It is in the anuga/pmesh directory.

From that directory run

python graphical_mesh_generator.py


As the name implies, you can create a msh file, but you can also open a
msh file.

Doesn't work well with large meshes though.

Cheers
Steve


On 11/02/12 06:03, Mathieu Lap?tre wrote:
> Hello,
>
> Is there a simple way to visualize the mesh from the .msh file instead
> of the .sww file on Windows?
> Cheers,
>
> Mathieu Lapotre
>
>
> ------------------------------------------------------------------------------
> Virtualization& Cloud Management Using Capacity Planning
> Cloud computing makes use of virtualization - but cloud computing
> also focuses on allowing computing to be delivered as a service.
> http://www.accelacomm.com/jaw/sfnl/114/51521223/
>
>
> _______________________________________________
> Anuga-user mailing list
> Anuga...@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/anuga-user


--
+--------------------------------------------------------------------+
| Steve Roberts | My Office: (61)(2) 6125 4445 |
| Department of Mathematics | Math Office: (61)(2) 6125 2908 |
| Mathematical Sciences Institute| Fax: (61)(2) 6125 4984 |
| John Dedman Building #27 | mailto:stephen...@anu.edu.au |
| Australian National University | http://www.maths.anu.edu.au/~steve|
| Canberra ACT 0200 AUSTRALIA | ANU CRICOS # 00120C |
+--------------------------------------------------------------------+

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 3
Date: Tue, 14 Feb 2012 09:28:11 +0100
From: "marcor...@libero.it" <marcor...@libero.it>
Subject: [Anuga-user] Adding an Interior Hole to geometry
To: anuga...@lists.sourceforge.net
Message-ID: <4F3A1B1B...@libero.it>
Content-Type: text/plain; charset="iso-8859-1"

Hi All,

I am trying to simulate the presence of a building within my inundation
area. I therefore tried to insert an internal hole, like this:

U = [[686072,923091], [686072,923084], [686082,923084], [686082,923091]]
# building perimeter
domain = anuga.create_domain_from_regions(project.bounding_polygon,
boundary_tags={'nord': [0], 'est': [1], 'sud': [2], 'ovest': [3]},
maximum_triangle_area=project.default_res,
mesh_filename=project.meshname, interior_holes=U,
use_cache=True,
verbose=True)

Unfortunately I get the following error message:

Traceback (most recent call last):
File "ammonite.py", line 47, in <module>
verbose=True)
File "/usr/lib/python2.6/dist-packages/anuga/__init__.py", line 254,
in create_domain_from_regions
compression=False)
File "/usr/lib/python2.6/dist-packages/anuga/caching/caching.py",
line 376, in cache
T = my_F(*args, **kwargs) # Built-in 'apply' deprecated in Py3K
File "/usr/lib/python2.6/dist-packages/anuga/__init__.py", line 292,
in _create_domain_from_regions
verbose=verbose)
File
"/usr/lib/python2.6/dist-packages/anuga/pmesh/mesh_interface.py", line
121, in create_mesh_from_regions
args, kwargs)
File
"/usr/lib/python2.6/dist-packages/anuga/pmesh/mesh_interface.py", line
241, in _create_mesh_from_regions
raise PolygonError, msg
PolygonError: Interior polygon [686072, 923091] is outside bounding
polygon: [[ 686007.82 923006.11]
[ 686145.11 923006.11]
[ 686145.11 923145.11]
[ 686003.91 923145.6 ]]

The building perimeter is actually inside, not outside, the bounding
polygon, but Anuga seems to think the opposite is true. Can you suggest
the reason of the error?
Thank you in advance!

Regards,
Marco

-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

Message: 4
Date: Tue, 14 Feb 2012 21:10:22 +1100
From: Stephen Roberts <stephen...@anu.edu.au>
Subject: Re: [Anuga-user] Adding an Interior Hole to geometry
To: "marcor...@libero.it" <marcor...@libero.it>,
anuga...@lists.sourceforge.net
Message-ID: <4F3A330E...@anu.edu.au>
Content-Type: text/plain; charset="iso-8859-1"

Hi Marco,

Sure seems that [686072, 923091] is inside the bounding polygon: [[
686007.82 923006.11]
[ 686145.11 923006.11]
[ 686145.11 923145.11]
[ 686003.91 923145.6 ]]

Could you provide me with your scripts. I would be interested to see
where the problem arises.

By the way, the interior_holes should be a list of polygons, so you
should pass the hole via

interior_holes=[U]

Ie make U into a list containing U, ie [U]

Cheers
Steve


On 14/02/2012 7:28 PM, marcor...@libero.it wrote:
>
> Hi All,
>
> I am trying to simulate the presence of a building within my
> inundation area. I therefore tried to insert an internal hole, like this:
>
> U = [[686072,923091], [686072,923084], [686082,923084],
> [686082,923091]] # building perimeter
> domain = anuga.create_domain_from_regions(project.bounding_polygon,
> boundary_tags={'nord': [0], 'est': [1], 'sud': [2], 'ovest': [3]},
> maximum_triangle_area=project.default_res,
> mesh_filename=project.meshname, interior_holes=U,
> use_cache=True,
> verbose=True)
>
> Unfortunately I get the following error message:
>
> Traceback (most recent call last):
> File "ammonite.py", line 47, in <module>
> verbose=True)
> File "/usr/lib/python2.6/dist-packages/anuga/__init__.py", line 254,
> in create_domain_from_regions
> compression=False)
> File "/usr/lib/python2.6/dist-packages/anuga/caching/caching.py",
> line 376, in cache
> T = my_F(*args, **kwargs) # Built-in 'apply' deprecated in Py3K
> File "/usr/lib/python2.6/dist-packages/anuga/__init__.py", line 292,
> in _create_domain_from_regions
> verbose=verbose)
> File
> "/usr/lib/python2.6/dist-packages/anuga/pmesh/mesh_interface.py", line
> 121, in create_mesh_from_regions
> args, kwargs)
> File
> "/usr/lib/python2.6/dist-packages/anuga/pmesh/mesh_interface.py", line
> 241, in _create_mesh_from_regions
> raise PolygonError, msg
> PolygonError: Interior polygon [686072, 923091] is outside bounding
> polygon: [[ 686007.82 923006.11]
> [ 686145.11 923006.11]
> [ 686145.11 923145.11]
> [ 686003.91 923145.6 ]]
>
> The building perimeter is actually inside, not outside, the bounding
> polygon, but Anuga seems to think the opposite is true. Can you
> suggest the reason of the error?
> Thank you in advance!
>
> Regards,
> Marco
>
>
>
> ------------------------------------------------------------------------------
> Keep Your Developer Skills Current with LearnDevNow!
> The most comprehensive online learning library for Microsoft developers
> is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5, CSS3, MVC3,
> Metro Style Apps, more. Free future releases when you subscribe now!
> http://p.sf.net/sfu/learndevnow-d2d
>
>
> _______________________________________________
> Anuga-user mailing list
> Anuga...@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/anuga-user


--
+--------------------------------------------------------------------+
| Steve Roberts | My Office: (61)(2) 6125 4445 |
| Department of Mathematics | Math Office: (61)(2) 6125 2908 |
| Mathematical Sciences Institute| Fax: (61)(2) 6125 4984 |
| John Dedman Building #27 | mailto:stephen...@anu.edu.au |
| Australian National University | http://www.maths.anu.edu.au/~steve|
| Canberra ACT 0200 AUSTRALIA | ANU CRICOS # 00120C |
+--------------------------------------------------------------------+


-------------- next part --------------
An HTML attachment was scrubbed...

------------------------------

------------------------------------------------------------------------------
Keep Your Developer Skills Current with LearnDevNow!
The most comprehensive online learning library for Microsoft developers
is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5, CSS3, MVC3,
Metro Style Apps, more. Free future releases when you subscribe now!
http://p.sf.net/sfu/learndevnow-d2d

------------------------------

_______________________________________________
Anuga-user mailing list
Anuga...@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/anuga-user


End of Anuga-user Digest, Vol 42, Issue 3
*****************************************

Reply all
Reply to author
Forward
0 new messages