1) Algunos de ustedes han manifestado que el algoritmo de EG les está
demorando mucho al momento de aplicarlo a la matriz B de la NASA. Esto
va a depender de la computadora que utilicen y efectivamente podría
demorar bastante o no (en mi computadora la EG junto con sustitución
hacia atrás ejecuta en 1 minuto). En todo caso recuerden que no es
obligación que implementen todos los algoritmos y pueden utilizar los
provistos por Octave/Python; que seguro serán más rápidos por estar
optimizados.
2) Tengan en cuenta que para la matriz C solamente les pedimos que
resuelvan con ordenamiento y Cholesky (no es obligatorio que apliquen EG
ni métodos iterativos).
3) Cuando se pide que reporten una matriz lo que esperamos es que nos
muestren la estructura de la matriz. Por ejemplo con un "gráfico" del
tipo que se obtiene con el comando spy() de Octave. Cuando pedimos que
reporten solución, no pretendemos que nos den los valores cuantitativos
de cada coordenada sino una medida de la distancia entre su solución
estimada y la solución "verdadera" alfa. Pueden tomar como solución
"verdadera" la que devuelve Octave/Python. En la parte 1, donde el
sistema es pequeño, deben reportar además los valores cuantitativos de
matrices y vector solución.
4) Pueden utilizar reordenamiento de Cuthill inverso, directo o ambos y
comparar.
5) Recuerden que para los algoritmos iterativos tipo Jacobi y
Gauss-Seidel les pedimos que incluyan en su informe un gráfico con la
evolución del error absoluto ||xk-alfa||. También pueden graficar la
evolución del residuo ||A*xk-b||.
En respuesta a Matías Valdés
Re: Algunas aclaraciones respecto al Obligatorio 1
de Diego Kiedanski -
Con respecto a la parte 2, a la hora de reportar la solución, aparte de lo que dijo Matías antes, hagan lo siguiente:
Escriban en el informe las primeras 5 coordenadas de la solución:
Por ejemplo: si la solución es un vector de 500 coordenadas:
[1, 2,4, 17, -13.06, 4, .... , 70, 11, 4]
En el informe pongan: las primeras 5 coordenadas de la solución son:
[1, 2, 4, 17, -13.06]