#vectorMathLibrary_2023FJun30_py.py
#by 	David at ColeCanada.com
#date: 	2021FJun30
#
#purpose: Library of functions to process a wave stored as a vector of
#   positive samples with or without timestamps
#
#
import math
def makeVectorVIX(nRange,type = "square"):
# Purpose: generates VIX waveforms in a complex list vector
# By: David@ColeCanada.com
# Date: 2023FJun29
# Tested: 2023FJun29
# NOTE: If the number of samples is N,
#       the list of samples will not include a sample at time 0.
#       eg. a sample #1 of a sawtooth (max value of 1) will NOT be 0,
#           it will be 1/N
#
# used to generate a list of samples for 1 cycle of a Wave in a Vector
# valid types are "toSin90", "toSin360","complex-toSin90", "complex-toSin360",
#               "square", "unary", "complex-square" and "complex-unary"
#               "sawtooth" and "complex-sawtooth".
# Version 1 added "sawtooth" waveform
#
# eg  VIXdegListA=makeVectorVIX(360,)
# produces [0,1,2,3,4,5,   . . . . 359]
    #print("waveformIn:",type)
    #if type=="toSin90":print ("type is:",type)
    typeList=[ \
                "sawtooth","complex-sawtooth", \
                "square","complex-square", \
                "toSin360","complex-toSin360", \
                "toSin90","complex-toSin90", \
                "unary","complex-unary"]
    if typeList.count(type)==0:
        print("Err001 wave type not in:",typeList)
        print("... set to unary waveform")
    #rType = type(nRange)
    #print("rType:",rType)
    #if type(nRange) == 'int' :
    #    print("Err003 range must be integer")
    #    print("  used 8")
    #    nRange = 8
    isRangeOK=True
    if nRange == 0 or nRange > 10000 : isRangeOK=False
    dohrange=1.0 * nRange
    if dohrange < 1.0 : isRangeOK=False
    intRange=1
    # intRange will be assigned next
    if isRangeOK==False:
        print ("Err002 nr line 51, range is invalid or exceeds 10000")
        print("  using nRange:8")
        #nRange = int(8)
        intRange= int(8)
    else:
        #print("at 56 nRange:",nRange)
        intRange=nRange
    #print("intRange:",intRange)
    thisList=[]     #empty list
    dX=1./(intRange)
    
    #print("in:",thisList)
    for iterN in range(intRange):
        iterN=iterN+1
        #print("degN:",degN)
        if type=="toSin90" or type=="complex-toSin90":
            maxDegree=90
            #print("set_maxDegree:",maxDegree)
            degrees=iterN*maxDegree/intRange
            rads=(degrees)*2*math.pi/360
            #print ("after for, type is:",type, "rads is",rads)
        #end if
        if type=="toSin360" or type=="complex-toSin360":
            maxDegree=360
            #print("set_maxDegree:",maxDegree)
            degrees=iterN*maxDegree/intRange
            rads=(degrees)*2*math.pi/360
        #end if

        #thisList.append(sin(rad)
        value=1.0
        if type=="toSin90" or type=="complex-toSin90":
            value=math.sin(rads)
        if type=="toSin360" or type=="complex-toSin360":
            value=math.sin(rads)
        if type=="unary" or type=="complex-unary":
            value=1.0
        if type=="sawtooth" or type=="complex-sawtooth":
            value=iterN*dX
        #end if
        if type=="square" or type=="complex-square":
            if iterN<intRange/2:
                value=1
            else:
                value=0
            #end if
        #end if
            
        isComplex=False
        if type=="complex-unary" or type=="complex-square" :
            isComplex=True
        if type=="complex-toSin90" or type=="complex-toSin360":
            isComplex=True
        if type=="complex-sawtooth":
            isComplex=True
        if isComplex :
            realN= (value*1000)//1
            realN=realN/1000
            dT=(dX*iterN*1000)//1
            dT=dT/1000
            itemN=(realN+dT*1j)
            #print("itemN:",itemN)
        else:
            realN= (value*1000)//1
            itemN=realN/1000
            #print("itemN",itemN)
        #end if
        thisList.append(itemN)
    #end for
    return thisList
# def makeVectorVIX.py

