El código está bien. El tema está en las comparaciones cuando hablamos de tipos en coma flotante, que es asunto de la máquina o del compilador, y pòco puedes hacer desde fuera.
Prueba incluyendo el archivo <math.h>, ahí puede que haya alguna definición o se active alguna cosa para el compilador que pueda ayudarte (aunque, francamente, lo dudo).
Prueba esta porción de código (el mismo que posteaste, pero simplificado), así puedes hacerte idea de lo que está pasando.
Código C:
Ver original#include <stdio.h>
#include <math.h>
int main (void)
{
double LIQST[11],c;
double eps = 1.0e-8;
int i,j,k;
char comp[11][15]={"CO2","H2S","C1","C2","C3","iC4","nC4","iC5","nC 5","C6","C7+"};
do
{
printf("DATOS DEL STOK TANK\nComposicion liquido:\n"); c = 0.0;
for(i=0; i<=10; i++)
{
printf("(%d) %s : ", (i
+1), comp
[i
]); c += LIQST[i];
printf("%24.20f *** %24.20f\n", LIQST
[i
], c
); }
}
while (fabs(c
- 1.0) > eps
); }
(está en C y no C++, no controlo C++). El formato de printf es para que muestre 20 cifras decimales.
El valor eps lo he establecido arbitrariamente para la comparación, si quieres seguir con la doble validación utiliza este sistema, prueba distintos valores. Dependerá de la precisión que necesites.
En cualquier caso, procura evitar las comparaciones de igualdad entre valores en coma flotante.
Si quieres probar también con tipos float, cambia el especificador de formato de scanf de %lf (double, 64 bits) a %f (float, 32 bits). Así podrás ver dónde empiezan a cambiar las cifras decimales.
El Análisis Numérico no es solamente resolver problemas matemáticos mediante un ordenador. Una de sus partes es, precisamente, minimizar los errores por redondeo. Fíjate en el que se te ha presentado con algo tan sencillo como una comprobación de unos pocos datos introducidos por teclado.
Saludos,