Changeset 297


Ignore:
Timestamp:
Sep 14, 2004, 7:20:30 PM (21 years ago)
Author:
ole
Message:

Adapted original pmesh_to_domain

Location:
inundation/ga/storm_surge/pyvolution
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • inundation/ga/storm_surge/pyvolution/HUSK.txt

    r281 r297  
    11
    2 Do datamanager for non-smooth data + tests. Smooth is OK (7/9)
     2Do data manager for non-smooth data + tests. Smooth is OK (7/9)
    33               
    44       
    55Perhaps control old/new style with minimum_allowed height
    6 In any case get rid of oldstyle.
     6In any case get rid of old style.
    77Then optimise protect_against.....
    88       
  • inundation/ga/storm_surge/pyvolution/mesh.py

    r291 r297  
    307307
    308308    def build_boundary_dictionary(self, boundary = None):
    309         """Build or input the dictionary of boundary tags.
     309        """Build or check the dictionary of boundary tags.
    310310         self.boundary is a dictionary of tags,
    311311         keyed by volume id and edge:
     
    448448
    449449
     450        #Check that all vertices have been registered
     451        for v in self.vertexlist:
     452            msg = 'Some points do not belong to an element.\n'
     453            msg += 'Make sure all points appear as element vertices!'
     454            if v is None:
     455                print 'WARNING (mesh.py): %s' %msg
     456                break
     457
     458            #Should this warning be an exception?
     459            #assert v is not None, msg
     460
     461               
    450462        #Check integrity of neighbour structure
    451463        for i in range(N):
  • inundation/ga/storm_surge/pyvolution/netherlands.py

    r288 r297  
    9393N = 600 #Size = 720000
    9494N = 220
     95N = 55
    9596
    9697
  • inundation/ga/storm_surge/pyvolution/quantity.py

    r292 r297  
    295295                #Go through all triangle, vertex pairs
    296296                #contributing to vertex k and register vertex value
    297                 contributions = []               
     297
     298                if L is None: continue #In case there are unused points
     299               
     300                contributions = []
    298301                for volume_id, vertex_id in L:
    299 
    300302                    v = self.vertex_values[volume_id, vertex_id]
    301303                    contributions.append(v)
  • inundation/ga/storm_surge/pyvolution/test_mesh.py

    r258 r297  
    8080        #General tests:
    8181
    82         #Test that points are arranged in a counter clock wise order
     82        #Test that points are arranged in a counter clock wise order etc
    8383        mesh.check_integrity()
    8484
  • inundation/ga/storm_surge/pyvolution/test_quantity.py

    r282 r297  
    2626        elements = [ [1,0,2], [1,2,4], [4,2,5], [3,1,4] ]
    2727
    28         self.mesh1 = Domain(points, [elements[0]])       
     28        self.mesh1 = Domain(points[:3], [elements[0]])       
    2929        self.mesh1.check_integrity()
    3030       
Note: See TracChangeset for help on using the changeset viewer.