Tecnologia
Página 1/6
Informe Final
Algoritmo de Gradientes Conjugados con precondicionamientoPieck, Darío
Introducción
El método de gradientes conjugados permite resolver de forma numérica sistemas de ecuaciones lineales, sin embargo puede fallar si no se cumplen ciertas condiciones, asaber, tener una matriz definida positiva y con bajo número de condición. Sin embargo si el número de condición K -definido como el cociente entre el eigen valor mas bajo y el mas alto- es grande perola matriz es del tipo diagonal dominante, entonces, el precondicionamiento automático permite llegar a una buena aproximación al mismo tiempo que reduce el número de iteraciones necesarias.
Elalgoritmo:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 /* Programa que hace el metodo de los gradientes conjugados con precondicionamiento automatico */ /* Autor:PIECK, Dario */ /* Fecha: 27/4/2007 */
Habilito ciertas librerias de funciones que necesito para correr el programa
#include #include #include // Declaro las funciones que voy a utilizar voidPMV(double*,double*,double*,int); void PVV(double*,double*,double*,int); void PVV_vectorial(double*,double*,double*,int); int main() { size_t n; double* matriz; double* semilla; double* B; double*Ax; double* alfad1; double* X1; double* Ax1; double* DK; double* RK1; double* aux2; double* XK;
Declaración de funciones que se utilizarán durante la ejecución, éstas tienen por objeto simplificarla escritura del programa, como asi tambien su lectura por el programador. Se escribe su nombre y entre paréntesis los valores que ingresa y devuelve. La notación * indica que está trabajando conpunteros.
Aquí ingreso la función principal, la main, esta contiene el bloque principal de instrucciones del programa Dentro de la función main lo primero que hacemos es definir las variables que se...
Regístrate para leer el documento completo.