def test_messagefederate_test_message_federate_send(mFed):
epid1 = h.helicsFederateRegisterEndpoint(mFed, "ep1", "")
epid2 = h.helicsFederateRegisterGlobalEndpoint(mFed, "ep2", "random")
h.helicsFederateSetTimeProperty(mFed, h.HELICS_PROPERTY_TIME_DELTA, 1.0)
h.helicsFederateEnterExecutingMode(mFed)
data = "random-data"
h.helicsEndpointSendEventRaw(epid1, "ep2", data, 1.0)
granted_time = h.helicsFederateRequestTime(mFed, 2.0)
assert granted_time == 1.0
res = h.helicsFederateHasMessage(mFed)
assert res is True
res = h.helicsEndpointHasMessage(epid1)
assert res is False
res = h.helicsEndpointHasMessage(epid2)
assert res is True
message = h.helicsEndpointGetMessage(epid2)
assert h.helicsMessageGetMessageID(message) == 55
assert h.helicsMessageIsValid(message) is True
assert h.helicsMessageGetString(message) == "random-data"
assert h.helicsMessageGetRawDataSize(message) == 11
assert h.helicsMessageGetOriginalDestination(message) == ""
assert h.helicsMessageGetOriginalSource(message) == "TestA Federate/ep1"
assert h.helicsMessageGetSource(message) == "TestA Federate/ep1"
assert h.helicsMessageGetTime(message) == 1.0
def test_filter_test_types_clone_test_broker_connections():
broker = createBroker(3)
sFed, fedinfo1 = createMessageFederate(1, "source", 1.0)
dFed, fedinfo2 = createMessageFederate(1, "dest", 1.0)
dcFed, fedinfo3 = createMessageFederate(1, "dest_clone", 1.0)
p1 = h.helicsFederateRegisterGlobalEndpoint(sFed, "src", "")
p2 = h.helicsFederateRegisterGlobalEndpoint(dFed, "dest", "")
p3 = h.helicsFederateRegisterGlobalEndpoint(dcFed, "cm", "")
f1 = h.helicsFederateRegisterGlobalCloningFilter(dcFed, "filt1")
h.helicsFilterAddDeliveryEndpoint(f1, "cm")
h.helicsBrokerAddSourceFilterToEndpoint(broker, "filt1", "src")
h.helicsBrokerAddSourceFilterToEndpoint(broker, "", "src")
h.helicsFederateEnterExecutingModeAsync(sFed)
h.helicsFederateEnterExecutingModeAsync(dcFed)
h.helicsFederateEnterExecutingMode(dFed)
h.helicsFederateEnterExecutingModeComplete(sFed)
h.helicsFederateEnterExecutingModeComplete(dcFed)
state = h.helicsFederateGetState(sFed)
assert state == h.HELICS_STATE_EXECUTION
data = "".join(["a" for _ in range(0, 500)]).encode()
h.helicsEndpointSendMessageRaw(p1, "dest", data)
h.helicsFederateRequestTimeAsync(sFed, 1.0)
h.helicsFederateRequestTimeAsync(dcFed, 1.0)
h.helicsFederateRequestTime(dFed, 1.0)
h.helicsFederateRequestTimeComplete(sFed)
h.helicsFederateRequestTimeComplete(dcFed)
assert h.helicsFederateHasMessage(dFed) is True
m2 = h.helicsEndpointGetMessage(p2)
assert h.helicsMessageGetSource(m2) == "src"
assert h.helicsMessageGetOriginalSource(m2) == "src"
assert h.helicsMessageGetDestination(m2) == "dest"
assert h.helicsMessageGetRawDataSize(m2) == len(data)
assert h.helicsFederateHasMessage(dcFed) is True
m2 = h.helicsEndpointGetMessage(p3)
assert h.helicsMessageGetSource(m2) == "src"
assert h.helicsMessageGetOriginalSource(m2) == "src"
assert h.helicsMessageGetDestination(m2) == "cm"
assert h.helicsMessageGetOriginalDestination(m2) == "dest"
assert h.helicsMessageGetRawDataSize(m2) == len(data)
h.helicsFederateFinalizeAsync(sFed)
h.helicsFederateFinalizeAsync(dFed)
h.helicsFederateFinalize(dcFed)
h.helicsFederateFinalizeComplete(sFed)
h.helicsFederateFinalizeComplete(dFed)
state = h.helicsFederateGetState(sFed)
assert state == h.HELICS_STATE_FINALIZE
destroyFederate(sFed, fedinfo1)
destroyFederate(dFed, fedinfo2)
destroyFederate(dcFed, fedinfo3)
destroyBroker(broker)
def run_cosim(fed, endid):
# The event queue ("eq") is the master list of events that the filter
# federates works on. In this simple filter federate, each event
# will be a dictionary with a few parameters:
# 'dest' - destination of message
# 'time' - time when the message should be sent on to its
# intended destination
# 'payload' - content of the message being sent
#
# When eq is empty, there are no messages being
# filtered by the federate. When there are events in the queue it
# indicates the filter federate has messages it is holding onto
# that it needs to forward on (at the indicated time).
#
eq = []
h.helicsFederateEnterExecutingMode(fed)
logger.info('Entered HELICS execution mode')
hours = 24 * 7 # one week
total_interval = int(60 * 60 * hours)
# Blocking call for a time request at max simulation time
fake_max_time = int(h.HELICS_TIME_MAXTIME / 1000)
starttime = fake_max_time
logger.debug(f'Requesting initial time {starttime}')
grantedtime = h.helicsFederateRequestTime(fed, starttime)
logger.debug(f'Granted time {grantedtime}')
while grantedtime < total_interval:
# In HELICS, when multiple messages arrive at an endpoint they
# queue up and are popped off one-by-one with the
# "helicsEndpointHasMessage" API call. When that API doesn't
# return a message, you've processed them all.
while h.helicsEndpointHasMessage(endid):
msg = h.helicsEndpointGetMessage(endid)
msg_str = h.helicsMessageGetString(msg)
source = h.helicsMessageGetOriginalSource(msg)
dest = h.helicsMessageGetOriginalDestination(msg)
time = h.helicsMessageGetTime(msg)
logger.debug(f'Received message from endpoint {source}'
f' to endpoint {dest}'
f' at time {grantedtime}'
f' with message {msg_str}')
event = {
"payload": msg_str,
"source": source,
"dest": dest,
"time": time
}
eq = filter_message(eq, event, 'delay')
if source == 'EVController_federate/endpoint':
eq = filter_message(eq, event, 'hack')
# Sort event queue to get it back in order
eq = sorted(eq, key=itemgetter('time'))
# Acting on any events that need to be dequeued
# Running interference filter. This filter has the ability to
# remove events from eq. We may not have any messages to send
# after interference runs
event = eq[0]
eq = filter_message(eq, event, 'interfere')
# After filtering, send all messages whose time has come (or past;
# in which case something has gone wrong)
while eq and eq[0]['time'] <= grantedtime:
h.helicsEndpointSendMessageRaw(endid, eq[0]['dest'],
eq[0]['payload'].encode())
logger.debug(f'Sent message from endpoint {endid}'
f' to endpoint {eq[0]["dest"]}'
f' at time {grantedtime}'
f' with message {eq[0]["payload"]}')
del eq[0]
if eq:
# Event queue not empty, need to schedule filter federate to
# run again when its time to deliver the next message in the
# queue
requested_time = eq[0]['time']
else:
# If no events in queue, schedule run for end of simulation.
# Filter federate will be granted an earlier time if a
# message is rerouted to the filter federate.
requested_time = fake_max_time
logger.debug(f'Requesting time {requested_time}')
grantedtime = h.helicsFederateRequestTime(fed, requested_time)
logger.debug(f'Granted time {grantedtime}')
def run_cosim(fed, endid, end_name, args):
# The event queue ("eq") is the master list of events that the filter
# federates works on. In this simple filter federate, each event
# will be a dictionary with a few parameters:
# 'dest' - destination of message
# 'time' - time when the message should be sent on to its
# intended destination
# 'payload' - content of the message being sent
#
# When eq is empty, there are no messages being
# filtered by the federate. When there are events in the queue it
# indicates the filter federate has messages it is holding onto
# that it needs to forward on (at the indicated time).
#
eq = []
# sub = h.helicsFederateRegisterSubscription(fed, "Charger/EV1_voltage", "")
logger.info('Attempting to enter execution mode')
h.helicsFederateEnterExecutingMode(fed)
logger.info('Entered HELICS execution mode')
hours = 24 * 7 # one week
total_interval = int(60 * 60 * hours)
# Blocking call for a time request at max simulation time
# fake_max_time = int(h.HELICS_TIME_MAXTIME / 1000)
fake_max_time = 100000000
starttime = fake_max_time
#starttime = 0
logger.debug(f'Requesting initial time {starttime}')
grantedtime = h.helicsFederateRequestTime(fed, starttime)
logger.debug(f'Granted time {grantedtime}')
while grantedtime < total_interval:
# value = h.helicsInputGetString(sub)
# logger.debug(f'Got message {value} from random sub at time {grantedtime}.')
# In HELICS, when multiple messages arrive at an endpoint they
# queue up and are popped off one-by-one with the
# "helicsEndpointHasMessage" API call. When that API doesn't
# return a message, you've processed them all.
while h.helicsEndpointHasMessage(endid):
msg = h.helicsEndpointGetMessage(endid)
msg_str = h.helicsMessageGetString(msg)
source = h.helicsMessageGetOriginalSource(msg)
dest = h.helicsMessageGetOriginalDestination(msg)
time = h.helicsMessageGetTime(msg)
logger.debug(f'Received message from endpoint {source}'
f' to endpoint {dest}'
f' for delivery at time {time}'
f' with payload \"{msg_str}\"')
msg_dict = {
'msg_obj': msg,
'payload': msg_str,
'source': source,
'dest': dest,
'time': time
}
# Adding messagge to end of eq as a reserved place for all
# filters to act on.
eq.append(msg_dict)
eq = filter_message(eq, 'drop_delay', args)
if source == 'Controller/ep':
eq = filter_message(eq, 'hack', args)
# Sort event queue to get it back in order
eq = sorted(eq, key=itemgetter('time'))
# Acting on any events that need to be dequeued
# Running interference filter. This filter has the ability to
# remove events from eq. We may not have any messages to send
# after interference runs. eq must be freshly sorted for this
# filter to work.
if len(eq) > 0:
eq = filter_message(eq, 'interfere', args)
# After filtering, send all messages whose time has come (or past;
# in which case something has gone wrong)
while eq and eq[0]['time'] <= grantedtime:
# Change destination to original destination before sending
# If you don't do this is sends the message back to the rerouted
# destination which, in this case, is the filter endpoint.
h.helicsMessageSetDestination(eq[0]['msg_obj'], eq[0]["dest"])
h.helicsEndpointSendMessage(endid, eq[0]['msg_obj'])
# h.helicsEndpointSendMessageRaw(endid, eq[0]['dest'],
# eq[0]['payload'].encode())
# logger.debug(eq[0]['msg_obj'])
logger.debug(f'Sent message from endpoint {end_name}'
f' appearing to come from {eq[0]["source"]}'
f' to endpoint {eq[0]["dest"]}'
f' at time {grantedtime}'
f' with payload \"{eq[0]["payload"]}\"')
del eq[0]
if eq:
# Event queue not empty, need to schedule filter federate to
# run again when its time to deliver the next message in the
# queue
requested_time = eq[0]['time']
else:
# Reachable if interference has removed all the messages
# from the event queue.
# No events in queue, schedule run for end of simulation.
# Filter federate will be granted an earlier time if a
# message is rerouted to the filter federate.
requested_time = fake_max_time
#requested_time = grantedtime + 5
else:
requested_time = fake_max_time
#requested_time = grantedtime + 5
logger.debug(f'Requesting time {requested_time}\n')
grantedtime = h.helicsFederateRequestTime(fed, requested_time)
logger.debug(f'Granted time {grantedtime}')
def test_misc_api():
fedInfo1 = h.helicsCreateFederateInfo()
h.helicsFederateInfoSetCoreInitString(fedInfo1, "-f 1")
h.helicsFederateInfoSetCoreName(fedInfo1, "core3")
h.helicsFederateInfoSetCoreType(fedInfo1, 3)
h.helicsFederateInfoSetCoreTypeFromString(fedInfo1, "zmq")
h.helicsFederateInfoSetFlagOption(fedInfo1, 1, True)
h.helicsFederateInfoSetTimeProperty(fedInfo1,
h.HELICS_PROPERTY_TIME_INPUT_DELAY,
1.0)
h.helicsFederateInfoSetIntegerProperty(fedInfo1,
h.HELICS_PROPERTY_INT_LOG_LEVEL, 1)
h.helicsFederateInfoSetIntegerProperty(
fedInfo1, h.HELICS_PROPERTY_INT_MAX_ITERATIONS, 100)
h.helicsFederateInfoSetTimeProperty(fedInfo1,
h.HELICS_PROPERTY_TIME_OUTPUT_DELAY,
1.0)
h.helicsFederateInfoSetTimeProperty(fedInfo1,
h.HELICS_PROPERTY_TIME_PERIOD, 1.0)
h.helicsFederateInfoSetTimeProperty(fedInfo1, h.HELICS_PROPERTY_TIME_DELTA,
1.0)
h.helicsFederateInfoSetTimeProperty(fedInfo1,
h.HELICS_PROPERTY_TIME_OFFSET, 0.1)
h.helicsFederateInfoFree(fedInfo1)
broker3 = h.helicsCreateBroker("zmq", "broker3",
"--federates 1 --loglevel 1")
fedInfo2 = h.helicsCreateFederateInfo()
coreInitString = "--federates 1"
h.helicsFederateInfoSetCoreInitString(fedInfo2, coreInitString)
h.helicsFederateInfoSetCoreTypeFromString(fedInfo2, "zmq")
h.helicsFederateInfoSetIntegerProperty(fedInfo2,
h.HELICS_PROPERTY_INT_LOG_LEVEL, 1)
h.helicsFederateInfoSetTimeProperty(fedInfo2, h.HELICS_PROPERTY_TIME_DELTA,
1.0)
fed1 = h.helicsCreateCombinationFederate("fed1", fedInfo2)
fed2 = h.helicsFederateClone(fed1)
_ = h.helicsGetFederateByName("fed1")
h.helicsFederateSetFlagOption(fed2, 1, False)
h.helicsFederateSetTimeProperty(fed2, h.HELICS_PROPERTY_TIME_INPUT_DELAY,
1.0)
h.helicsFederateSetIntegerProperty(fed1, h.HELICS_PROPERTY_INT_LOG_LEVEL,
1)
h.helicsFederateSetIntegerProperty(fed2,
h.HELICS_PROPERTY_INT_MAX_ITERATIONS,
100)
h.helicsFederateSetTimeProperty(fed2, h.HELICS_PROPERTY_TIME_OUTPUT_DELAY,
1.0)
h.helicsFederateSetTimeProperty(fed2, h.HELICS_PROPERTY_TIME_PERIOD, 0.0)
h.helicsFederateSetTimeProperty(fed2, h.HELICS_PROPERTY_TIME_DELTA, 1.0)
_ = h.helicsFederateRegisterCloningFilter(fed1, "fed1/Ep1")
fed1DestinationFilter = h.helicsFederateRegisterFilter(
fed1, h.HELICS_FILTER_TYPE_DELAY, "fed1DestinationFilter")
h.helicsFilterAddDestinationTarget(fed1DestinationFilter, "Ep2")
ep1 = h.helicsFederateRegisterEndpoint(fed1, "Ep1", "string")
ep2 = h.helicsFederateRegisterGlobalEndpoint(fed1, "Ep2", "string")
pub1 = h.helicsFederateRegisterGlobalPublication(fed1, "pub1",
h.HELICS_DATA_TYPE_DOUBLE,
"")
pub2 = h.helicsFederateRegisterGlobalTypePublication(
fed1, "pub2", "complex", "")
sub1 = h.helicsFederateRegisterSubscription(fed1, "pub1")
sub2 = h.helicsFederateRegisterSubscription(fed1, "pub2")
h.helicsInputAddTarget(sub2, "Ep2")
pub3 = h.helicsFederateRegisterPublication(fed1, "pub3",
h.HELICS_DATA_TYPE_STRING, "")
pub1KeyString = h.helicsPublicationGetKey(pub1)
pub1TypeString = h.helicsPublicationGetType(pub1)
pub1UnitsString = h.helicsPublicationGetUnits(pub1)
sub1KeyString = h.helicsSubscriptionGetKey(sub1)
sub1UnitsString = h.helicsInputGetUnits(sub1)
assert "pub1" == pub1KeyString
assert "double" == pub1TypeString
assert "" == pub1UnitsString
assert "pub1" == sub1KeyString
assert "" == sub1UnitsString
fed1SourceFilter = h.helicsFederateRegisterFilter(
fed1, h.HELICS_FILTER_TYPE_DELAY, "fed1SourceFilter")
h.helicsFilterAddSourceTarget(fed1SourceFilter, "Ep2")
h.helicsFilterAddDestinationTarget(fed1SourceFilter, "fed1/Ep1")
h.helicsFilterRemoveTarget(fed1SourceFilter, "fed1/Ep1")
h.helicsFilterAddSourceTarget(fed1SourceFilter, "Ep2")
h.helicsFilterRemoveTarget(fed1SourceFilter, "Ep2")
fed1SourceFilterNameString = h.helicsFilterGetName(fed1SourceFilter)
assert fed1SourceFilterNameString == "fed1/fed1SourceFilter"
sub3 = h.helicsFederateRegisterSubscription(fed1, "fed1/pub3", "")
pub4 = h.helicsFederateRegisterTypePublication(fed1, "pub4", "int", "")
sub4 = h.helicsFederateRegisterSubscription(fed1, "fed1/pub4", "")
pub5 = h.helicsFederateRegisterGlobalTypePublication(
fed1, "pub5", "boolean", "")
sub5 = h.helicsFederateRegisterSubscription(fed1, "pub5", "")
pub6 = h.helicsFederateRegisterGlobalPublication(fed1, "pub6",
h.HELICS_DATA_TYPE_VECTOR,
"")
sub6 = h.helicsFederateRegisterSubscription(fed1, "pub6", "")
pub7 = h.helicsFederateRegisterGlobalPublication(
fed1, "pub7", h.HELICS_DATA_TYPE_NAMED_POINT, "")
sub7 = h.helicsFederateRegisterSubscription(fed1, "pub7", "")
assert """helics.HelicsPublication(name = "pub1", type = "double", units = "", info = "")""" in repr(
pub1)
assert """helics.HelicsPublication(name = "pub2", type = "complex", units = "", info = "")""" in repr(
pub2)
assert """helics.HelicsPublication(name = "fed1/pub3", type = "string", units = "", info = "")""" in repr(
pub3)
assert """helics.HelicsPublication(name = "fed1/pub4", type = "int", units = "", info = "")""" in repr(
pub4)
assert """helics.HelicsPublication(name = "pub5", type = "boolean", units = "", info = "")""" in repr(
pub5)
assert """helics.HelicsPublication(name = "pub6", type = "double_vector", units = "", info = "")""" in repr(
pub6)
assert """helics.HelicsPublication(name = "pub7", type = "named_point", units = "", info = "")""" in repr(
pub7)
assert (
"""helics.HelicsInput(name = "_input_18", units = "", injection_units = "", publication_type = "", type = "", target = "pub7", info = "")"""
in repr(sub7))
h.helicsInputSetDefaultBoolean(sub5, False)
h.helicsInputSetDefaultComplex(sub2, -9.9 + 2.5j)
h.helicsInputSetDefaultDouble(sub1, 3.4)
h.helicsInputSetDefaultInteger(sub4, 6)
h.helicsInputSetDefaultNamedPoint(sub7, "hollow", 20.0)
h.helicsInputSetDefaultString(sub3, "default")
sub6Default = [3.4, 90.9, 4.5]
h.helicsInputSetDefaultVector(sub6, sub6Default)
h.helicsEndpointSubscribe(ep2, "fed1/pub3")
h.helicsFederateEnterInitializingModeAsync(fed1)
rs = h.helicsFederateIsAsyncOperationCompleted(fed1)
if rs == 0:
time.sleep(0.500)
rs = h.helicsFederateIsAsyncOperationCompleted(fed1)
if rs == 0:
time.sleep(0.500)
rs = h.helicsFederateIsAsyncOperationCompleted(fed1)
if rs == 0:
assert True is False
h.helicsFederateEnterInitializingModeComplete(fed1)
h.helicsFederateEnterExecutingModeAsync(fed1)
h.helicsFederateEnterExecutingModeComplete(fed1)
assert (
"""helics.HelicsInput(name = "_input_18", units = "", injection_units = "", publication_type = "named_point", type = "", target = "pub7", info = "")"""
in repr(sub7))
mesg1 = h.helicsFederateCreateMessage(fed1)
h.helicsMessageSetString(mesg1, "Hello")
h.helicsMessageSetSource(mesg1, "fed1/Ep1")
h.helicsMessageSetOriginalSource(mesg1, "fed1/Ep1")
h.helicsMessageSetDestination(mesg1, "Ep2")
h.helicsMessageSetOriginalDestination(mesg1, "Ep2")
h.helicsEndpointSendMessage(ep1, mesg1)
mesg1 = h.helicsFederateCreateMessage(fed1)
h.helicsMessageSetString(mesg1, "There")
h.helicsMessageSetSource(mesg1, "fed1/Ep1")
h.helicsMessageSetOriginalSource(mesg1, "fed1/Ep1")
h.helicsMessageSetDestination(mesg1, "Ep2")
h.helicsMessageSetOriginalDestination(mesg1, "Ep2")
h.helicsEndpointSendMessage(ep1, mesg1)
h.helicsEndpointSetDefaultDestination(ep2, "fed1/Ep1")
ep1NameString = h.helicsEndpointGetName(ep1)
ep1TypeString = h.helicsEndpointGetType(ep1)
assert ep1NameString == "fed1/Ep1"
assert ep1TypeString == "string"
_ = h.helicsFederateGetCoreObject(fed1)
fed1Time = h.helicsFederateGetCurrentTime(fed1)
assert fed1Time == 0.0
fed1EndpointCount = h.helicsFederateGetEndpointCount(fed1)
assert fed1EndpointCount == 2
fed1NameString = h.helicsFederateGetName(fed1)
assert fed1NameString == "fed1"
fed1State = h.helicsFederateGetState(fed1)
assert fed1State == 2
fed1PubCount = h.helicsFederateGetPublicationCount(fed1)
assert fed1PubCount == 7
fed1SubCount = h.helicsFederateGetInputCount(fed1)
assert fed1SubCount == 7
h.helicsPublicationPublishBoolean(pub5, True)
h.helicsPublicationPublishComplex(pub2, 5.6 + -0.67j)
h.helicsPublicationPublishDouble(pub1, 457.234)
h.helicsPublicationPublishInteger(pub4, 1)
h.helicsPublicationPublishNamedPoint(pub7, "Blah Blah", 20.0)
h.helicsPublicationPublishString(pub3, "Mayhem")
pub6Vector = [4.5, 56.5]
h.helicsPublicationPublishVector(pub6, pub6Vector)
time.sleep(0.500)
h.helicsFederateRequestTimeAsync(fed1, 1.0)
returnTime = h.helicsFederateRequestTimeComplete(fed1)
assert returnTime == 1.0
ep2MsgCount = h.helicsEndpointPendingMessages(ep2)
assert ep2MsgCount == 2
ep2HasMsg = h.helicsEndpointHasMessage(ep2)
assert ep2HasMsg == 1
msg2 = h.helicsEndpointGetMessage(ep2)
assert h.helicsMessageGetTime(msg2) == 1.0
assert h.helicsMessageGetString(msg2) == "Hello"
assert h.helicsMessageGetOriginalSource(msg2) == "fed1/Ep1"
assert h.helicsMessageGetSource(msg2) == "fed1/Ep1"
assert h.helicsMessageGetDestination(msg2) == "Ep2"
assert h.helicsMessageGetOriginalDestination(msg2) == "Ep2"
fed1MsgCount = h.helicsFederatePendingMessages(fed1)
assert fed1MsgCount == 1
assert h.helicsFederateHasMessage(fed1) == 1
msg3 = h.helicsFederateGetMessage(fed1)
assert h.helicsMessageGetTime(msg3) == 1.0
assert h.helicsMessageGetString(msg3) == "There"
assert h.helicsMessageGetOriginalSource(msg3) == "fed1/Ep1"
assert h.helicsMessageGetSource(msg3) == "fed1/Ep1"
assert h.helicsMessageGetDestination(msg3) == "Ep2"
assert h.helicsMessageGetOriginalDestination(msg3) == "Ep2"
sub1Updated = h.helicsInputIsUpdated(sub1)
# TODO: figure out why this test is broken
assert sub1Updated is False
# TODO: figure out why this test is broken
assert h.helicsInputLastUpdateTime(sub2) == 0.0
# assert h.helicsInputGetComplex(sub2) == 5.6 - 0.67j
# assert h.helicsInputGetDouble(sub1) == 457.234
# assert h.helicsInputGetInteger(sub4) == 1
sub7PointString, sub7DoubleValue = h.helicsInputGetNamedPoint(sub7)
# assert sub7PointString == "Blah Blah"
assert sub7DoubleValue == 20.0
# assert h.helicsInputGetBoolean(sub5) == True
# assert h.helicsInputGetString(sub3) == "Mayhem"
sub3ValueSize = h.helicsInputGetRawValueSize(sub3)
# assert sub3ValueSize == 6
# assert h.helicsInputGetVector(sub6) == [4.5, 56.5]
h.helicsFederateFinalize(fed1)
h.helicsFederateFinalize(fed2)
h.helicsFederateFree(fed1)
h.helicsFederateFinalize(fed2)
h.helicsFederateFree(fed2)
h.helicsFederateInfoFree(fedInfo2)
h.helicsBrokerDisconnect(broker3)
h.helicsBrokerFree(broker3)
h.helicsCleanupLibrary()
h.helicsCloseLibrary()
def run_cosim(fed, endid, end_name, args):
# The event queue ("eq") is the master list of events that the filter
# federates works on. In this simple filter federate, each event
# will be a dictionary with a few parameters:
# 'dest' - destination of message
# 'time' - time when the message should be sent on to its
# intended destination
# 'payload' - content of the message being sent
#
# When eq is empty, there are no messages being
# filtered by the federate. When there are events in the queue it
# indicates the filter federate has messages it is holding onto
# that it needs to forward on (at the indicated time).
#
eq = []
# sub = h.helicsFederateRegisterSubscription(fed, "Charger/EV1_voltage", "")
logger.info('Attempting to enter execution mode')
h.helicsFederateEnterExecutingMode(fed)
logger.info('Entered HELICS execution mode')
hours = 24 * 7 # one week
total_interval = int(60 * 60 * hours)
# Blocking call for a time request at max simulation time
# fake_max_time = int(h.HELICS_TIME_MAXTIME / 1000)
fake_max_time = 100000000
starttime = fake_max_time
#starttime = 0
logger.debug(f'Requesting initial time {starttime}')
grantedtime = h.helicsFederateRequestTime(fed, starttime)
logger.debug(f'Granted time {grantedtime}')
while grantedtime < total_interval:
# value = h.helicsInputGetString(sub)
# logger.debug(f'Got message {value} from random sub at time {grantedtime}.')
# In HELICS, when multiple messages arrive at an endpoint they
# queue up and are popped off one-by-one with the
# "helicsEndpointHasMessage" API call. When that API doesn't
# return a message, you've processed them all.
while h.helicsEndpointHasMessage(endid):
msg = h.helicsEndpointGetMessage(endid)
# msg_str = h.helicsMessageGetString(msg)
# source = h.helicsMessageGetOriginalSource(msg)
dest = h.helicsMessageGetOriginalDestination(msg)
# time = h.helicsMessageGetTime(msg)
# logger.debug(f'Received message from endpoint {source}'
# f' to endpoint {dest}'
# f' for delivery at time {time}'
# f' with payload \"{msg_str}\"')
h.helicsMessageSetDestination(msg, dest)
h.helicsEndpointSendMessage(endid, msg)
# logger.debug(f'Sent message from endpoint {end_name}'
# f' appearing to come from {source}'
# f' to endpoint {dest}'
# f' at time {grantedtime}'
# f' with payload \"{msg_str}\"')
logger.debug(f'Requesting time {fake_max_time}\n')
grantedtime = h.helicsFederateRequestTime(fed, fake_max_time)
logger.debug(f'Granted time {grantedtime}')