tengo la linea memset(C, 0, N*M); , donde C es una matriz de dobles. (Estoy usando WSL.)
Sin embargo, si miro en gdb, ciertos elementos en la matriz se establecen en -nan(0xffffffffffff8) y otros se establecen en, por ejemplo 9.0096750001652956e-314 .
El primero no da ningún error, pero un += no parece cambiar nada (o al menos no parece hacer desaparecer la cosa nan), mientras que el segundo es un problema si el elemento no se cambia. o solo tiene += 0 , ya que la comparación if (0 == C[i][j]) luego falla.
Si establezco los valores en 0 manualmente, estos problemas no surgen en absoluto.
¿Es esto una cosa de WSL, o hay algo sobre memset que no entiendo?
No inicializa completamente la matriz: memset() espera una cantidad de bytes. Suponiendo que el diseño de la matriz sea lineal, ya sea 1D o 2D, debe borrar el sizeof(double) * N * M bytes.
Si su matriz se define como una matriz 2D, podría escribir:
#define N 10 #define M 20 double C[N][M]; memset(C, 0, sizeof C);Si la matriz se recibe como un argumento de función, en realidad obtiene un puntero, por lo que debe tener más cuidado:
void clear_matrix(double C[N][M]) { memset(C, 0, sizeof(*C) * N); }o posiblemente más legible:
void clear_matrix(double C[N][M]) { memset(C, 0, sizeof(C[0][0]) * N * M); }o simplemente, como sugiere Lundin, pero es probable que se rompa si cambia el tipo de elemento de matriz:
void clear_matrix(double C[N][M]) { memset(C, 0, sizeof(double[N][M]); } Sin embargo, tenga en cuenta que memset() borrará los datos de la matriz a todos los bits cero, lo que establece los valores double en +0.0 si el sistema usa la representación IEEE-754, pero no es completamente portátil. Una versión portátil usaría un bucle anidado y un buen compilador generará la misma llamada de memset o código en línea si es apropiado para el sistema de destino:
#include <stddef.h> #define N 10 #define M 20 void clear_matrix(double C[N][M]) { for (size_t i = 0; i < N; i++) { for (size_t j = 0; j < M; j++) { C[i][j] = 0.0; } } }El M*N es el número de elementos (dobles) pero memset espera el número de bytes . Para una matriz double C[M][N] , debe hacer memset(C, 0, sizeof(C)); .