def test_ReadEnsightGeo(self): parts = [1, 2] partRec = [] THIS_FOLDER = os.path.dirname(os.path.abspath(__file__)) my_file = os.path.join(THIS_FOLDER, 'testEnsightReader.geo') mesh = EnsightReader2.readEnsightGeo(my_file, parts, partRec) # mesh = EnsightReader2.readEnsightGeo('MMPTestCase_v1.geo', parts, partRec) timeUnits = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0]) temperatureUnits = PQ.PhysicalUnit('K', 1., [0, 0, 0, 0, 1, 0, 0, 0, 0]) # f = EnsightReader2.readEnsightField('fld_TEMPERATURE.escl', parts, partRec, 1, FieldID.FID_Temperature, mesh, temperatureUnits, 0) self.assertTrue(1 == 1)
def test_cmp(self): timeUnits = PQ.PhysicalUnit('s', 1., [0,0,1,0,0,0,0,0,0]) timeUnits2 = PQ.PhysicalUnit('s', 2., [0, 0, 1, 0, 0, 0, 0, 0, 0]) timeUnits3 = PQ.PhysicalUnit('s', -1., [0, 0, 1, 0, 0, 0, 0, 0, 0]) timeUnits4 = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0]) temperatureUnits = PQ.PhysicalUnit('K', 1., [0,0,0,0,1,0,0,0,0]) try: a = timeUnits.__cmp__(temperatureUnits) except TypeError: pass except Exception as e: self.fail('Unexpected exception raised:', e) else: self.fail('Exception not raised') self.assertTrue(timeUnits.__cmp__(timeUnits2) < 0) self.assertTrue(timeUnits.__cmp__(timeUnits3) > 0) self.assertTrue(timeUnits.__cmp__(timeUnits4) == 0) self.assertTrue(timeUnits.__cmp__(5) == -1)
#!/usr/bin/env python from __future__ import print_function from builtins import str import sys sys.path.append('../../..') from mupif import * import logging log = logging.getLogger() import mupif.Physics.PhysicalQuantities as PQ timeUnits = PQ.PhysicalUnit('s', 1., [0,0,1,0,0,0,0,0,0]) class application1(Application.Application): """ Simple application that generates a property with a value equal to actual time """ def __init__(self, file): super(application1,self).__init__(self,file) def getProperty(self, propID, time, objectID=0): if (propID == PropertyID.PID_Velocity): return Property.ConstantProperty(self.value, PropertyID.PID_Velocity, ValueType.Scalar, 'm/s', time, 0) else: raise APIError.APIError ('Unknown property ID') def solveStep(self, tstep, stageID=0, runInBackground=False): time = tstep.getTime().inUnitsOf(timeUnits).getValue() self.value=1.0*time def getCriticalTimeStep(self):
import sys sys.path.append('../../..') from mupif import * import mupif.Physics.PhysicalQuantities as PQ timeUnits = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0]) temperatureUnits = PQ.PhysicalUnit('K', 1., [0, 0, 0, 0, 1, 0, 0, 0, 0]) class Celsian(Application.Application): def __init__(self, file): super(Celsian, self).__init__(file) self.mesh = None def getField(self, fieldID, time): parts = [1, 2] partRec = [] if (self.mesh == None): self.mesh = EnsightReader2.readEnsightGeo( 'paraview/MMPTestCase_v1.geo', parts, partRec) f = EnsightReader2.readEnsightField('paraview/fld_TEMPERATURE.escl', parts, partRec, 1, FieldID.FID_Temperature, self.mesh, temperatureUnits, time) return f def solveStep(self, tstep, stageID=0, runInBackground=False): time = tstep.getTime() self.value = 1.0 * time
#!/usr/bin/env python from builtins import range # # This example requires pyvtk module, install it using # pip install pyvtk # Note: pyvtk is only available in Python 2.x (Sept, 2014) # import sys import logging sys.path.append('../../..') from mupif import * import mupif.Physics.PhysicalQuantities as PQ log = logging.getLogger() temperatureUnit = PQ.PhysicalUnit('K', 1., [0, 0, 0, 0, 1, 0, 0, 0, 0]) def main(): # generate field and corresponding mesh mesh = Mesh.UnstructuredMesh() vertices = [] values1 = [] values2 = [] num = 0 for i in range(40): for j in range(15): vertices.append( Vertex.Vertex((i * 15) + j, (i * 15) + j + 1, (float(i), float(j), 0.0))) values1.append((num, ))
import sys sys.path.append('../../..') # for mupif import import os import platform from subprocess import * import socket import time as timeTime import pyvtk from mupif import * import logging import mupif.Physics.PhysicalQuantities as PQ dummyUnits = PQ.PhysicalUnit('dummy', 1., [0, 0, 0, 0, 0, 0, 0, 0, 0]) timeUnits = PQ.PhysicalUnit('s', 1., [0, 0, 1, 0, 0, 0, 0, 0, 0]) import micressConfig as mConf #Micress import #import Property import MICPropertyID import MICFieldID VtkReader2.pyvtk_monkeypatch() class micress(Application.Application): """ MICRESS application interface class