Ignore:
Timestamp:
Jan 13, 2009, 5:42:30 PM (15 years ago)
Author:
rwilson
Message:

Change Numeric imports to general form - ready to change to NumPy?.

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r5252 r6157  
    3030"""
    3131
     32import Numeric as num
     33
     34
    3235def earthquake_tsunami(ns,NSMAX,length, width, strike, depth,\
    3336                       dip, xi, yi,z0, slip, rake,\
    … …  
    3639    from anuga.abstract_2d_finite_volumes.quantity import Quantity
    3740    from math import sin, radians
    38     from Numeric import zeros, Float
    3941
    4042    #zrec0 = Quantity(domain)
    … …  
    4345    zrec=zrec0.get_vertex_values(xy=True)
    4446   
    45     x0= zeros(ns,Float)
    46     y0= zeros(ns,Float)
     47    x0= num.zeros(ns,num.Float)
     48    y0= num.zeros(ns,num.Float)
    4749    if ns ==1:
    4850        x0[0]=xi
    … …  
    130132        produced by a submarine mass failure.
    131133        """
    132         from Numeric import zeros, Float
    133134        from string import replace,strip
    134135        from math import sin, cos, radians, exp, cosh
    135         from Numeric import zeros, Float
    136136        #ensure vectors x and y have the same length
    137137        N = len(x)
    … …  
    151151        zrec=self.zrec
    152152        #initialization
    153         disp0=zeros(3,Float)
    154         strain0=zeros(6,Float)
    155         tilt0 = zeros(2,Float)
    156         dislocations=zeros(ns,Float)
    157         depths=zeros(ns,Float)
    158         strikes= zeros(ns,Float)
    159         lengths= zeros(ns,Float)
    160         slips= zeros(ns,Float)
    161         rakes= zeros(ns,Float)
    162         widths= zeros(ns,Float)
    163         dips= zeros(ns,Float)
    164         strikes= zeros(ns,Float)
    165         strikes= zeros(ns,Float)
    166         strain = zeros((N,6),Float)
    167         disp = zeros((N,3),Float)
    168         tilt = zeros((N,2),Float)
    169         xs =zeros(ns,Float)
    170         ys =zeros(ns,Float)
     153        disp0=num.zeros(3,num.Float)
     154        strain0=num.zeros(6,num.Float)
     155        tilt0 = num.zeros(2,num.Float)
     156        dislocations=num.zeros(ns,num.Float)
     157        depths=num.zeros(ns,num.Float)
     158        strikes= num.zeros(ns,num.Float)
     159        lengths= num.zeros(ns,num.Float)
     160        slips= num.zeros(ns,num.Float)
     161        rakes= num.zeros(ns,num.Float)
     162        widths= num.zeros(ns,num.Float)
     163        dips= num.zeros(ns,num.Float)
     164        strikes= num.zeros(ns,num.Float)
     165        strikes= num.zeros(ns,num.Float)
     166        strain = num.zeros((N,6),num.Float)
     167        disp = num.zeros((N,3),num.Float)
     168        tilt = num.zeros((N,2),num.Float)
     169        xs =num.zeros(ns,num.Float)
     170        ys =num.zeros(ns,num.Float)
    171171        z=[]
    172172        if ns==1:
    … …  
    376376#        COMMON /C1/DUMMY(8),R,dumm(15)                                   
    377377#          DIMENSION  U(12),DUA(12),DUB(12),DUC(12)
    378           from Numeric import zeros, Float
    379378         
    380379          F0=0.0                                                   
    381           U=zeros((12,1),Float)
    382           DUA=zeros((12,1),Float)
    383           DUB=zeros((12,1),Float)
    384           DUC=zeros((12,1),Float)
     380          U=num.zeros((12,1),num.Float)
     381          DUA=num.zeros((12,1),num.Float)
     382          DUB=num.zeros((12,1),num.Float)
     383          DUC=num.zeros((12,1),num.Float)
    385384         
    386385         
    … …  
    483482#    *           UY,VY,WY,UZ,VZ,WZ
    484483
    485        from Numeric import zeros, Float
    486484         
    487485                                                             
    … …  
    528526       WZ=self.WZ
    529527
    530        DUA=zeros((12,1),Float)
    531        DU=zeros((12,1),Float)
    532        U=zeros((12,1),Float)
     528       DUA=num.zeros((12,1),num.Float)
     529       DU=num.zeros((12,1),num.Float)
     530       U=num.zeros((12,1),num.Float)
    533531#-----                                                               
    534532       for I in range(0,12):                                                   
    … …  
    637635#     *        /0.D0,1.D0,2.D0,3.D0,4.D0,5.D0,8.D0,9.D0/                 
    638636#      DATA PI2/6.283185307179586D0/
    639         from Numeric import zeros, Float
    640637       
    641         DUB=zeros((12,1),Float)
    642         DU=zeros((12,1),Float)
    643         U=zeros((12,1),Float)
     638        DUB=num.zeros((12,1),num.Float)
     639        DU=num.zeros((12,1),num.Float)
     640        U=num.zeros((12,1),num.Float)
    644641       
    645642        F0=0.0
    … …  
    814811       
    815812
    816       from Numeric import zeros, Float
    817        
    818       DUC=zeros((12,1),Float)
    819       DU=zeros((12,1),Float)
    820       U=zeros((12,1),Float)
     813      DUC=num.zeros((12,1),num.Float)
     814      DU=num.zeros((12,1),num.Float)
     815      U=num.zeros((12,1),num.Float)
    821816             
    822817      F0=0.0
    … …  
    10061001 #     DIMENSION  XI(2),ET(2),KXI(2),KET(2)                             
    10071002 #     DIMENSION  U(12),DU(12),DUA(12),DUB(12),DUC(12)                 
    1008       from Numeric import zeros, Float
    10091003      from math import sqrt
    10101004      F0 =0.0
    10111005      EPS=0.000001
    10121006       
    1013       XI=zeros(2,Float)
    1014       ET=zeros(2,Float)
    1015       KXI=zeros(2,Float)
    1016       KET=zeros(2,Float)
    1017       U=zeros(12,Float)
    1018       DU=zeros(12,Float)
    1019       DUA=zeros(12,Float)
    1020       DUB=zeros(12,Float)
    1021       DUC=zeros(12,Float)
     1007      XI=num.zeros(2,num.Float)
     1008      ET=num.zeros(2,num.Float)
     1009      KXI=num.zeros(2,num.Float)
     1010      KET=num.zeros(2,num.Float)
     1011      U=num.zeros(12,num.Float)
     1012      DU=num.zeros(12,num.Float)
     1013      DUA=num.zeros(12,num.Float)
     1014      DUB=num.zeros(12,num.Float)
     1015      DUC=num.zeros(12,num.Float)
    10221016
    10231017#-----                                                                 
    … …  
    12171211
    12181212
    1219       from Numeric import zeros, Float
    1220        
    1221       U=zeros(12,Float)
    1222       DU=zeros(12,Float)
    1223       DUA=zeros(12,Float)
     1213      U=num.zeros(12,num.Float)
     1214      DU=num.zeros(12,num.Float)
     1215      DUA=num.zeros(12,num.Float)
    12241216      F0 =0.0
    12251217      F2=2.0
    … …  
    13501342#      DATA  F0,F1,F2,PI2/0.D0,1.D0,2.D0,6.283185307179586D0/
    13511343
    1352       from Numeric import zeros, Float
    1353        
    1354       DUB=zeros(12,Float)
    1355       DU=zeros(12,Float)
    1356       U=zeros(12,Float)
     1344      DUB=num.zeros(12,num.Float)
     1345      DU=num.zeros(12,num.Float)
     1346      U=num.zeros(12,num.Float)
    13571347     
    13581348      F0=0.0
    … …  
    15171507 #     DATA F0,F1,F2,F3,PI2/0.D0,1.D0,2.D0,3.D0,6.283185307179586D0/
    15181508
    1519       from Numeric import zeros, Float
    1520        
    1521       DUC=zeros(12,Float)
    1522       DU=zeros(12,Float)
    1523       U=zeros(12,Float)
     1509      DUC=num.zeros(12,num.Float)
     1510      DU=num.zeros(12,num.Float)
     1511      U=num.zeros(12,num.Float)
    15241512     
    15251513      F0=0.0
Note: See TracChangeset for help on using the changeset viewer.