Changeset 7509


Ignore:
Timestamp:
Sep 11, 2009, 4:04:14 PM (11 years ago)
Author:
ole
Message:

Fiddled with ticket:223 and discovered problem with Nans and 'filler'
Test is disabled again.

Location:
anuga_core/source/anuga/shallow_water
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/data_manager.py

    r7346 r7509  
    35603560    give a different final boundary, or crash.
    35613561    """
    3562 
     3562   
    35633563    from Scientific.IO.NetCDF import NetCDFFile
    35643564    from shallow_water import Domain
     
    36443644        coordinates, volumes, boundary = weed(coordinates, volumes,boundary)
    36453645
     3646     
     3647   
    36463648    try:
    36473649        domain = Domain(coordinates, volumes, boundary)
  • anuga_core/source/anuga/shallow_water/test_data_manager.py

    r7342 r7509  
    37163716        ##################################################################
    37173717
     3718        #FIXME: See ticket 223
    37183719
    37193720        from mesh_factory import rectangular
     
    37623763        filename = domain.datadir + os.sep + domain.get_name() + '.sww'
    37633764
    3764         #Fail because NaNs are present
     3765        # Fail because NaNs are present
     3766        #domain2 = sww2domain(filename,
     3767        #                     boundary,
     3768        #                     fail_if_NaN=True,
     3769        #                     verbose=self.verbose)       
    37653770        try:
    3766             domain2 = sww2domain(filename,boundary,fail_if_NaN=True,verbose=self.verbose)
    3767         except:
    3768             #Now import it, filling NaNs to be 0
    3769             filler = 0
    3770             domain2 = sww2domain(filename,None,fail_if_NaN=False,NaN_filler = filler,verbose=self.verbose)
    3771 
    3772         #Clean up
     3771            domain2 = sww2domain(filename,
     3772                                 boundary,
     3773                                 fail_if_NaN=True,
     3774                                 verbose=self.verbose)
     3775        except DataDomainError:
     3776            # Now import it, filling NaNs to be -9999
     3777            filler = -9999
     3778            domain2 = sww2domain(filename,
     3779                                 None,
     3780                                 fail_if_NaN=False,
     3781                                 NaN_filler=filler,
     3782                                 verbose=self.verbose)
     3783        else:
     3784            raise Exception, 'should have failed'
     3785
     3786           
     3787        # Now import it, filling NaNs to be 0
     3788        filler = -9999
     3789        domain2 = sww2domain(filename,
     3790                             None,
     3791                             fail_if_NaN=False,
     3792                             NaN_filler=filler,
     3793                             verbose=self.verbose)           
     3794                             
     3795        import sys; sys.exit()
     3796           
     3797        # Clean up
    37733798        os.remove(filename)
    37743799
     
    37863811            assert num.allclose(eval('domain.'+bit),eval('domain2.'+bit))
    37873812
    3788         #print filler
    3789         #print max(max(domain2.get_quantity('xmomentum').get_values()))
     3813        print
     3814        print
     3815        print domain2.get_quantity('xmomentum').get_values()
     3816        print
     3817        print domain2.get_quantity('stage').get_values()
     3818        print
     3819             
     3820        print 'filler', filler
     3821        print max(domain2.get_quantity('xmomentum').get_values().flat)
    37903822       
    37913823        assert max(max(domain2.get_quantity('xmomentum').get_values()))==filler
     
    1167911711
    1168011712if __name__ == "__main__":
     11713    #suite = unittest.makeSuite(Test_Data_Manager, 'test_sww2domain2')
    1168111714    suite = unittest.makeSuite(Test_Data_Manager, 'test')
     11715   
    1168211716   
    1168311717   
Note: See TracChangeset for help on using the changeset viewer.