Ignore:
Timestamp:
Aug 24, 2005, 1:49:47 PM (19 years ago)
Author:
ole
Message:

Embedded caching functionality within quantity.set_values and modified validation example lwru2.py to illustrate the advantages that can be gained from supervised caching.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/pyvolution/test_shallow_water.py

    r1671 r1753  
    12381238        val3 = 2.
    12391239
    1240         domain.set_quantity('stage', [val0, val1, val2, val3], 'centroids')
     1240        domain.set_quantity('stage', [val0, val1, val2, val3],
     1241                            location='centroids')
    12411242        L = domain.quantities['stage'].vertex_values
    12421243
     
    12741275            return x**2
    12751276
    1276         domain.set_quantity('stage', stage, 'centroids')
     1277        domain.set_quantity('stage', stage, location='centroids')
    12771278
    12781279        a, b = domain.quantities['stage'].compute_gradients()
     
    13111312            return x**4+y**2
    13121313
    1313         domain.set_quantity('stage', stage, 'centroids')
     1314        domain.set_quantity('stage', stage, location='centroids')
    13141315        #print domain.quantities['stage'].centroid_values
    13151316
     
    13551356
    13561357        domain.set_quantity('elevation', slope)
    1357         domain.set_quantity('stage', stage, 'centroids')
     1358        domain.set_quantity('stage', stage, location='centroids')
    13581359
    13591360        #print domain.quantities['elevation'].centroid_values
     
    14431444        domain.set_quantity('elevation', slope)
    14441445        domain.set_quantity('stage',
    1445                             [0.01298164, 0.00365611, 0.01440365, -0.0381856437096],
    1446                             'centroids')
     1446                            [0.01298164, 0.00365611,
     1447                             0.01440365, -0.0381856437096],
     1448                            location='centroids')
    14471449        domain.set_quantity('xmomentum',
    1448                             [0.00670439, 0.01263789, 0.00647805, 0.0178180740668],
    1449                             'centroids')
     1450                            [0.00670439, 0.01263789,
     1451                             0.00647805, 0.0178180740668],
     1452                            location='centroids')
    14501453        domain.set_quantity('ymomentum',
    1451                             [-7.23510980e-004, -6.30413883e-005, 6.30413883e-005, 0.000200907255866],
    1452                             'centroids')
     1454                            [-7.23510980e-004, -6.30413883e-005,
     1455                             6.30413883e-005, 0.000200907255866],
     1456                            location='centroids')
    14531457
    14541458        E = domain.quantities['elevation'].vertex_values
     
    14661470        #print Y[1,:]
    14671471
    1468         assert allclose(L[1,:], [-0.00825735775384, -0.00801881482869, 0.0272445025825])
    1469         assert allclose(X[1,:], [0.0143507718962, 0.0142502147066, 0.00931268339717])
    1470         assert allclose(Y[1,:], [-0.000117062180693, 7.94434448109e-005, -0.000151505429018])
     1472        assert allclose(L[1,:], [-0.00825735775384,
     1473                                 -0.00801881482869,
     1474                                 0.0272445025825])
     1475        assert allclose(X[1,:], [0.0143507718962,
     1476                                 0.0142502147066,
     1477                                 0.00931268339717])
     1478        assert allclose(Y[1,:], [-0.000117062180693,
     1479                                 7.94434448109e-005,
     1480                                 -0.000151505429018])
    14711481
    14721482
     
    22112221
    22122222
    2213                 domain.set_quantity('stage', L, 'centroids')
    2214                 domain.set_quantity('xmomentum', X, 'centroids')
    2215                 domain.set_quantity('ymomentum', Y, 'centroids')
     2223                domain.set_quantity('stage', L, location='centroids')
     2224                domain.set_quantity('xmomentum', X, location='centroids')
     2225                domain.set_quantity('ymomentum', Y, location='centroids')
    22162226
    22172227                domain.check_integrity()
Note: See TracChangeset for help on using the changeset viewer.