Example #1
0
    def __init__(self,coeffFileName):

        Integrator.__init__(self,coeffFileName)
        i=self.initVars
        i.yinit=np.array([i.yval])
        i.nVars=len(i.yinit)
Example #2
0
try:
    from lab5_funs import Integrator
except ImportError:
    import os,sys
    libdir=os.path.abspath('../')
    sys.path.append(libdir)
    from lab5_funs import Integrator

if __name__=="__main__":
    import numpy as np
    import scipy as sp
    import matplotlib.pyplot as plt
 
    theSolver=Integrator('daisy.ini')
    timeVals,yVals,errorList=theSolver.timeloop5Err()
    whiteDaisies=[frac[0] for frac in yVals]
    
    thefig=plt.figure(1)
    thefig.clf()
    theAx=thefig.add_subplot(111)

    points=theAx.plot(timeVals,whiteDaisies,'b+')
    points[0].set_markersize(12)
    theLines=theAx.plot(timeVals,yVals)
    theLines[0].set_marker('+')
    theLines[1].set_linestyle('--')
    theLines[1].set_color('k')
    theAx.set_title('lab 5 interactive 6')
    theAx.set_xlabel('time')
    theAx.set_ylabel('fractional coverage')
    theAx.legend(theLines,('white daisies','black daisies'),loc='best')
Example #3
0
    def __init__(self, coeffFileName):

        Integrator.__init__(self, coeffFileName)
        i = self.initVars
        i.yinit = np.array([i.yval])
        i.nVars = len(i.yinit)
Example #4
0
try:
    from lab5_funs import Integrator
except ImportError:
    import os,sys
    libdir=os.path.abspath('../')
    sys.path.append(libdir)
    from lab5_funs import Integrator

if __name__=="__main__":
    import numpy as np
    import scipy as sp
    import matplotlib.pyplot as plt
 
    theSolver=Integrator('daisy.ini')
    timeVals,yVals,errorList = theSolver.timeloop5Err()

    thefig=plt.figure(1)
    thefig.clf()
    theAx=thefig.add_subplot(111)
    theLines=theAx.plot(timeVals,yVals)
    theLines[1].set_linestyle('--')
    theLines[1].set_color('k')
    theAx.set_title('lab 5 interactive 4')
    theAx.set_xlabel('time')
    theAx.set_ylabel('fractional coverage')
    theAx.legend(theLines,('white daisies','black daisies'),loc='center right')
    plt.show()