Ignore:
Timestamp:
Nov 8, 2006, 5:35:23 PM (17 years ago)
Author:
ole
Message:

One large step towards major cleanup. This has mainly to do with
the way vertex coordinates are handled internally.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/utilities/quad.py

    r3566 r3945  
    139139        if len(args) == 2:
    140140            point_id = int(args[1])
    141             x = self.__class__.mesh.coordinates[point_id][0]
    142             y = self.__class__.mesh.coordinates[point_id][1]
     141            x, y = self.__class__.mesh.get_nodes()[point_id]
    143142
    144143            #print point_id, x, y
     
    394393    #Make root cell
    395394    #print mesh.coordinates
    396    
    397     xmin = min(mesh.coordinates[:,0])
    398     xmax = max(mesh.coordinates[:,0])
    399     ymin = min(mesh.coordinates[:,1])
    400     ymax = max(mesh.coordinates[:,1])
     395
     396    nodes = mesh.get_nodes()
     397    xmin = min(nodes[:,0])
     398    xmax = max(nodes[:,0])
     399    ymin = min(nodes[:,1])
     400    ymax = max(nodes[:,1])
    401401
    402402   
     
    425425   
    426426    #Insert indices of all vertices
    427     root.insert( range(len(mesh.coordinates)) )
     427    root.insert( range(mesh.number_of_nodes) )
    428428
    429429    #Build quad tree and return
Note: See TracChangeset for help on using the changeset viewer.