77
88def countCrossings ( plot , thresh ):
99 vec = moose .element ( plot ).vector
10- # print( vec )
10+ print ( vec )
1111 # print (vec[:-1] <= thresh)
1212 return sum ( (vec [:- 1 ] < thresh ) * (vec [1 :] >= thresh ) )
1313
1414def test_xreac2 ():
1515 mfile = os .path .join ( os .path .dirname ( __file__ ), 'OSC_diff_vols.g' )
1616 runtime = 4000.0
1717 modelId = moose .loadModel ( mfile , 'model' , 'ee' )
18- print (modelId )
19-
2018 kin = moose .element ( '/model/kinetics' )
2119 compt1 = moose .element ( '/model/compartment_1' )
2220 compt1 .x1 += kin .x1
@@ -29,26 +27,8 @@ def test_xreac2():
2927 ds1 = moose .Dsolve ( '/model/kinetics/dsolve' )
3028 s1 = moose .Stoich ( '/model/kinetics/stoich' )
3129 s1 .compartment = moose .element ( '/model/kinetics' )
32- print (s1 .compartment )
3330 s1 .ksolve = ks1
34- print (s1 .ksolve )
3531 s1 .dsolve = ds1
36- print (s1 .dsolve )
37-
38-
39- for r in moose .wildcardFind ('/model/##[TYPE=Reac]' ):
40- print ('REAC ====' )
41- print (r )
42- print ('subs' , r .neighbors ['sub' ])
43- print ('prds' , r .neighbors ['prd' ])
44-
45- # ps = moose.wildcardFind('/model/##')
46- # print(len(ps), "Total paths")
47- # for p in sorted(ps, key=lambda k:k.name):
48- # print(p)
49- # quit()
50-
51- print ('s1 path: ' , s1 .path )
5232 s1 .path = '/model/kinetics/##'
5333
5434 ks2 = moose .Ksolve ( '/model/compartment_1/ksolve' )
@@ -57,13 +37,10 @@ def test_xreac2():
5737 s2 .compartment = moose .element ( '/model/compartment_1' )
5838 s2 .ksolve = ks2
5939 s2 .dsolve = ds2
60-
61-
6240 s2 .path = '/model/compartment_1/##'
6341
6442 ds2 .buildMeshJunctions ( ds1 )
6543
66-
6744 moose .reinit ()
6845 moose .start ( runtime )
6946 # I don't have an analytic way to assess oscillations
0 commit comments