Resolución de Permanente
[Métodos de resolución]

Diagrama de colaboración para Resolución de Permanente:

Funciones
int	resolucion_permanente (basicos *p_DATOS)
	Función que se encarga de la resolución del sistema permanente del circuito a analizar.

---

Descripción detallada

Grupo encargado de la resolución del sistema permanente del circuito que se esté analizando

---

Documentación de las funciones

int resolucion_permanente

(

basicos *

p_DATOS

)

Función que se encarga de la resolución del sistema permanente del circuito a analizar.

Parámetros:

Puntero

a estructura basicos

Definición en la línea 135 del archivo estructuras.cpp.

{
Vec X,W; Mat T;
PetscErrorCode ierr;
KSP ksp; 
PetscViewer viewer, escribeDATOS;
//Creamos el Vector que guardará la solución del sistema estacionario
p_DATOS->x_perm=p_DATOS->creax_perm();
VecAssemblyBegin(p_DATOS->x_perm); VecAssemblyEnd(p_DATOS->x_perm);
//Cargamos los parametros de preacondicionamiento
PetscOptionsInsertFile("opciones.txt");
ierr = KSPCreate(PETSC_COMM_WORLD,&ksp);CHKERRQ(ierr);
//Copias de los elementos con los que trabajaremos T,W y X
T=p_DATOS->T; W=p_DATOS->W; X=p_DATOS->x_perm;
//Resolución del sistema lineal
ierr = KSPSetOperators(ksp,T,T,DIFFERENT_NONZERO_PATTERN);CHKERRQ(ierr);
KSPSetFromOptions(ksp);
ierr = KSPSolve(ksp,W,X);CHKERRQ(ierr);
//llamada a funcion que excriba los valores del sistema
PetscViewerASCIIOpen(PETSC_COMM_WORLD, "Datos_x_perm.data", &escribeDATOS);
MatView(p_DATOS->T, escribeDATOS);
VecView(p_DATOS->W, escribeDATOS);
//Volcado del resultado en la Estructura DATOS BASICOS
p_DATOS->x_perm=X;
ierr = MatDestroy(T);CHKERRQ(ierr);
ierr = VecDestroy(W);CHKERRQ(ierr);
ierr = VecDestroy(X);CHKERRQ(ierr);

//Escritura de la solución en fichero
PetscViewerASCIIOpen(PETSC_COMM_WORLD, "x_perm.data", &viewer);
VecView(p_DATOS->x_perm, viewer);
return 0;
}