Abort pri NAN nebo INF

[Jiri Dobes]

Diky moc za rady!!!

dosud jsem to mel cestecne vyresene

if(!finite(residual_norm)) abort();

ale to kontrolovalo NAN jen jednou za cely cyklus. 
Pomoci fenv.h (soucast glibc 2.1) se to da vyresit (odkaz 
od Jirho Masika). Ted jsem nemel cas to zkouset, ale 
pristi tyden... :-)

		Jirka


-- 
*****************************************************************
** Jiri Dobes <dobes na centrum.cz> http://www.bigfoot.com/~bobes 
** Mobil: ++32 (0) 48 59 21 859, ++420 (0) 728 14 32 60 
** Von Karman Institute for Fluid Dynamics   
** Chaussee de Waterloo 72, B-1640 Rhode-St-Genese
*****************************************************************