def makeJSON(dt,wholeDays,wholeSeconds,fracSeconds,sRate, hertzEstimate,
#              Date       hh:mm:ss      4568Hz   100r/s       12KHz   
#             (secs)       (secs)       (secs)  readings/sec (secs)
#          65743402        36899       .00067                .00005/sec
unixTimeStart="0",
numSamp=60,     #to fix  ................................................
sampMax=1,
sampBase=0x10000,
equivUnits="3v3",
name="<no-name>.tix",
description="boiler #45 pressure",
auth="D@CC",
gps="45.8/-78.0/1000.1",
):
#                  variable name    extended name     explanation/source
#                  -------------    --------------    ------------------------------
    print("01 dt:",dt)              #time interval    from name
    print("08 nS:",vReadCount)      #vectorCount      number of readings
    print("09 sM:",sampMax)         #maxRead          max. value read
    print("05 sB:",sampBase)        #fmtRead          special format if any
    print("06 hz:",hertzEstimate)   #vectorHz         Hz (cps)
    print("07 uT:",vUnixTime)       #unixDateTime     unix time (# of sec since ue)
    print("02 wD:",intDaysInSec)    #from unixDate    int. days
    print("03 wS:",intSec)          #from unixDate    int. secs since midnight
    print("04 fS:",fracSec)         #from unixDate    fractional secs
    print("10 sR:",npsReads)        #readsPerSec      number of readings/sec
    print("11 eU:",vUnits)          #vectorUnits      units (e.g. volts)
    print("12 na:",vName)           #vectorName       actual file name w extension
    print("13 de:",vDescr)          #vectorDesc       vector description
    print("14 au:",vAuthor)         #vectorAuthor     vector author (person)
    print("15 gp:",vGPSA)            #vectorGPSA        vector coordinates (lat,long,alt)
    print("16 wp:",vPeriod          #vectorPeriod     vector period (1/Hz)
    
#hertzEstimate assuming 1 cycle per full vector
#sampMax:      assuming 1 = sampMax when sampBase is 0x10000
#sampBase      is 0x10000 assuming ADC is a RPi Pico
#equivUnits    assuming sampMax is 1 is 3.3 Vcc which is 3v3   
    return dt,unixTimeStart,dt/numSamp,sampMax,equivUnits,  name,  description, auth, gps
#      65743402.234598      .001   2.2   1.0volts  Sampset37A, boiler #45 pressure
# derived info: period, timeTrigger, nSamples, sampMax ,  sampMin, Hertz
#          .0000001125      6894323     100     3.1        0.5    12kHz
#where . . . 
#unixTimeStart is unixTime (GMT) in microseconds
#dt/samp       is time between readings
#sampMax       is actual maximum reading (as in vector)
#equivUnits    is multiplier eg 1V  is 1 volt
#name          is *.vix file name of this data vector, *.tix is its timestamp
#description   is phrase describing this set of readings
#auth          is operator/author/engineer/technician eg D@CC
#GPS           is latitude/longitude/altitude(m)
# test program follows:
nSamp= 60
#print(type(nSamp))
waveform="complex-sawtooth"
myVector=makeVectorVIX(nSamp,waveform) # for "sin" nSamp sb 90 or 360 etc
print(myVector)
code='myVector=makeVectorVIX('+str(nSamp)+',"'+waveform+'")'
print(code)
#newVector=myVector*8
print("type(myVector):",type(myVector[1]))
myJSON=makeJSON(0.01,65743402.234598,1000,0,6000,0.001,65743402.2,60,1,0x10000,"1.0V","sampSet37A10ms.vix", "boiler #45 pressure", \
                 "D@CC","45.33857665474693,-75.52479136378064,35m")

#
#  The following calculations need more work and JSON format !!!!!!!!!!!!!!!!!!!!
#
print("myJSON:",myJSON)
myJSONfileName="sampSet37A10ms.json" # which is ["sampSet37A10ms.vix","sampSet37A10ms.tix"]
print("myJSONfileName: ",myJSONfileName,)
#vectorMathLibrary_2021FJun09

#vectorMathLibrary stored in 203/vectorMathLibrary_2023FJun30_py.txt

#vectorMathLibrary_2021FJun09.py
#vectorMathLibrary_2023FJun30_py.py