fix bug in the betar

A-006 · A-006 · commit bbe292635095 · 2025-07-12T11:56:19.000+08:00
diff --git a/source/source_cell/pseudo.cpp b/source/source_cell/pseudo.cpp
@@ -16,7 +16,8 @@ void pseudo::check_betar()
 		int flag = mesh+1;
 		for (int ir = 0; ir < mesh; ir++)
 		{
-			if (std::abs(betar(ib, ir)) < 1.0e-30)
+			// to check is non-normal number and it shoule not be zero
+			if ((std::abs(betar(ib, ir)) < 1.0e-30) && (std::abs(betar(ib, ir))!=0))
 			{
 				flag = ir;
 				break;

Original file line number	Diff line number	Diff line change
`@@ -16,7 +16,8 @@ void pseudo::check_betar()`
`16`	`16`	`int flag = mesh+1;`
`17`	`17`	`for (int ir = 0; ir < mesh; ir++)`
`18`	`18`	`{`
`19`		`- if (std::abs(betar(ib, ir)) < 1.0e-30)`
	`19`	`+ // to check is non-normal number and it shoule not be zero`
	`20`	`+ if ((std::abs(betar(ib, ir)) < 1.0e-30) && (std::abs(betar(ib, ir))!=0))`
`20`	`21`	`{`
`21`	`22`	`flag = ir;`
`22`	`23`	`break;`