acoular 23.11__py3-none-any.whl → 24.5__py3-none-any.whl

acoular/tests/run_tests.sh

@@ -1,18 +0,0 @@
-#!/bin/bash
-
-#remove cache data before testing
-rm -rf ./cache/*
-rm -rf ./__pycache__
-
-#build a test suite object which runs the tests in this folder
-python -m unittest discover -v -p "test_*.py"
-
-VAL=$?
-
-#remove cache data after testing
-rm -rf ./cache/*
-rm -rf ./__pycache__
-
-exit $VAL
-
-

acoular/tests/run_tests_osx.sh

@@ -1,16 +0,0 @@
-#!/bin/bash
-
-#remove cache data before testing
-rm -rf ./cache/*
-rm -rf ./__pycache__
-
-#build a test suite object which runs the tests in this folder
-python3 -m unittest discover -v -p "test_*.py"
-
-VAL=$?
-
-#remove cache data after testing
-rm -rf ./cache/*
-rm -rf ./__pycache__
-
-exit $VAL

acoular/tests/test_beamformer_results.py

@@ -1,204 +0,0 @@
-# -*- coding: utf-8 -*-
-#pylint: disable-msg=E0611, E1101, C0103, R0901, R0902, R0903, R0904, W0232
-#------------------------------------------------------------------------------
-# Copyright (c) Acoular Development Team.
-#------------------------------------------------------------------------------
-"""Implements testing of frequency beamformers.
-"""
-
-import unittest
-
-from os.path import join
-
-import numpy as np
-
-#acoular imports
-import acoular
-acoular.config.global_caching = 'none' # to make sure that nothing is cached
-
-from acoular import Calib, MicGeom, PowerSpectra, \
-RectGrid, BeamformerBase, BeamformerEig, BeamformerOrth, BeamformerCleansc, \
-MaskedTimeSamples, BeamformerCMF, \
-BeamformerCapon, BeamformerMusic, BeamformerDamas, BeamformerClean, \
-BeamformerFunctional, BeamformerDamasPlus, BeamformerGIB, BeamformerGridlessOrth,\
-SteeringVector, Environment
-
-# if this flag is set to True
-WRITE_NEW_REFERENCE_DATA = False
-# new beamformer results are generated for comparison during testing. 
-# Should always be False. Only set to True if it is necessary to 
-# recalculate the data due to intended changes of the Beamformers.
-
-
-#load exampledata
-datafile = join('..','..','examples','example_data.h5')
-calibfile = join('..','..','examples','example_calib.xml')
-micgeofile = join( '..','xml','array_56.xml')
-
-#frequencies to test
-cfreqs = 1000,8000
-
-# values from example 1
-t1 = MaskedTimeSamples(name=datafile)
-t1.start = 0 # first sample, default
-t1.stop = 16000 # last valid sample = 15999
-invalid = [1,7] # list of invalid channels (unwanted microphones etc.)
-t1.invalid_channels = invalid 
-t1.calib = Calib(from_file=calibfile)
-m = MicGeom(from_file=micgeofile)
-m.invalid_channels = invalid
-g = RectGrid(x_min=-0.6, x_max=-0.0, y_min=-0.3, y_max=0.3, z=0.68,
-             increment=0.05 )
-env=Environment(c=346.04)
-st = SteeringVector(grid=g, mics=m, env=env)
-f = PowerSpectra(time_data=t1, 
-               window='Hanning', overlap='50%', block_size=128, #FFT-parameters
-               cached = False )  
-
-# produces a tuple of beamformer objects to test
-# because we need new objects for each test we have to call this more than once
-def fbeamformers():
-    bb = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = False)
-    be = BeamformerEig(freq_data=f, steer=st, r_diag=True, n=54, cached = False)
-
-    #frequency beamformers to test
-    bbase = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = False)
-    bc = BeamformerCapon(freq_data=f, steer=st, cached=False)
-    beig = BeamformerEig(freq_data=f, steer=st, r_diag=True, n=54, cached = False)
-    bm = BeamformerMusic(freq_data=f, steer=st, n=6, cached = False)
-    bd = BeamformerDamas(beamformer=bb, n_iter=10, cached = False)
-    bdp = BeamformerDamasPlus(beamformer=bb, n_iter=100, cached = False)
-    bo = BeamformerOrth(freq_data=f, steer=st, r_diag=True, eva_list=list(range(38,54)), cached = False)
-    bs = BeamformerCleansc(freq_data=f, steer=st, r_diag=True, cached = False)
-    bcmf = BeamformerCMF(freq_data=f, steer=st, method='LassoLarsBIC', cached = False)
-    bl = BeamformerClean(beamformer=bb, n_iter=10, cached = False)
-    bf = BeamformerFunctional(freq_data=f, steer=st, r_diag=False, gamma=3, cached = False)
-    bgib = BeamformerGIB(freq_data=f, steer=st, method= 'LassoLars', n=2, cached = False)
-    bgo = BeamformerGridlessOrth(freq_data=f, steer=st, r_diag=False, n=1,  shgo={'n':16}, cached = False)
-    return (bbase, bc, beig, bm, bl, bo, bs, bd, bcmf, bf, bdp, bgib, bgo)
-
-class acoular_beamformer_test(unittest.TestCase):
-
-    def test_beamformer_freq_results(self):
-        # we expect the results to computed
-        acoular.config.global_caching = 'none'
-        for b in fbeamformers():
-            with self.subTest(b.__class__.__name__+" global_caching = none"):
-                name = join('reference_data',f'{b.__class__.__name__}.npy')
-                # stack all frequency band results together
-                actual_data = np.array([b.synthetic(cf,1) for cf in cfreqs],dtype=np.float32)
-                if WRITE_NEW_REFERENCE_DATA:
-                    np.save(name,actual_data)
-                ref_data = np.load(name)
-                np.testing.assert_allclose(actual_data, ref_data, rtol=5e-5, atol=5e-8)
-        # we expect the results to be computed and written to cache
-        acoular.config.global_caching = 'individual'
-        for b in fbeamformers():
-            b.cached = True
-            with self.subTest(b.__class__.__name__+" global_caching = individual"):
-                name = join('reference_data',f'{b.__class__.__name__}.npy')
-                actual_data = np.array([b.synthetic(cf,1) for cf in cfreqs],dtype=np.float32)
-                ref_data = np.load(name)
-                np.testing.assert_allclose(actual_data, ref_data, rtol=5e-5, atol=5e-8)
-        # we expect the results to be read from cache
-        acoular.config.global_caching = 'all'
-        for b in fbeamformers():
-            b.cached = True
-            with self.subTest(b.__class__.__name__+" global_caching = all"):
-                name = join('reference_data',f'{b.__class__.__name__}.npy')
-                actual_data = np.array([b.synthetic(cf,1) for cf in cfreqs],dtype=np.float32)
-                ref_data = np.load(name)
-                np.testing.assert_allclose(actual_data, ref_data, rtol=5e-5, atol=5e-8)
-        # we expect the cached results to be overwritten
-        acoular.config.global_caching = 'overwrite'
-        for b0,b1 in zip(fbeamformers(),fbeamformers()):
-            b0.cached = True
-            b1.cached = True
-            with self.subTest(b0.__class__.__name__+" global_caching = overwrite"):
-                if hasattr(b0,'beamformer'): # BeamformerClean, BeamformerDamas, BeamformerDamasplus 
-                                            # do not pass because the .beamformer result is not take from cache
-                    continue                         # nor recalculated   
-                b0.result[:] = 0
-                self.assertFalse(np.any(b0.result))
-                name = join('reference_data',f'{b1.__class__.__name__}.npy')
-                actual_data = np.array([b1.synthetic(cf,1) for cf in cfreqs],dtype=np.float32)
-                ref_data = np.load(name)
-                np.testing.assert_allclose(actual_data, ref_data, rtol=5e-5, atol=5e-8)
-
-    def test_beamformer_caching(self):
-        # within each subcase, we need new beamformer objects because result is not updated when
-        # global_caching or cached changes
-        with self.subTest("global_caching = 'none'"):
-            acoular.config.global_caching = 'none'
-            b0 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            b1 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            self.assertNotEqual(id(b0.result),id(b1.result))
-        with self.subTest("global_caching = 'individual'"):
-            acoular.config.global_caching = 'individual'
-            b0 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            b1 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            self.assertEqual(id(b0.result),id(b1.result))
-            b1 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = False)
-            self.assertNotEqual(id(b0.result),id(b1.result))
-        with self.subTest("global_caching = 'all'"):
-            acoular.config.global_caching = 'all'
-            b0 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            b1 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = False)
-            self.assertEqual(id(b0.result),id(b1.result))
-        with self.subTest("global_caching = 'readonly'"):
-            acoular.config.global_caching = 'readonly'
-            b0 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            b1 = BeamformerBase(freq_data=f, steer=st, r_diag=True, cached = True)
-            self.assertEqual(id(b0.result),id(b1.result))
-
-class Test_PowerSpectra(unittest.TestCase):
-
-    def test_csm(self):
-        """ test that csm result has not changed over different releases"""
-        name = join('reference_data',f'{f.__class__.__name__}_csm.npy')
-        # test only two frequencies
-        actual_data = np.array(f.csm[(16,32),:,:],dtype=np.complex64)
-        if WRITE_NEW_REFERENCE_DATA:
-            np.save(name,actual_data)
-        ref_data = np.load(name)
-        np.testing.assert_allclose(actual_data, ref_data, rtol=1e-5, atol=1e-8)
-
-    def test_ev(self):
-        """ test that eve and eva result has not changed over different releases"""
-        name = join('reference_data',f'{f.__class__.__name__}_ev.npy')
-        # test only two frequencies
-        actual_data = np.array((f.eve*f.eva[:,:,np.newaxis])[(16,32),:,:],dtype=np.complex64)
-        if WRITE_NEW_REFERENCE_DATA:
-            np.save(name,actual_data)
-        ref_data = np.load(name)
-        np.testing.assert_allclose(actual_data, ref_data, rtol=1e-5, atol=1e-8)
-
-
-class TestSteerFormulation(unittest.TestCase):
-
-    def test_all_steer_formulation(self):
-        """ tests all variants of beamformerFreq subroutines
-        """
-        st = SteeringVector(grid=g, mics=m, env=env)
-        b0 = BeamformerBase(freq_data=f, steer=st, cached = False)
-        b1 = BeamformerEig(freq_data=f, steer=st, n=54, cached = False)
-        for ki,kind in enumerate(('classic', 'inverse','true level', 'true location')):
-            st.steer_type = kind
-            for b in (b0,b1):
-                for dr in (True,False):
-                    b.r_diag = dr
-                    with self.subTest(f"{b.__class__.__name__} r_diag:{dr} steer:{kind}"):
-                        name = join('reference_data',f"{b.__class__.__name__}{dr}{ki+1}.npy")
-                        actual_data = np.array([b.synthetic(cf,1) for cf in cfreqs],dtype=np.float32)
-                        if WRITE_NEW_REFERENCE_DATA:
-                            np.save(name,actual_data)
-                        ref_data = np.load(name)
-                        np.testing.assert_allclose(actual_data, ref_data, rtol=1e-5, atol=1e-8)
-
-
-if __name__ == '__main__':
-    unittest.main() #exit=False
-
-
-
-

acoular/tests/test_classes.py

@@ -1,60 +0,0 @@
-# -*- coding: utf-8 -*-
-#pylint: disable-msg=E0611, E1101, C0103, R0901, R0902, R0903, R0904, W0232
-#------------------------------------------------------------------------------
-# Copyright (c) Acoular Development Team.
-#------------------------------------------------------------------------------
-"""Implements testing of frequency beamformers.
-"""
-
-import unittest
-import acoular as ac
-from traits.api import HasTraits, Int, Float, Bool, Range, TraitEnum, Enum
-
-class Test_Instancing(unittest.TestCase):
-    """Test that ensures that digest of Acoular classes changes correctly on
-    changes of CArray and List attributes.
-    """
-
-    def test_instancing(self):
-        """ test that all Acoular classes can be instatiated """
-        # iterate over all Acoular definitions labels
-        for i in dir(ac):
-            with self.subTest(i):
-                j = getattr(ac,i) # class, function or variable
-                if isinstance(j,type): # is this a class ?
-                    j() # this is an instance of the class
-
-    def test_set_traits(self):
-        """ test that important traits can be set"""
-        # iterate over all Acoular definitions labels
-        for i in dir(ac):
-            j = getattr(ac,i) # class, function or variable
-            # HasTraits derived class ?
-            if isinstance(j,type) \
-                and issubclass(j,HasTraits) \
-                and ('digest' in j.class_traits().keys()):
-                do = j.class_traits()['digest'].depends_on
-                if do:
-                    obj = j()
-                    for k in do:
-                        with self.subTest(i+'.'+k):
-                            if k in j.class_trait_names():
-                                tr = j.class_traits()[k]
-                                # handling different Trait types
-                                # todo: use hypothesis based setattr
-                                if tr.is_trait_type(Int):
-                                    setattr(obj,k,1)
-                                elif tr.is_trait_type(Float):
-                                    setattr(obj,k,0.1)
-                                elif tr.is_trait_type(Bool):
-                                    setattr(obj,k,False)
-                                elif tr.is_trait_type(Range):
-                                    low = tr.handler._low
-                                    high = tr.handler._high
-                                    setattr(obj,k,(high+low)/2)
-                                elif tr.is_trait_type(TraitEnum) or tr.is_trait_type(Enum):
-                                    v = tr.handler.values
-                                    setattr(obj,k,v[len(v)//2])
-
-if __name__ == '__main__':
-    unittest.main()

acoular/tests/test_digest.py

@@ -1,125 +0,0 @@
-from acoular import UniformFlowEnvironment,SlotJet,OpenJet,RotatingFlow,MicGeom, SteeringVector,\
-    BeamformerOrth, RectGrid, MaskedTimeSamples, Sector,\
-        UncorrelatedNoiseSource, SourceMixer, SamplesGenerator, BeamformerTimeTraj, BeamformerTimeSqTraj,\
-            BeamformerCleantTraj, BeamformerCleantSqTraj, IntegratorSectorTime, MaskedTimeInOut, ChannelMixer,\
-                SpatialInterpolator, SpatialInterpolatorRotation, SpatialInterpolatorConstantRotation, Mixer,\
-                    WriteWAV, MergeGrid, FiltWNoiseGenerator, SphericalHarmonicSource, PointSource
-from numpy import array
-from unittest import TestCase
-
-# a dictionary containing all classes that should change their digest on
-# changes of the following trait types:
-#   * List
-#   * CArray 
-UNEQUAL_DIGEST_TEST_DICT = {
-#    "MicGeom.mpos_tot item assignment" : (MicGeom(mpos_tot=[[1.,2.,3.]]), "obj.mpos_tot[:] = 0."),
-    "MicGeom.mpos_tot new array assignment" : (MicGeom(mpos_tot=[[1.,2.,3.]]), "obj.mpos_tot = array([0.])"),
-#    "MicGeom.invalid_channels item assignment" : (MicGeom(mpos_tot=[[1.,2.,3.]],invalid_channels=[1]), "obj.invalid_channels[0] = 0"),
-    "MicGeom.invalid_channels new list assignment" : (MicGeom(mpos_tot=[[1.,2.,3.]],invalid_channels=[1]), "obj.invalid_channels = [0]"),
-    # environments.py
-#    "UniformFlowEnvironment.fdv item assignment": (UniformFlowEnvironment(), "obj.fdv[0] = 0."),
-    "UniformFlowEnvironment.fdv array assignment": (UniformFlowEnvironment(), "obj.fdv = array((0., 0., 0.))"),
-#    "SlotJet.origin item assignment": (SlotJet(), "obj.origin[0] = 1."),
-    "SlotJet.origin array assignment": (SlotJet(), "obj.origin = array((1., 0., 0.))"),
-#    "SlotJet.flow item assignment": (SlotJet(), "obj.flow[0] = 0."),
-    "SlotJet.flow array assignment": (SlotJet(), "obj.flow = array((0., 0., 0.))"),    
-#    "SlotJet.plane item assignment": (SlotJet(), "obj.plane[0] = 1."),
-    "SlotJet.plane array assignment": (SlotJet(), "obj.plane = array((1., 0., 0.))"),
-#    "OpenJet.origin item assignment": (OpenJet(), "obj.origin[0] = 1."),
-    "OpenJet.origin array assignment": (OpenJet(), "obj.origin = array((1., 0., 0.))"),   
-#    "RotatingFlow.origin item assignment": (RotatingFlow(), "obj.origin[0] = 1."),
-    "RotatingFlow.origin array assignment": (RotatingFlow(), "obj.origin = array((1., 0., 0.))"),  
-    #fbeamform.py
-#    "SteeringVector.ref item assignment": (SteeringVector(), "obj.ref[0] = 1."),  
-    "SteeringVector.ref array assignment": (SteeringVector(), "obj.ref = array((1., 1., 1.))"),  
-#    "BeamformerOrth.eva_list item assignment": (BeamformerOrth(eva_list=array((0, 1))), "obj.eva_list[0] = 2"),  
-    "BeamformerOrth.eva_list array assignment": (BeamformerOrth(eva_list=array((0, 1))), "obj.eva_list = array((2))"),  
-    #grids.py
-    "MergeGrid.grids item assignment": (MergeGrid(grids=[RectGrid()]), "obj.grids[0] = RectGrid()"),  
-    "MergeGrid.grids list assignment": (MergeGrid(), "obj.grids = [RectGrid()]"),
-    # signals.py
-#    "FiltWNoiseGenerator.ar item assignment": (FiltWNoiseGenerator(ar=[1.,2.,3.]), "obj.ar[0] = 0."),
-    "FiltWNoiseGenerator.ar array assignment": (FiltWNoiseGenerator(ar=[1.,2.,3.]), "obj.ar = array((0., 0., 0.))"),   
-#    "FiltWNoiseGenerator.ma item assignment": (FiltWNoiseGenerator(ma=[1.,2.,3.]), "obj.ma[0] = 0."),
-    "FiltWNoiseGenerator.mar array assignment": (FiltWNoiseGenerator(ma=[1.,2.,3.]), "obj.ma = array((0., 0., 0.))"), 
-    #sources.py 
-    "MaskedTimeSamples.invalid_channels item assignment": (MaskedTimeSamples(invalid_channels=[1]), "obj.invalid_channels[0] = 0"),  
-    "MaskedTimeSamples.invalid_channels list assignment": (MaskedTimeSamples(), "obj.invalid_channels = [0]"),
-#    "SphericalHarmonicSource.alpha item assignment": (SphericalHarmonicSource(alpha=array((0, 1))), "obj.alpha[0] = 1."),  
-    "SphericalHarmonicSource.alpha array assignment": (SphericalHarmonicSource(alpha=array((0, 1))), "obj.alpha = array((1., 1., 1.))"),  
-#    "UncorrelatedNoiseSource.seed item assignment": (UncorrelatedNoiseSource(seed=array((1, 2))), "obj.seed[0] = 3"),  
-    "UncorrelatedNoiseSource.seed array assignment": (UncorrelatedNoiseSource(seed=array((1, 2))), "obj.seed = array((3,4))"),  
-    "SourceMixer.sources item assignment": (SourceMixer(sources=[SamplesGenerator()]), "obj.sources[0] = PointSource()"),  
-    "SourceMixer.sources list assignment": (SourceMixer(sources=[SamplesGenerator()]), "obj.sources = [PointSource()]"),  
-    # tbeamform.py
-#    "BeamformerTimeTraj.rvec item assignment": (BeamformerTimeTraj(), "obj.rvec[0] = 1."),
-    "BeamformerTimeTraj.rvec array assignment": (BeamformerTimeTraj(), "obj.rvec = array((1., 0., 0.))"),
-#    "BeamformerTimeSqTraj.rvec item assignment": (BeamformerTimeSqTraj(), "obj.rvec[0] = 1."),
-    "BeamformerTimeSqTraj.rvec array assignment": (BeamformerTimeSqTraj(), "obj.rvec = array((1., 0., 0.))"),
-#    "BeamformerCleantTraj.rvec item assignment": (BeamformerCleantTraj(), "obj.rvec[0] = 1."),
-    "BeamformerCleantTraj.rvec array assignment": (BeamformerCleantTraj(), "obj.rvec = array((1., 0., 0.))"),
-#    "BeamformerCleantSqTraj.rvec item assignment": (BeamformerCleantSqTraj(), "obj.rvec[0] = 1."),
-    "BeamformerCleantSqTraj.rvec array assignment": (BeamformerCleantSqTraj(), "obj.rvec = array((1., 0., 0.))"),
-    "IntegratorSectorTime.sectors item assignment": (IntegratorSectorTime(sectors=[Sector()]), "obj.sectors[0] = Sector()"),  
-    "IntegratorSectorTime.sectors list assignment": (IntegratorSectorTime(sectors=[Sector()]), "obj.sectors = [Sector()]"),  
-    # tprocess.py
-    "MaskedTimeInOut.invalid_channels item assignment": (MaskedTimeInOut(invalid_channels=[1]), "obj.invalid_channels[0] = 0"),  
-    "MaskedTimeInOut.invalid_channels list assignment": (MaskedTimeInOut(), "obj.invalid_channels = [0]"),
-#    "ChannelMixer.weights item assignment" : (ChannelMixer(weights=[[1.,2.,3.]]), "obj.weights[0] = 0."),
-    "ChannelMixer.weights new array assignment" : (ChannelMixer(weights=[[1.,2.,3.]]), "obj.weights = array([0.])"),
-#    "SpatialInterpolator.Q item assignment": (SpatialInterpolator(), "obj.Q[0] = 0."),
-    "SpatialInterpolator.Q array assignment": (SpatialInterpolator(), "obj.Q = array([[0.,0.,0.],[0.,0.,0.],[0.,0.,0.]])"),
-#    "SpatialInterpolatorRotation.Q item assignment": (SpatialInterpolatorRotation(), "obj.Q[0] = 0."),
-    "SpatialInterpolatorRotation.Q array assignment": (SpatialInterpolatorRotation(), "obj.Q = array([[0.,0.,0.],[0.,0.,0.],[0.,0.,0.]])"),
-#    "SpatialInterpolatorConstantRotation.Q item assignment": (SpatialInterpolatorConstantRotation(), "obj.Q[0] = 0."),
-    "SpatialInterpolatorConstantRotation.Q array assignment": (SpatialInterpolatorConstantRotation(), "obj.Q = array([[0.,0.,0.],[0.,0.,0.],[0.,0.,0.]])"),
-    "Mixer.sources item assignment": (Mixer(source=SamplesGenerator(), sources=[SamplesGenerator()]), "obj.sources[0] = SamplesGenerator(numchannels=1)"),  
-    "Mixer.sources list assignment": (Mixer(sources=[SamplesGenerator()]), "obj.sources = [SamplesGenerator(), SamplesGenerator()]"),  
-    "WriteWAV.channels item assignment": (WriteWAV(channels=[1]), "obj.channels[0] = 0"),  
-    "WriteWAV.channels list assignment": (WriteWAV(), "obj.channels = [0]"),
-    }
-
-
-class Test_DigestChange(TestCase):
-    """Test that ensures that digest of Acoular classes changes correctly on 
-    changes of CArray and List attributes.
-    """
-
-    def get_digests(self,obj,statement):
-        """A function that collects the digest of the obj before (d1) and
-        after (d2) executing the statement that should yield to a change
-        of the object digest.
-
-        Parameters
-        ----------
-        obj : instance
-            class instance that has an attribute `digest`
-        statement : "str"
-            a string that can be executed 
-
-        Returns
-        -------
-        str, str
-            digest before and after statement execution
-        """
-        d1 = obj.digest
-        exec(statement)
-        d2 = obj.digest
-        return d1,d2
-
-    def test_digest_changes_on_assignment(self):
-        """ test that object digest change on statement execution """
-        for test_label,(obj,statement) in UNEQUAL_DIGEST_TEST_DICT.items():
-            # use subtest to get individual test result.
-            # tests will not stop if one subtest fails!
-            with self.subTest(test_label):
-                digest_before_statement,digest_after_statement = self.get_digests(obj, statement)
-                self.assertNotEqual(digest_before_statement,digest_after_statement)        
-
-
-
-if __name__ == '__main__':
-    import unittest 
-    unittest.main()
-
-

acoular/tests/test_environments.py

@@ -1,73 +0,0 @@
-# -*- coding: utf-8 -*-
-#pylint: disable-msg=E0611, E1101, C0103, R0901, R0902, R0903, R0904, W0232
-#------------------------------------------------------------------------------
-# Copyright (c) Acoular Development Team.
-#------------------------------------------------------------------------------
-"""Implements testing of environments and helpers.
-"""
-
-import unittest
-
-from os.path import join
-
-import numpy as np
-#acoular imports
-import acoular
-acoular.config.global_caching = 'none' # to make sure that nothing is cached
-
-from acoular import MicGeom, RectGrid3D, \
-    SlotJet, OpenJet, RotatingFlow, \
-    Environment, UniformFlowEnvironment, GeneralFlowEnvironment
-
-# if this flag is set to True
-WRITE_NEW_REFERENCE_DATA = False
-# results are generated for comparison during testing.
-# Should always be False. Only set to True if it is necessary to
-# recalculate the data due to intended changes of the Beamformers.
-
-m = MicGeom()
-m.mpos_tot = ((0.5,0.5,0),(0,0,0),(-0.5,-0.5,0))
-mc = m.mpos
-g = RectGrid3D(x_min=-0.2, x_max=0.2, y_min=-0.2, y_max=0.2, z_min=0.5, z_max=0.9,
-             increment=0.2 )
-gc = g.gpos
-
-flows = [
-    SlotJet(v0=70.0, origin=(-0.7,0,0.7)),
-    OpenJet(v0=70.0, origin=(-0.7,0,0.7)),
-    RotatingFlow(v0=70.0, rpm=1000.0)
-]
-
-envs = [
-    Environment(),
-    UniformFlowEnvironment(ma=0.3),
-    GeneralFlowEnvironment(ff=OpenJet(v0=70.0, origin=(-0.7,0,0.7)))
-]
-
-
-class acoular_env_test(unittest.TestCase):
-
-    def test_flow_results(self):
-        for fl in flows:
-            with self.subTest(fl.__class__.__name__):
-                name = join('reference_data',f'{fl.__class__.__name__}.npy')
-                # stack all results
-                actual_data = np.array([np.vstack(fl.v(x)) for x in gc.T])
-                if WRITE_NEW_REFERENCE_DATA:
-                    np.save(name,actual_data)
-                ref_data = np.load(name)
-                np.testing.assert_allclose(actual_data, ref_data, rtol=1e-5, atol=1e-8)
-
-    def test_env_results(self):
-        for env in envs:
-            with self.subTest(env.__class__.__name__):
-                name = join('reference_data',f'{env.__class__.__name__}.npy')
-                # stack all results
-                actual_data = np.vstack((env._r(gc,mc).T,env._r(gc).T))
-                if WRITE_NEW_REFERENCE_DATA:
-                    np.save(name,actual_data)
-                ref_data = np.load(name)
-                np.testing.assert_allclose(actual_data, ref_data, rtol=1e-5, atol=1e-8)
-
-if __name__ == '__main__':
-    unittest.main()