Skip to content

Commit 73df17b

Browse files
committed
Round branching numbers/fractions to 6 significant digits
1 parent fe62d54 commit 73df17b

1 file changed

Lines changed: 5 additions & 5 deletions

File tree

src/AutomaticMechanismAnalysis.jl

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -389,14 +389,14 @@ Generate a string report from the analysis objects
389389
function getrxnanalysisstring(sim,ra;branchingcutoff=1e-2,radbranchfract=0.01)
390390
spcname = sim.names[ra.spcind]
391391
rstr = getrxnstr(sim.domain.phase.reactions[ra.rxnind])
392-
sens = ra.sens
392+
sens = round(ra.sens,sigdigits=6)
393393
s = "Analyzing $spcname sensitivity to $rstr at a value of $sens \n"
394394
s *= "\n"
395395
for branch in ra.branchings
396396
sname = sim.names[branch.spcind]
397397
s *= "Key branching for $sname \n"
398398
for i = 1:length(branch.rxninds)+1
399-
br = branch.branchingratios[i]
399+
br = round(branch.branchingratios[i],sigdigits=6)
400400
if br < branchingcutoff
401401
break
402402
else
@@ -414,22 +414,22 @@ function getrxnanalysisstring(sim,ra;branchingcutoff=1e-2,radbranchfract=0.01)
414414
s *= "Associated key reaction path in $spname loss direction \n"
415415
for i = 1:length(rp.rxninds)
416416
rstr = getrxnstr(sim.reactions[rp.rxninds[i]])
417-
br = rp.branchfracts[i]
417+
br = round(rp.branchfracts[i],sigdigits=6)
418418
s *= "$rstr at path branching of $br \n"
419419
end
420420
else
421421
s *= "Associated key reaction path in $spname production direction \n"
422422
revinds = reverse(rp.rxninds)
423423
for i = 1:length(rp.rxninds)
424424
rstr = getrxnstr(sim.reactions[revinds[i]])
425-
br = rp.branchfracts[i]
425+
br = round(rp.branchfracts[i],sigdigits=6)
426426
s *= "$rstr at path step branching of $br \n"
427427
end
428428
end
429429
s *= "\n"
430430
end
431431
if abs(ra.radprodlossfract) > radbranchfract
432-
radfract = abs(ra.radprodlossfract)
432+
radfract = round(abs(ra.radprodlossfract),sigdigits=6)
433433
if ra.radprodlossfract > 0
434434
s *= "Reaction accounts for $radfract of the net radical production \n"
435435
else

0 commit comments

Comments
 (0)