Skip to content

Instantly share code, notes, and snippets.

@jlblancoc
Created June 7, 2018 10:36
Show Gist options
  • Save jlblancoc/5c933c83a42933d012f07fdd4512a819 to your computer and use it in GitHub Desktop.
Save jlblancoc/5c933c83a42933d012f07fdd4512a819 to your computer and use it in GitHub Desktop.
Ayuda a depuración de cálculo matricial de estructuras en C++
% Comprobación en MATLAB de que una matriz 'K.txt' es correcta:
% Dentro de MATLAB, irse al directorio donde esté el fichero K.txt:
K=load('K.txt');
% "K" debe ser simétrica, por lo que K menos su transpuesta debe ser todo ceros.
errores = K - K';
% Inspeccionar la matriz "errores" en el visor de variables, y en caso de haber
% elementos distintos de cero (ignorad los errores numéricos, valores muy pequeños)
% revisar el ensamblaje en las coordenadas donde ocurra el error.
// Estas anotaciones están pensadas para depurar errores en la práctica propuesta en:
// https://github.com/ingmec-ual/practicas-calculo-matricial-estructuras
// Añadir el siguiente include:
#include <fstream>
// ...
// Después de tener la matriz K ensamblada,
// se puede guardar a un fichero de texto con:
std::ofstream("K.txt") << K;
// Si no se añade un directorio, se guarda en el "directorio de trabajo"
// del ejecutable, que se puede ver haciendo botón derecho en el
// nombre del proyecto (barra lateral de Visual Studio) -> "Abrir carpeta".
// Después, vamos a MATLAB para comprobar que la matriz es correcta.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment