There seems to be a bug in bicgstab.cpp in meep-1.2.
Without the following patch physical.cpp test gives nan ratio, although
make check passed.
Best regards
Munehiro
--- bicgstab.cpp.orig 2012-12-18 22:52:19.159110261 +0900
+++ bicgstab.cpp 2012-12-18 22:53:44.655994596 +0900
@@ -84,7 +84,7 @@
long double sum = 0;
for (i = 0; i < n; ++i) {
double xabs = fabs(x[i]);
- if (xmax > xabs) xmax = xabs;
+ if (xmax < xabs) xmax = xabs;
}
xmax = max_to_all(xmax);
if (xmax == 0) return 0;
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
