Esempio n. 1
0
    def checkEvents(self, maxevtpts, eventActive, eventDir, eventTerm,
                    eventInt, eventDelay, eventTol, maxbisect, eventNearCoef):

        if not self.initBasic:
            raise InitError('You must initialize the integrator before ' + \
                            'checking events. (initBasic)')

        verify_nonneg('maxevtpts', maxevtpts, _all_int)

        list_len_val = (self.nEvents, 'nEvents')

        verify_intbool('eventActive', eventActive,
                      list_ok=True, list_len=list_len_val)

        verify_values('eventDir', eventDir, [-1, 0, 1], list_ok=True,
                      list_len=list_len_val)

        verify_intbool('eventTerm', eventTerm,
                       list_ok=True, list_len=list_len_val)

        verify_nonneg('eventInt', eventInt, _real_types,
                      list_ok=True, list_len=list_len_val)

        verify_nonneg('eventDelay', eventDelay, _real_types,
                      list_ok=True, list_len=(self.nEvents, 'nEvents'))

        verify_pos('eventTol', eventTol, _real_types,
                   list_ok=True, list_len=list_len_val)

        verify_nonneg('maxbisect', maxbisect, _all_int,
                      list_ok=True, list_len=list_len_val)

        verify_pos('eventNearCoef', eventNearCoef, _real_types)
Esempio n. 2
0
 def checkBasic(self, rhs, phaseDim, paramDim, nAux, nEvents, nExtInputs,
                hasJac, hasJacP, hasMass, extraSpace):
     # Check that inputs to this function are correct
     try:
         if not isinstance(rhs, str):
             raise TypeError("right hand side rhs must be a string")
         verify_nonneg('phaseDim', phaseDim, _all_int)
         verify_nonneg('paramDim', paramDim, _all_int)
         verify_intbool('hasJac', hasJac)
         verify_intbool('hasJacP', hasJacP)
         verify_intbool('hasMass', hasMass)
         verify_nonneg('nAux', nAux, _all_int)
         verify_nonneg('nEvents', nEvents, _all_int)
         verify_nonneg('nExtInputs', nExtInputs, _all_int)
         verify_nonneg('extraSpace', extraSpace, _all_int)
     except:
         print sys.exc_info()[0], sys.exc_info()[1]
         raise InitError('Integrator initialization failed!')