validseq,plot=3,bangbang=True,clearall

Author: markgarnold

nintan/cstim.py

@@ -1,6 +1,12 @@
 import concore
 from ast import literal_eval
 import time
+
+try:
+    bangbang = concore.params['bang']
+except:
+    bangbang = False
+
 concore.delay = 0.01
 init_simtime_u = "[0.0, 0.0, 0.0]"
 init_simtime_ym = "[0.0, 0.0, 0.0]"
@@ -12,11 +18,17 @@
 while(concore.simtime<concore.maxtime):
     while concore.unchanged():
         ym = concore.read(1,"ym",init_simtime_ym)
-    try:
-        u[0] = float(literal_eval(input()))
-    except:
-        print("bad input, using .5 instead")
-        u[0] = 0.5        
+    if bangbang:
+        if ym[0] > 0:
+            u[0] = 0.9
+        else:
+            u[0] = 0.1
+    else:
+        try:
+            u[0] = float(literal_eval(input()))
+        except:
+            print("bad input, using .5 instead")
+            u[0] = 0.5        
     print("ym="+str(ym[0])+" u="+str(u[0]));
     concore.write(1,"u",u);
     wallclock2 = time.perf_counter()

nintan/pacq.py

@@ -17,7 +17,7 @@
 try:
     showPlot = concore.params['plot']
 except:
-    showPlot = False
+    showPlot = 0
 
 try:
     tsamp = concore.params['tsamp']
@@ -26,6 +26,11 @@
 except:
     tsamp = 1 
 
+def validseq(t):
+    for i in range(0,len(t)-1):
+        if t[i] >= t[i+1]:
+             print('i='+str(i)+': '+str(t[i])+'>='+str(t[i+1]))
+
 def extract(amplifierData):
     s = 0.0
     for data in amplifierData:
@@ -51,6 +56,16 @@ def readUint16(array, arrayIndex):
     arrayIndex = arrayIndex + 2
     return variable, arrayIndex
 
+def clearall():
+    # Clear TCP data output to ensure no TCP channels are enabled
+    scommand.sendall(b'execute clearalldataoutputs')
+    time.sleep(0.1)
+    # Send TCP commands to set up TCP Data Output Enabled for wide
+    # band of channel A-010
+    scommand.sendall(b'set a-010.tcpdataoutputenabled true')
+    time.sleep(0.1)
+
+
 def acq():
     # Run controller for tsamp second
     scommand.sendall(b'set runmode run')
@@ -127,14 +142,6 @@ def acq():
 # Calculate timestep from sample rate
 timestep = 1 / sampleRate
 
-# Clear TCP data output to ensure no TCP channels are enabled
-scommand.sendall(b'execute clearalldataoutputs')
-time.sleep(0.1)
-
-# Send TCP commands to set up TCP Data Output Enabled for wide
-# band of channel A-010
-scommand.sendall(b'set a-010.tcpdataoutputenabled true')
-time.sleep(0.1)
 
 # Calculations for accurate parsing
 # At 30 kHz with 1 channel, 1 second of wideband waveform data (including magic number, timestamps, and amplifier data) is 181,420 bytes
@@ -148,6 +155,7 @@ def acq():
 waveformBytesPerFrame = 4 + 2
 waveformBytesPerBlock = framesPerBlock * waveformBytesPerFrame + 4
 
+clearall()
 
 concore.delay = 0.01
 init_simtime_u = "[0.0, 0.0, 0.0]"
@@ -169,16 +177,23 @@ def acq():
     oldAmpD = amplifierData
     amplifierTimestamps = []
     amplifierData = []
-    acq_thread.start()
-    nxtacq_thread = threading.Thread(target=acq, daemon=True)
-    ym[0] = extract(oldAmpD)
-    print("ym="+str(ym[0]))
-    if showPlot:
+    if showPlot&1==1:
         plt.plot(oldAmpT,oldAmpD)
-        plt.title("ym="+str(ym[0]))
+        plt.title("ym="+str(extract(oldAmpD))
         plt.xlabel('Time (s)')
         plt.ylabel('Voltage (uV)')
         plt.show()
+    if showPlot&2==2:
+        plt.plot(range(0,len(oldAmpD)),oldAmpD)
+        plt.title("ym="+str(extract(oldAmpD))
+        plt.xlabel('index')
+        plt.ylabel('Voltage (uV)')
+        plt.show()
+    acq_thread.start()
+    nxtacq_thread = threading.Thread(target=acq, daemon=True)
+    ym[0] = extract(oldAmpD)
+    validseq(oldAmpT)
+    print("ym="+str(ym[0]))
     concore.write(1,"ym",ym)
 acq_thread.join()
 print("retry="+str(concore.retrycount))