Changeset 5951


Ignore:
Timestamp:
Nov 12, 2008, 3:16:51 PM (16 years ago)
Author:
rwilson
Message:

More NumPy? changes.

Location:
anuga_core/source_numpy_conversion/anuga
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • anuga_core/source_numpy_conversion/anuga/load_mesh/test_loadASCII.py

    r5910 r5951  
    179179        #print "*(*("
    180180       
    181         self.failUnless(numpy.array(meshDict['vertices'])  ==
    182                         numpy.array(loadedDict['vertices']),
    183                          'test_export_mesh_file failed. Test 1')
    184         self.failUnless(numpy.array(meshDict['triangles'])  ==
    185                         numpy.array(loadedDict['triangles']),
    186                          'test_export_mesh_file failed. Test 2')
    187         self.failUnless(numpy.array(meshDict['segments'])  ==
    188                         numpy.array(loadedDict['segments']),
    189                          'test_export_mesh_file failed. Test 3')
    190         self.failUnless(numpy.array(meshDict['triangle_tags'])  ==
    191                         numpy.array(loadedDict['triangle_tags']),
     181        self.failUnless(numpy.alltrue(numpy.array(meshDict['vertices'])  ==
     182                                      numpy.array(loadedDict['vertices'])
     183                                     ),
     184                        'test_export_mesh_file failed. Test 1')
     185        self.failUnless(numpy.alltrue(numpy.array(meshDict['triangles'])  ==
     186                                      numpy.array(loadedDict['triangles'])),
     187                        'test_export_mesh_file failed. Test 2')
     188        self.failUnless(numpy.alltrue(numpy.array(meshDict['segments'])  ==
     189                                      numpy.array(loadedDict['segments'])
     190                                     ),
     191                        'test_export_mesh_file failed. Test 3')
     192        self.failUnless(numpy.alltrue(numpy.array(meshDict['triangle_tags'])  ==
     193                                      numpy.array(loadedDict['triangle_tags'])
     194                                     ),
    192195                         'test_export_mesh_file failed. Test 4')
    193196       
    194197        self.failUnless(meshDict['vertex_attributes']  ==
    195198                        loadedDict['vertex_attributes'],
    196                          'test_export_mesh_file failed. Test 5')
    197         self.failUnless(numpy.array(meshDict['triangle_neighbors'])  ==
    198                         numpy.array(loadedDict['triangle_neighbors']),
    199                          'test_export_mesh_file failed. Test 6')
    200         self.failUnless(numpy.array(meshDict['segment_tags'])  ==
    201                         numpy.array(loadedDict['segment_tags']),
    202                          'test_export_mesh_file failed. Test 7')
    203         self.failUnless(numpy.array(meshDict['vertex_attribute_titles'])  ==
    204                         numpy.array(loadedDict['vertex_attribute_titles']),
    205                          'test_export_mesh_file failed. Test 8')
    206         self.failUnless(numpy.array(meshDict['geo_reference'])  ==
    207                         numpy.array(loadedDict['geo_reference']),
    208                          'test_export_mesh_file failed. Test 9')
     199                        'test_export_mesh_file failed. Test 5')
     200        self.failUnless(numpy.alltrue(numpy.array(meshDict['triangle_neighbors'])  ==
     201                                      numpy.array(loadedDict['triangle_neighbors'])
     202                                     ),
     203                        'test_export_mesh_file failed. Test 6')
     204        self.failUnless(numpy.alltrue(numpy.array(meshDict['segment_tags'])  ==
     205                                      numpy.array(loadedDict['segment_tags'])
     206                                     ),
     207                        'test_export_mesh_file failed. Test 7')
     208        self.failUnless(numpy.alltrue(numpy.array(meshDict['vertex_attribute_titles'])  ==
     209                                      numpy.array(loadedDict['vertex_attribute_titles'])
     210                                     ),
     211                        'test_export_mesh_file failed. Test 8')
     212        self.failUnless(numpy.alltrue(numpy.array(meshDict['geo_reference'])  ==
     213                                      numpy.array(loadedDict['geo_reference'])
     214                                     ),
     215                        'test_export_mesh_file failed. Test 9')
    209216           
    210217        os.remove(fileName)

  • anuga_core/source_numpy_conversion/anuga/mesh_engine/test_generate_mesh.py

    r5911 r5951  
    7171                              pointattlist,segattlist, mode, points)
    7272
    73         self.failUnless(data['generatedtrianglelist'] ==[(1, 0, 2), (2, 3, 1)],
    74                         'trianglelist is wrong!')
    75         self.failUnless(data['generatedsegmentlist'] ==[(0, 1), (1, 3),
    76                                                         (3, 2), (2, 0)],
    77                         'segmentlist is wrong!')
    78         self.failUnless(data['generatedpointlist'] ==[(0.0, 0.0), (0.0, 10.0),
    79                                                       (3.0, 0.0), (3.0, 10.0)],
     73        self.failUnless(numpy.alltrue(data['generatedtrianglelist'] ==
     74                                      [(1, 0, 2), (2, 3, 1)]
     75                                     ),
     76                        'trianglelist is wrong!')
     77        self.failUnless(numpy.alltrue(data['generatedsegmentlist'] ==
     78                                      [(0, 1), (1, 3), (3, 2), (2, 0)]
     79                                     ),
     80                        'segmentlist is wrong!')
     81        self.failUnless(numpy.alltrue(data['generatedpointlist'] ==
     82                                      [(0.0, 0.0), (0.0, 10.0), (3.0, 0.0), (3.0, 10.0)]
     83                                     ),
    8084                        ' is wrong!')
    8185     
    … …  
    98102                              pointattlist,segattlist, mode, points)
    99103
    100         self.failUnless(data['generatedtrianglelist'] ==[(1, 0, 2), (2, 3, 1)],
    101                         'trianglelist is wrong!')
    102         self.failUnless(data['generatedsegmentlist'] ==[(0, 1), (1, 3),
    103                                                         (3, 2), (2, 0)],
    104                         'segmentlist is wrong!')
    105         self.failUnless(data['generatedpointlist'] ==[(0.0, 0.0), (0.0, 10.0),
    106                                                       (3.0, 0.0), (3.0, 10.0)],
     104        self.failUnless(numpy.alltrue(data['generatedtrianglelist'] ==
     105                                      [(1, 0, 2), (2, 3, 1)]),
     106                        'trianglelist is wrong!')
     107        self.failUnless(numpy.alltrue(data['generatedsegmentlist'] ==
     108                                      [(0, 1), (1, 3), (3, 2), (2, 0)]),
     109                        'segmentlist is wrong!')
     110        self.failUnless(numpy.alltrue(data['generatedpointlist'] ==
     111                                      [(0.0, 0.0), (0.0, 10.0), (3.0, 0.0), (3.0, 10.0)]),
    107112                        ' is wrong!')
    108113       
    … …  
    172177                              pointattlist,segattlist, mode, points)
    173178
    174         self.failUnless(data['generatedtrianglelist'] ==[(1, 0, 2), (2, 3, 1)],
    175                         'trianglelist is wrong!')
    176         self.failUnless(data['generatedsegmentlist'] ==[(0, 1), (1, 3),
    177                                                         (3, 2), (2, 0)],
    178                         'segmentlist is wrong!')
    179         self.failUnless(data['generatedpointlist'] ==[(0.0, 0.0), (0.0, 10.0),
    180                                                       (3.0, 0.0), (3.0, 10.0)],
     179        self.failUnless(numpy.alltrue(data['generatedtrianglelist'] ==
     180                                      [(1, 0, 2), (2, 3, 1)]),
     181                        'trianglelist is wrong!')
     182        self.failUnless(numpy.alltrue(data['generatedsegmentlist'] ==
     183                                      [(0, 1), (1, 3), (3, 2), (2, 0)]),
     184                        'segmentlist is wrong!')
     185        self.failUnless(numpy.alltrue(data['generatedpointlist'] ==
     186                                      [(0.0, 0.0), (0.0, 10.0), (3.0, 0.0), (3.0, 10.0)]),
    181187                        'generatedpointlist is wrong!')
    182         self.failUnless(data['generatedsegmentmarkerlist'] ==[1,2,3,4],
     188        self.failUnless(numpy.alltrue(data['generatedsegmentmarkerlist'] == [1,2,3,4]),
    183189                        'generatedsegmentmarkerlist is wrong!')
    184190
    … …  
    351357                              pointattlist,segattlist, mode, points)
    352358
    353         self.failUnless(data['generatedtriangleattributelist'] ==[[77.0], [77.0], [77.0], [77.0]],
     359        self.failUnless(numpy.alltrue(data['generatedtriangleattributelist'] ==
     360                                      [[77.0], [77.0], [77.0], [77.0]]),
    354361                        'triangleattributelist is wrong!')
    355362
    … …  
    373380                              pointattlist,segattlist, mode, points)
    374381
    375         self.failUnless(data['generatedtriangleattributelist'] ==[[77.0], [77.0], [77.0], [77.0]],
     382        self.failUnless(numpy.alltrue(data['generatedtriangleattributelist'] ==
     383                                      [[77.0], [77.0], [77.0], [77.0]]),
    376384                        'triangleattributelist is wrong!')
    377385
    … …  
    599607        #print "data", data
    600608
    601         self.failUnless(data['generatedtrianglelist'] ==[(1, 0, 2), (2, 3, 1)],
    602                         'trianglelist is wrong!')
    603         self.failUnless(data['generatedsegmentlist'] ==[(0, 1), (1, 3),
    604                                                         (3, 2), (2, 0)],
    605                         'segmentlist is wrong!')
    606         self.failUnless(data['generatedpointlist'] ==[(0.0, 0.0), (0.0, 10.0),
    607                                                       (3.0, 0.0), (3.0, 10.0)],
    608                         ' is wrong!')
    609         self.failUnless(data['generatedtriangleattributelist'] == [[tri_tag],
    610                                                                    [tri_tag]],
    611                         ' is wrong!')
    612        
    613         self.failUnless(data['generatedsegmentlist'] == seglist,
    614                         ' is wrong!')
    615         self.failUnless(data['generatedsegmentmarkerlist'] == segattlist,
     609        self.failUnless(numpy.alltrue(data['generatedtrianglelist'] ==
     610                                      [(1, 0, 2), (2, 3, 1)]
     611                                     ),
     612                        'trianglelist is wrong!')
     613        self.failUnless(numpy.alltrue(data['generatedsegmentlist'] ==
     614                                      [(0, 1), (1, 3), (3, 2), (2, 0)]
     615                                     ),
     616                        'segmentlist is wrong!')
     617        self.failUnless(numpy.alltrue(data['generatedpointlist'] ==
     618                                      [(0.0, 0.0), (0.0, 10.0), (3.0, 0.0), (3.0, 10.0)]
     619                                     ),
     620                        ' is wrong!')
     621        self.failUnless(numpy.alltrue(data['generatedtriangleattributelist'] ==
     622                                      [[tri_tag], [tri_tag]]
     623                                     ),
     624                        ' is wrong!')
     625       
     626        self.failUnless(numpy.alltrue(data['generatedsegmentlist'] == seglist),
     627                        ' is wrong!')
     628        self.failUnless(numpy.alltrue(data['generatedsegmentmarkerlist'] == segattlist),
    616629                        ' is wrong!')
    617630        # I copied these answers from the output, so bad test..
    618         self.failUnless(data['generatedtriangleneighborlist'] == \
    619                         [(-1, 1, -1), (-1, 0, -1)],
     631        self.failUnless(numpy.alltrue(data['generatedtriangleneighborlist'] ==
     632                                      [(-1, 1, -1), (-1, 0, -1)]
     633                                     ),
    620634                        ' is wrong!')
    621635       
    … …  
    635649        data = generate_mesh(points,seglist,holelist,regionlist,
    636650                              pointattlist,segattlist, mode, points)
    637         self.failUnless(data['generatedpointattributelist'] == [[0.0],[0.0],
    638                                                                 [10],[10]],
     651        self.failUnless(numpy.alltrue(data['generatedpointattributelist'] ==
     652                                      [[0.0],[0.0],[10],[10]]
     653                                     ),
    639654                        ' is wrong!')
    640655       
    … …  
    644659        data = generate_mesh(points,seglist,holelist,regionlist,
    645660                              pointattlist,segattlist, mode, points)
    646         self.failUnless(data['generatedpointattributelist'] == [[0.0],[0.0],
    647                                                                 [10],[10]],
     661        self.failUnless(numpy.alltrue(data['generatedpointattributelist'] ==
     662                                      [[0.0],[0.0],[10],[10]]),
    648663                        ' is wrong!')
    649664        pointattlist = [[0.,1],[0.,1],[10.,20],[10.,20]]
    … …  
    652667                              pointattlist,segattlist, mode, points)
    653668        #print "data", data
    654         self.failUnless(data['generatedpointattributelist'] == pointattlist,
     669        self.failUnless(numpy.alltrue(data['generatedpointattributelist'] == pointattlist),
    655670                        ' is wrong!')
    656671           

  • anuga_core/source_numpy_conversion/anuga/pmesh/test_mesh.py

    r5912 r5951  
    426426        vert = m.getMeshVerticeAttributes()
    427427       
    428         self.failUnless(vert[0] == [12.0, 2.0] and
    429                         vert[1] == [9.0, 7.0] and
    430                         vert[2] == [14.0,3.0] and
    431                         vert[3] == [12.232233047033631, 4.4142135623730949] and
    432                         vert[4] == [13.0, 2.5] ,
     428        self.failUnless(numpy.alltrue(vert[0] == [12.0, 2.0]) and
     429                        numpy.alltrue(vert[1] == [9.0, 7.0]) and
     430                        numpy.alltrue(vert[2] == [14.0,3.0]) and
     431                        numpy.alltrue(vert[3] == [12.232233047033631, 4.4142135623730949]) and
     432                        numpy.alltrue(vert[4] == [13.0, 2.5]),
    433433                        'vertex attributes are wrong!')
    434434
    … …  
    16531653        #print "dict['vertices']",dict['vertices']
    16541654       
    1655         self.failUnless( dict['vertices'] == answer,
    1656                          'test_Mesh2IOTriangulationDict failed. test 2')
    1657 
    1658         self.failUnless( dict['vertices'] == verts,
    1659                          'test_Mesh2IOTriangulationDict failed. test vert')
    1660         self.failUnless( dict['vertex_attributes'] == vert_as,
    1661                          'test_Mesh2IOTriangulationDict failed. test vert ats')
    1662 
    1663         self.failUnless( dict['segments'][0] == [0,1],
     1655        self.failUnless(numpy.alltrue(dict['vertices'] == answer),
     1656                        'test_Mesh2IOTriangulationDict failed. test 2')
     1657
     1658        self.failUnless(numpy.alltrue(dict['vertices'] == verts),
     1659                        'test_Mesh2IOTriangulationDict failed. test vert')
     1660        self.failUnless(numpy.alltrue(dict['vertex_attributes'] == vert_as),
     1661                        'test_Mesh2IOTriangulationDict failed. test vert ats')
     1662
     1663        self.failUnless(numpy.alltrue(dict['segments'][0] == [0,1]),
    16641664                        'test_Mesh2IODict failed. test 3')
    16651665       
    1666         self.failUnless( dict['segment_tags'] == seg_tags,
     1666        self.failUnless(numpy.alltrue(dict['segment_tags'] == seg_tags),
    16671667                        'test_Mesh2IODict failed. test 3')
    16681668        #print " dict['triangles'][0]", dict['triangles'][0]
    1669         self.failUnless( dict['triangles'][0] == [3,2,4],
     1669        self.failUnless(numpy.alltrue(dict['triangles'][0] == [3,2,4]),
    16701670                        'test_Mesh2IODict failed. test 5')
    1671         self.failUnless( dict['triangle_neighbors'][0] == [-1,2,3],
     1671        self.failUnless(numpy.alltrue(dict['triangle_neighbors'][0] == [-1,2,3]),
    16721672                        'test_Mesh2IODict failed. test 6')
    16731673        #print "dict['triangle_tags'][0]", dict['triangle_tags'][0]
    … …  
    17141714        #print "dict['vertices']",dict['vertices']
    17151715       
    1716         self.failUnless( dict['vertices'] == answer,
     1716        self.failUnless(numpy.alltrue(dict['vertices'] == answer),
    17171717                         'test_Mesh2IOTriangulationDict failed. test 2')
    17181718
    1719         self.failUnless( dict['vertices'] == verts,
     1719        self.failUnless(numpy.alltrue(dict['vertices'] == verts),
    17201720                         'test_Mesh2IOTriangulationDict failed. test vert')
    1721         self.failUnless( dict['vertex_attributes'] == vert_as,
     1721        self.failUnless(numpy.alltrue(dict['vertex_attributes'] == vert_as),
    17221722                         'test_Mesh2IOTriangulationDict failed. test vert ats')
    17231723
    1724         self.failUnless( dict['segments'][0] == [0,1],
     1724        self.failUnless(numpy.alltrue(dict['segments'][0] == [0,1]),
    17251725                        'test_Mesh2IODict failed. test 3')
    17261726       
    1727         self.failUnless( dict['segment_tags'] == seg_tags,
     1727        self.failUnless(numpy.alltrue(dict['segment_tags'] == seg_tags),
    17281728                        'test_Mesh2IODict failed. test 3')
    17291729        #print " dict['triangles'][0]", dict['triangles'][0]
    1730         self.failUnless( dict['triangles'][0] == [3,2,4],
     1730        self.failUnless(numpy.alltrue(dict['triangles'][0] == [3,2,4]),
    17311731                        'test_Mesh2IODict failed. test 5')
    1732         self.failUnless( dict['triangle_neighbors'][0] == [-1,2,3],
     1732        self.failUnless(numpy.alltrue(dict['triangle_neighbors'][0] == [-1,2,3]),
    17331733                        'test_Mesh2IODict failed. test 6')
    17341734        #print "dict['triangle_tags'][0]", dict['triangle_tags'][0]
    1735         self.failUnless( dict['triangle_tags'][0] == "1.3",
    1736                          'test_Mesh2IODict failed. test 7')
     1735        self.failUnless(dict['triangle_tags'][0] == "1.3",
     1736                        'test_Mesh2IODict failed. test 7')
    17371737
    17381738        seg = m.getUserSegments()
    … …  
    21782178                        'FAILED!')
    21792179        vert= m.get_user_vertices(absolute=True)
    2180         self.failUnless(vert==points_ab,
    2181                         'FAILED!')
     2180        self.failUnless(numpy.alltrue(vert==points_ab), 'FAILED!')
    21822181
    21832182   
    … …  
    21942193                        'FAILED!')
    21952194        vert= m.get_user_vertices(absolute=True)
    2196         self.failUnless(vert==points,
    2197                         'FAILED!')
     2195        self.failUnless(numpy.alltrue(vert==points), 'FAILED!')
    21982196   
    21992197    def test_add_verticesII(self):
    … …  
    22112209                        'FAILED!')
    22122210        vert= m.get_user_vertices(absolute=True)
    2213         self.failUnless(vert==points_ab,
     2211        self.failUnless(numpy.alltrue(vert==points_ab),
    22142212                        'FAILED!')
    22152213
    … …  
    22272225                        'FAILED!')
    22282226        vert= m.get_user_vertices(absolute=True)
    2229         self.failUnless(vert==points_ab,
    2230                         'FAILED!')
     2227        self.failUnless(numpy.alltrue(vert==points_ab), 'FAILED!')
    22312228
    22322229        #geo =  Geo_reference(56,23,21)
    … …  
    22482245                        'FAILED!')
    22492246        vert= m.get_user_vertices(absolute=True)
    2250         self.failUnless(vert==points_ab,
     2247        self.failUnless(numpy.alltrue(vert==points_ab),
    22512248                        'FAILED!')
    22522249        vert= m.get_user_vertices(absolute=False)

  • anuga_core/source_numpy_conversion/anuga/shallow_water/data_manager.py

    r5947 r5951  
    602602        from Scientific.IO.NetCDF import NetCDFFile
    603603
    604         self.precision = numpy.float #Use full precision
     604        self.precision = Float #Use full precision
    605605
    606606        Data_format.__init__(self, domain, 'sww', mode)
    … …  
    10801080
    10811081    M = size  #Number of lines
    1082     xx = numpy.zeros((M,3), numpy.float)
    1083     yy = numpy.zeros((M,3), numpy.float)
    1084     zz = numpy.zeros((M,3), numpy.float)
     1082    xx = numpy.zeros((M,3), Float)
     1083    yy = numpy.zeros((M,3), Float)
     1084    zz = numpy.zeros((M,3), Float)
    10851085
    10861086    for i in range(M):
    … …  
    11261126
    11271127    M = len(lines)  #Number of lines
    1128     x = numpy.zeros((M,3), numpy.float)
    1129     y = numpy.zeros((M,3), numpy.float)
    1130     z = numpy.zeros((M,3), numpy.float)
     1128    x = numpy.zeros((M,3), Float)
     1129    y = numpy.zeros((M,3), Float)
     1130    z = numpy.zeros((M,3), Float)
    11311131
    11321132    ##i = 0
    … …  
    14431443            newcols = lenv # ncols_in_bounding_box
    14441444
    1445         telev = numpy.zeros(newcols, numpy.float)
    1446         tpoints = numpy.zeros((newcols, 2), numpy.float)
     1445        telev = numpy.zeros(newcols, Float)
     1446        tpoints = numpy.zeros((newcols, 2), Float)
    14471447
    14481448        local_index = 0
    … …  
    19931993        #the temporal dimension
    19941994        if verbose: print 'Reducing quantity %s' %quantity
    1995         q_reduced = numpy.zeros( number_of_points, numpy.float )
     1995        q_reduced = numpy.zeros( number_of_points, Float )
    19961996       
    19971997        if timestep is not None:
    … …  
    20592059    assert len(vertex_points.shape) == 2
    20602060
    2061     grid_points = numpy.zeros ( (ncols*nrows, 2), numpy.float )
     2061    grid_points = numpy.zeros ( (ncols*nrows, 2), Float )
    20622062
    20632063
    … …  
    23582358        # the temporal dimension
    23592359        if verbose: print 'Reducing quantity %s' %quantity
    2360         q_reduced = numpy.zeros( number_of_points, numpy.float )
     2360        q_reduced = numpy.zeros( number_of_points, Float )
    23612361
    23622362        for k in range(number_of_points):
    … …  
    26672667    from Scientific.IO.NetCDF import NetCDFFile
    26682668
    2669     precision = numpy.float
     2669    precision = Float
    26702670
    26712671    msg = 'Must use latitudes and longitudes for minlat, maxlon etc'
    … …  
    28212821    #Cleanup
    28222822    missing = (amplitudes == nan_ha)
    2823     if sometrue(missing):
     2823    if numpy.sometrue(missing):
    28242824        if fail_on_NaN:
    28252825            msg = 'NetCDFFile %s contains missing values'\
    … …  
    29252925    sww.store_header(outfile, times, number_of_volumes,
    29262926                     number_of_points, description=description,
    2927                      verbose=verbose,sww_precision=numpy.float)
     2927                     verbose=verbose,sww_precision=Float)
    29282928
    29292929    # Store
    29302930    from anuga.coordinate_transforms.redfearn import redfearn
    2931     x = numpy.zeros(number_of_points, numpy.float)  #Easting
    2932     y = numpy.zeros(number_of_points, numpy.float)  #Northing
     2931    x = numpy.zeros(number_of_points, Float)  #Easting
     2932    y = numpy.zeros(number_of_points, Float)  #Northing
    29332933
    29342934
    … …  
    29872987    outfile.variables['z'][:] = z             #FIXME HACK for bacwards compat.
    29882988    outfile.variables['elevation'][:] = z
    2989     outfile.variables['volumes'][:] = volumes.astype(numpy.int32) #For Opteron 64
     2989    outfile.variables['volumes'][:] = volumes.astype(Int) #For Opteron 64
    29902990
    29912991
    … …  
    31233123    d = len(q)
    31243124
    3125     T = numpy.zeros(N, numpy.float)       #Time
    3126     Q = numpy.zeros((N, d), numpy.float)  #Values
     3125    T = numpy.zeros(N, Float)       #Time
     3126    Q = numpy.zeros((N, d), Float)  #Values
    31273127
    31283128    for i, line in enumerate(lines):
    … …  
    35703570        if verbose: print 'Processing row %d of %d' %(i, nrows_new)
    35713571        lower_index = global_index
    3572         telev =  numpy.zeros(ncols_new, numpy.float)
     3572        telev =  numpy.zeros(ncols_new, Float)
    35733573        local_index = 0
    35743574        trow = i * cellsize_ratio
    … …  
    36683668    from anuga.coordinate_transforms.redfearn import redfearn
    36693669
    3670     precision = numpy.float # So if we want to change the precision its done here
     3670    precision = Float # So if we want to change the precision its done here
    36713671
    36723672    # go in to the bath dir and load the only file,
    … …  
    37953795    #Store
    37963796    from anuga.coordinate_transforms.redfearn import redfearn
    3797     x = numpy.zeros(number_of_points, numpy.float)  #Easting
    3798     y = numpy.zeros(number_of_points, numpy.float)  #Northing
     3797    x = numpy.zeros(number_of_points, Float)  #Easting
     3798    y = numpy.zeros(number_of_points, Float)  #Northing
    37993799
    38003800    if verbose: print 'Making triangular grid'
    … …  
    38593859    outfile.variables['z'][:] = z # FIXME (Ole): Remove once viewer has been recompiled and changed to use elevation instead of z
    38603860    outfile.variables['elevation'][:] = z 
    3861     outfile.variables['volumes'][:] = volumes.astype(numpy.int32) # On Opteron 64
     3861    outfile.variables['volumes'][:] = volumes.astype(Int) # On Opteron 64
    38623862
    38633863    stage = outfile.variables['stage']
    … …  
    40354035lat_name = 'LAT'
    40364036time_name = 'TIME'
    4037 precision = numpy.float # So if we want to change the precision its done here       
     4037precision = Float # So if we want to change the precision its done here       
    40384038class Write_nc:
    40394039    """
    … …  
    42604260    lonlatdep = p_array.array('f')      ###NumPy?
    42614261    lonlatdep.read(mux_file, columns * points_num)
    4262     lonlatdep = numpy.array(lonlatdep, typecode=numpy.float)   
     4262    lonlatdep = numpy.array(lonlatdep, dtype=Float)   
    42634263    lonlatdep = numpy.reshape(lonlatdep, (points_num, columns))
    42644264   
    … …  
    42674267    lon_sorted.sort()
    42684268
    4269     if not lon == lon_sorted:
     4269    if not numpy.alltrue(lon == lon_sorted):
    42704270        msg = "Longitudes in mux file are not in ascending order"
    42714271        raise IOError, msg
    … …  
    42734273    lat_sorted.sort()
    42744274
    4275     if not lat == lat_sorted:
     4275    if not numpy.alltrue(lat == lat_sorted):
    42764276        msg = "Latitudes in mux file are not in ascending order"
    42774277   
    … …  
    42874287        hz_p_array = p_array.array('f')
    42884288        hz_p_array.read(mux_file, points_num)
    4289         hz_p = numpy.array(hz_p_array, typecode=numpy.float)
     4289        hz_p = numpy.array(hz_p_array, dtype=Float)
    42904290        hz_p = numpy.reshape(hz_p, (len(lon), len(lat)))
    42914291        hz_p = numpy.transpose(hz_p) #mux has lat varying fastest, nc has long v.f.
    … …  
    43654365    QUANTITY = 2
    43664366
    4367     long_lat_dep = ensure_numeric(long_lat_dep, numpy.float)
     4367    long_lat_dep = ensure_numeric(long_lat_dep, Float)
    43684368   
    43694369    num_points = long_lat_dep.shape[0]
    … …  
    44134413    # use numeric transpose, after reshaping the quantity vector
    44144414#    quantity = zeros(len(long_lat_dep), Float)
    4415     quantity = numpy.zeros(num_points, numpy.float)
     4415    quantity = numpy.zeros(num_points, Float)
    44164416   
    44174417#    print 'num',num_points
    … …  
    47954795           
    47964796    points_utm=ensure_numeric(points_utm)
    4797     assert ensure_numeric(mesh_dic['generatedpointlist']) == \
    4798            ensure_numeric(points_utm)
     4797    assert numpy.alltrue(ensure_numeric(mesh_dic['generatedpointlist']) ==
     4798                         ensure_numeric(points_utm))
    47994799   
    48004800    volumes = mesh_dic['generatedtrianglelist']
    … …  
    48124812    sww = Write_sww()
    48134813    sww.store_header(outfile, times, len(volumes), len(points_utm),
    4814                      verbose=verbose,sww_precision=numpy.float)
     4814                     verbose=verbose,sww_precision=Float)
    48154815    outfile.mean_stage = mean_stage
    48164816    outfile.zscale = zscale
    … …  
    48354835                                 xmomentum=xmomentum,
    48364836                                 ymomentum=ymomentum,
    4837                                  sww_precision=numpy.float)
     4837                                 sww_precision=Float)
    48384838        j += 1
    48394839    if verbose: sww.verbose_quantities(outfile)
    … …  
    48714871    numSrc=len(filenames)
    48724872
    4873     file_params=-1*numpy.ones(3,numpy.float) #[nsta,dt,nt]
     4873    file_params=-1*numpy.ones(3,Float) #[nsta,dt,nt]
    48744874   
    48754875    # Convert verbose to int C flag
    … …  
    48804880       
    48814881    if permutation is None:
    4882         permutation = ensure_numeric([], numpy.float)   
     4882        permutation = ensure_numeric([], Float)   
    48834883
    48844884    # Call underlying C implementation urs2sts_ext.c   
    … …  
    49304930             
    49314931    times=dt*numpy.arange(parameters_index)   
    4932     latitudes=numpy.zeros(number_of_selected_stations, numpy.float)
    4933     longitudes=numpy.zeros(number_of_selected_stations, numpy.float)
    4934     elevation=numpy.zeros(number_of_selected_stations, numpy.float)
    4935     quantity=numpy.zeros((number_of_selected_stations, parameters_index), numpy.float)
     4932    latitudes=numpy.zeros(number_of_selected_stations, Float)
     4933    longitudes=numpy.zeros(number_of_selected_stations, Float)
     4934    elevation=numpy.zeros(number_of_selected_stations, Float)
     4935    quantity=numpy.zeros((number_of_selected_stations, parameters_index), Float)
    49364936   
    49374937   
    … …  
    50405040    if weights is None:
    50415041        # Default is equal weighting
    5042         weights=numpy.ones(numSrc,numpy.float)/numSrc
     5042        weights=numpy.ones(numSrc,Float)/numSrc
    50435043    else:
    50445044        weights = ensure_numeric(weights)
    … …  
    51745174    # 0 to number_of_points-1
    51755175    if permutation is None:
    5176         permutation = numpy.arange(number_of_points, typecode=numpy.int)
     5176        permutation = numpy.arange(number_of_points, dtype=Int)
    51775177   
    51785178   
    … …  
    51915191                     description=description,
    51925192                     verbose=verbose,
    5193                      sts_precision=numpy.float)
     5193                     sts_precision=Float)
    51945194   
    51955195    # Store
    51965196    from anuga.coordinate_transforms.redfearn import redfearn
    5197     x = numpy.zeros(number_of_points, numpy.float)  # Easting
    5198     y = numpy.zeros(number_of_points, numpy.float)  # Northing
     5197    x = numpy.zeros(number_of_points, Float)  # Easting
     5198    y = numpy.zeros(number_of_points, Float)  # Northing
    51995199
    52005200    # Check zone boundaries
    … …  
    52295229
    52305230    elevation = numpy.resize(elevation,outfile.variables['elevation'][:].shape)
    5231     outfile.variables['permutation'][:] = permutation.astype(numpy.int32) # On Opteron 64
     5231    outfile.variables['permutation'][:] = permutation.astype(Int) # On Opteron 64
    52325232    outfile.variables['x'][:] = x - geo_ref.get_xllcorner()
    52335233    outfile.variables['y'][:] = y - geo_ref.get_yllcorner()
    … …  
    55315531        outfile.variables['z'][:] = elevation
    55325532        outfile.variables['elevation'][:] = elevation  #FIXME HACK
    5533         outfile.variables['volumes'][:] = volumes.astype(numpy.int32) #On Opteron 64
     5533        outfile.variables['volumes'][:] = volumes.astype(Int) #On Opteron 64
    55345534
    55355535        q = 'elevation'
    … …  
    55395539
    55405540
    5541     def store_quantities(self, outfile, sww_precision=numpy.float32,
     5541    def store_quantities(self, outfile, sww_precision=Float32,
    55425542                         slice_index=None, time=None,
    55435543                         verbose=False, **quant):
    … …  
    57065706                     number_of_points,
    57075707                     description='Converted from URS mux2 format',
    5708                      sts_precision=numpy.float32,
     5708                     sts_precision=Float32,
    57095709                     verbose=False):
    57105710        """
    … …  
    57655765
    57665766        # Doing sts_precision instead of Float gives cast errors.
    5767         outfile.createVariable('time', Floatr, ('number_of_timesteps',))
     5767        outfile.createVariable('time', Float, ('number_of_timesteps',))
    57685768
    57695769        for q in Write_sts.sts_quantities:
    … …  
    58745874        outfile.variables[q+Write_sts.RANGE][1] = max(elevation)
    58755875
    5876     def store_quantities(self, outfile, sts_precision=numpy.float32,
     5876    def store_quantities(self, outfile, sts_precision=Float32,
    58775877                         slice_index=None, time=None,
    58785878                         verbose=False, **quant):
    … …  
    59645964        lonlatdep = p_array.array('f')
    59655965        lonlatdep.read(mux_file, columns * self.points_num)
    5966         lonlatdep = numpy.array(lonlatdep, typecode=numpy.float)   
     5966        lonlatdep = numpy.array(lonlatdep, dtype=Float)   
    59675967        lonlatdep = numpy.reshape(lonlatdep, (self.points_num, columns))
    59685968        #print 'lonlatdep',lonlatdep
    … …  
    59935993        hz_p_array = p_array.array('f')
    59945994        hz_p_array.read(self.mux_file, self.points_num)
    5995         hz_p = numpy.array(hz_p_array, typecode=numpy.float)
     5995        hz_p = numpy.array(hz_p_array, dtype=Float)
    59965996        self.iter_time_step += 1
    59975997       
    … …  
    61086108    #array to store data, number in there is to allow float...
    61096109    #i'm sure there is a better way!
    6110     data=numpy.array([],typecode=numpy.float)
     6110    data=numpy.array([],dtype=Float)
    61116111    data=numpy.resize(data,((len(lines)-1),len(header_fields)))
    61126112#    print 'number of fields',range(len(header_fields))
    … …  
    66856685   
    66866686        # Get the relevant quantities (Convert from single precison)
    6687         elevation = numpy.array(fid.variables['elevation'][:], numpy.float)
    6688         stage = numpy.array(fid.variables['stage'][:], numpy.float)
     6687        elevation = numpy.array(fid.variables['elevation'][:], Float)
     6688        stage = numpy.array(fid.variables['stage'][:], Float)
    66896689   
    66906690   

  • anuga_core/source_numpy_conversion/anuga/shallow_water/test_data_manager.py

    r5947 r5951  
    49644964                         'failed')
    49654965
    4966         self.failUnless(m2d == [[5,6],[9,10]],
    4967                          'failed')
     4966        self.failUnless(numpy.alltrue(m2d == [[5,6],[9,10]]),
     4967                        'failed')
    49684968
    49694969    def test_get_min_max_indexes_lat_ascending(self):
    … …  
    50305030        self.failUnless(latitudes_new == [-30, -35, -40] and \
    50315031                        longitudes_new == [148, 149,150],
    5032                          'failed')
    5033         self.failUnless(m2d == [[0,1,2],[4,5,6],[8,9,10]],
    5034                          'failed')
     5032                        'failed')
     5033        self.failUnless(numpy.alltrue(m2d == [[0,1,2],[4,5,6],[8,9,10]]),
     5034                        'failed')
    50355035
    50365036    def test_get_min_max_indexes3(self):

  • anuga_core/source_numpy_conversion/anuga/shallow_water/tsunami_okada.py

    r5901 r5951  
    3131
    3232import numpy
     33from anuga.config import Float
     34
    3335
    3436def earthquake_tsunami(ns,NSMAX,length, width, strike, depth,\
    … …  
    4446    zrec=zrec0.get_vertex_values(xy=True)
    4547   
    46     x0= numpy.zeros(ns,numpy.float)
    47     y0= numpy.zeros(ns,numpy.float)
     48    x0= numpy.zeros(ns,Float)
     49    y0= numpy.zeros(ns,Float)
    4850    if ns ==1:
    4951        x0[0]=xi
    … …  
    150152        zrec=self.zrec
    151153        #initialization
    152         disp0=numpy.zeros(3,numpy.float)
    153         strain0=numpy.zeros(6,numpy.float)
    154         tilt0 = numpy.zeros(2,numpy.float)
    155         dislocations=numpy.zeros(ns,numpy.float)
    156         depths=numpy.zeros(ns,numpy.float)
    157         strikes= numpy.zeros(ns,numpy.float)
    158         lengths= numpy.zeros(ns,numpy.float)
    159         slips= numpy.zeros(ns,numpy.float)
    160         rakes= numpy.zeros(ns,numpy.float)
    161         widths= numpy.zeros(ns,numpy.float)
    162         dips= numpy.zeros(ns,numpy.float)
    163         strikes= numpy.zeros(ns,numpy.float)
    164         strikes= numpy.zeros(ns,numpy.float)
    165         strain = numpy.zeros((N,6),numpy.float)
    166         disp = numpy.zeros((N,3),numpy.float)
    167         tilt = numpy.zeros((N,2),numpy.float)
    168         xs =numpy.zeros(ns,numpy.float)
    169         ys =numpy.zeros(ns,numpy.float)
     154        disp0=numpy.zeros(3,Float)
     155        strain0=numpy.zeros(6,Float)
     156        tilt0 = numpy.zeros(2,Float)
     157        dislocations=numpy.zeros(ns,Float)
     158        depths=numpy.zeros(ns,Float)
     159        strikes= numpy.zeros(ns,Float)
     160        lengths= numpy.zeros(ns,Float)
     161        slips= numpy.zeros(ns,Float)
     162        rakes= numpy.zeros(ns,Float)
     163        widths= numpy.zeros(ns,Float)
     164        dips= numpy.zeros(ns,Float)
     165        strikes= numpy.zeros(ns,Float)
     166        strikes= numpy.zeros(ns,Float)
     167        strain = numpy.zeros((N,6),Float)
     168        disp = numpy.zeros((N,3),Float)
     169        tilt = numpy.zeros((N,2),Float)
     170        xs =numpy.zeros(ns,Float)
     171        ys =numpy.zeros(ns,Float)
    170172        z=[]
    171173        if ns==1:
    … …  
    213215            xrec=y
    214216            yrec=x
     217##            print 'xrec=%s' % str(xrec)
     218##            print 'yrec=%s' % str(yrec)
     219##            print 'zrec=%s' % str(zrec)
    215220            for i in range(0,len(zrec[0])):
    216                 if zrec[0][i]==yrec and zrec[1][i]==xrec:
     221                if numpy.alltrue(zrec[0][i]==yrec) and \
     222                   numpy.alltrue(zrec[1][i]==xrec):
     223##                    print 'Setting Z!'
    217224                    Z=zrec[2][i]
    218225                    Z=0.001*Z
    … …  
    377384         
    378385          F0=0.0                                                   
    379           U=numpy.zeros((12,1),numpy.float)
    380           DUA=numpy.zeros((12,1),numpy.float)
    381           DUB=numpy.zeros((12,1),numpy.float)
    382           DUC=numpy.zeros((12,1),numpy.float)
     386          U=numpy.zeros((12,1),Float)
     387          DUA=numpy.zeros((12,1),Float)
     388          DUB=numpy.zeros((12,1),Float)
     389          DUC=numpy.zeros((12,1),Float)
    383390         
    384391         
    … …  
    523530       WZ=self.WZ
    524531
    525        DUA=numpy.zeros((12,1),numpy.float)
    526        DU=numpy.zeros((12,1),numpy.float)
    527        U=numpy.zeros((12,1),numpy.float)
     532       DUA=numpy.zeros((12,1),Float)
     533       DU=numpy.zeros((12,1),Float)
     534       U=numpy.zeros((12,1),Float)
    528535#-----                                                               
    529536       for I in range(0,12):                                                   
    … …  
    633640#      DATA PI2/6.283185307179586D0/
    634641       
    635         DUB=numpy.zeros((12,1),numpy.float)
    636         DU=numpy.zeros((12,1),numpy.float)
    637         U=numpy.zeros((12,1),numpy.float)
     642        DUB=numpy.zeros((12,1),Float)
     643        DU=numpy.zeros((12,1),Float)
     644        U=numpy.zeros((12,1),Float)
    638645       
    639646        F0=0.0
    … …  
    807814#      DATA PI2/6.283185307179586D0/
    808815       
    809       DUC=numpy.zeros((12,1),numpy.float)
    810       DU=numpy.zeros((12,1),numpy.float)
    811       U=numpy.zeros((12,1),numpy.float)
     816      DUC=numpy.zeros((12,1),Float)
     817      DU=numpy.zeros((12,1),Float)
     818      U=numpy.zeros((12,1),Float)
    812819             
    813820      F0=0.0
    … …  
    10011008      EPS=0.000001
    10021009       
    1003       XI=numpy.zeros(2,numpy.float)
    1004       ET=numpy.zeros(2,numpy.float)
    1005       KXI=numpy.zeros(2,numpy.float)
    1006       KET=numpy.zeros(2,numpy.float)
    1007       U=numpy.zeros(12,numpy.float)
    1008       DU=numpy.zeros(12,numpy.float)
    1009       DUA=numpy.zeros(12,numpy.float)
    1010       DUB=numpy.zeros(12,numpy.float)
    1011       DUC=numpy.zeros(12,numpy.float)
     1010      XI=numpy.zeros(2,Float)
     1011      ET=numpy.zeros(2,Float)
     1012      KXI=numpy.zeros(2,Float)
     1013      KET=numpy.zeros(2,Float)
     1014      U=numpy.zeros(12,Float)
     1015      DU=numpy.zeros(12,Float)
     1016      DUA=numpy.zeros(12,Float)
     1017      DUB=numpy.zeros(12,Float)
     1018      DUC=numpy.zeros(12,Float)
    10121019
    10131020#-----                                                                 
    … …  
    12061213#      DATA F0,F2,PI2/0.D0,2.D0,6.283185307179586D0/
    12071214
    1208       U=numpy.zeros(12,numpy.float)
    1209       DU=numpy.zeros(12,numpy.float)
    1210       DUA=numpy.zeros(12,numpy.float)
     1215      U=numpy.zeros(12,Float)
     1216      DU=numpy.zeros(12,Float)
     1217      DUA=numpy.zeros(12,Float)
    12111218      F0 =0.0
    12121219      F2=2.0
    … …  
    13371344#      DATA  F0,F1,F2,PI2/0.D0,1.D0,2.D0,6.283185307179586D0/
    13381345
    1339       DUB=numpy.zeros(12,numpy.float)
    1340       DU=numpy.zeros(12,numpy.float)
    1341       U=numpy.zeros(12,numpy.float)
     1346      DUB=numpy.zeros(12,Float)
     1347      DU=numpy.zeros(12,Float)
     1348      U=numpy.zeros(12,Float)
    13421349     
    13431350      F0=0.0
    … …  
    15021509 #     DATA F0,F1,F2,F3,PI2/0.D0,1.D0,2.D0,3.D0,6.283185307179586D0/
    15031510
    1504       DUC=numpy.zeros(12,numpy.float)
    1505       DU=numpy.zeros(12,numpy.float)
    1506       U=numpy.zeros(12,numpy.float)
     1511      DUC=numpy.zeros(12,Float)
     1512      DU=numpy.zeros(12,Float)
     1513      U=numpy.zeros(12,Float)
    15071514     
    15081515      F0=0.0
Note: See TracChangeset for help on using the changeset viewer.