Changeset 3757


Ignore:
Timestamp:
Oct 11, 2006, 5:05:39 PM (17 years ago)
Author:
ole
Message:

Work towards automated validation of ANUGA using okushiri

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/automated_validation_tests/okushiri_tank_validation/validate_okushiri.py

    r3722 r3757  
    2222import project
    2323
     24try:
     25    import pylab
     26except:
     27    plotting = False
     28else:
     29    plotting = True
     30   
     31
    2432finaltime = 22.5
    2533timestep = 0.05
     
    2836# Create Domain
    2937#-------------------------
    30 #N = 50
    31 N = 10
     38N = 50
     39#N = 10
    3240points, vertices, boundary = rectangular_cross(N, N/5*3,
    3341                                               len1=5.448, len2=3.402)
     
    110118
    111119
    112 predicted_gauge_values = [[], [],[],[]]
     120predicted_gauge_values = [[], [], [], []]
    113121
    114122
     
    126134    stage = domain.get_quantity('stage')
    127135    for i, (x, y) in enumerate(gauge_locations):
    128         w = stage.get_values(interpolation_points=[[x,y]])
     136        w = stage.get_values(interpolation_points=[[x,y]])[0]
    129137        predicted_gauge_values[i].append(w)
    130138
    131         print 'difference', x, y, w - reference_gauge_values[i][j]
     139        #print 'difference', x, y,\
     140        #      predicted_gauge_values[i][j] - reference_gauge_values[i][j]
    132141   
    133142
     
    138147    predicted_gauge_values[i] = ensure_numeric(predicted_gauge_values[i])
    139148
     149    print predicted_gauge_values[i].shape
     150    print reference_gauge_values[i].shape
     151    #print predicted_gauge_values.shape
     152    #print reference_gauge_values.shape   
    140153   
    141     sqsum = sum((predicted_gauge_values[i] - reference_gauge_values)**2)
    142     denom = sum(reference_gauge_values**2)
     154    sqsum = sum((predicted_gauge_values[i].flat - \
     155                 reference_gauge_values[i].flat)**2)
     156    denom = sum(reference_gauge_values[i].flat**2)
    143157
    144158    print i, sqsum, sqsum/denom
    145159
    146     msg 'Rms error is %f' %(sqsum/denom)
     160    msg = 'Rms error is %f' %(sqsum/denom)
     161    print msg
     162   
    147163    assert sqsum/denom < 0.01, msg
     164   
    148165
    149     try:
     166    if plotting is True:
     167       
    150168        from pylab import *
    151     except:
    152         pass
    153     else:
     169
    154170        ion()
    155171        hold(False)
     
    165181       
    166182
    167 try:
    168     from pylab import *
    169 except:
    170     pass
    171 else:
     183if plotting is True:
    172184    show()
    173185   
Note: See TracChangeset for help on using the changeset viewer.