@@ -876,70 +876,4 @@ def _default_time_vector(sys, N=None, tfinal=None):
876876 if N is None :
877877 N = int (np .clip (tfinal / ideal_dt , N_min_ct , N_max )) # N<-[N_min, N_max]
878878
879- return np .linspace (0 , tfinal , N , endpoint = False )
880-
881- def junk (sys , N , tfinal ):
882- N_max = 5000
883- N_min = 10
884- if tfinal is None :
885- ideal_tfinal , ideal_dt = _ideal_tfinal_and_dt (sys )
886- N = int (np .clip (ideal_tfinal / ideal_dt , N_min , N_max ))# N->[N_min, N_max]
887- if isdtime (sys , strict = True ):
888- # for discrete time, change tfinal if N gets too large/small
889- tfinal = sys .dt * N
890- else :
891- # for continuous time, simulate to tfinal but change dt if N large
892- tfinal = ideal_tfinal
893- else :
894- if isdtime (sys , strict = True ):
895- N = int (tfinal / sys .dt )
896- else :
897- _ , ideal_dt = _ideal_tfinal_and_dt (sys )
898- N = int (np .clip (tfinal / sys .dt , N_min , N_max )) # N->[N_min, N_max]
899-
900-
901- if isdtime (sys , strict = True ):
902- if tfinal is None :
903- tfinal , dt = _ideal_tfinal_and_dt (sys )
904- # sanity bounds on tfinal
905- N = int (np .clip (tfinal / sys .dt , N_min , N_max )) # [N_min, N_max]
906- tfinal = sys .dt * N
907- else :
908- N = int (tfinal / sys .dt )
909- else :
910- if tfinal is None :
911- tfinal , dt = _ideal_tfinal_and_dt (sys )
912- N = int (np .clip (tfinal / dt , N_min , N_max )) # [N_min, N_max]
913- dt = tfinal / N
914- return np .linspace (0 , tfinal , N , endpoint = False )
915-
916- if 1 :
917- # run to N steps, or to ideal tfinal if N not given
918- if N is None :
919- # run to tfinal
920- if tfinal is None :
921- tfinal , dt = _ideal_tfinal_and_dt (sys )
922- N = int (np .clip (tfinal / sys .dt , N_min , N_max )) # [N_min, N_max]
923- tfinal = sys .dt * N
924- else :
925- tfinal = sys .dt * N
926- else :
927- # continuous-time
928- ideal_tfinal , ideal_dt = _ideal_tfinal_and_dt (sys )
929-
930-
931- # tfinal
932- if tfinal is None :
933- ideal_tfinal , ideal_dt = _ideal_tfinal_and_dt (sys )
934-
935- if N is None :
936- ideal_tfinal , ideal_dt = _ideal_tfinal_and_dt (sys )
937- if tfinal is None :
938- tfinal = ideal_tfinal
939- N = int (np .clip (tfinal / dt , N_min , N_max )) # limit to [N_min, N_max]
940- if isdtime (sys , strict = True ):
941- tfinal = sys .dt * N
942- else :
943- dt = tfinal / N
944- return np .arange (0 , tfinal , dt )
945- return np .linspace (0 , tfinal , N , endpoint = False )
879+ return np .linspace (0 , tfinal , N , endpoint = False )
0 commit comments