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

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

Funciones
Vec	creavectorauxiliar (basicos *p_DATOS)
	Función que se encarga crear un vector auxiliar de cálculo del sistema a resolver.
void	num_iteraciones (basicos *p_DATOS)
	Función que se encarga calcular el número de iteraciones necesarias en el método iterativo directo utilizado.
int	resolucion_transitorio (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 transitorio del circuito que se esté analizando

---

Documentación de las funciones

Vec creavectorauxiliar

(

basicos *

p_DATOS

)

Función que se encarga crear un vector auxiliar de cálculo del sistema a resolver.

Parámetros:

Puntero

a estructura basicos return Vector tipo PETSc

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

{
Vec z;
VecCreate(PETSC_COMM_SELF,&z);
VecSetSizes(z,PETSC_DECIDE,p_DATOS->filas());
VecSetFromOptions(z);
VecAssemblyBegin(z); VecAssemblyEnd(z); 
return z;
}

void num_iteraciones

(

basicos *

p_DATOS

)

Función que se encarga calcular el número de iteraciones necesarias en el método iterativo directo utilizado.

Parámetros:

Puntero

a estructura basicos

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

{
p_DATOS->numv=(p_DATOS->tf-p_DATOS->t0)/p_DATOS->h;

}

int resolucion_transitorio

(

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 177 del archivo estructuras.cpp.

{
Vec BBCC=p_DATOS->BBCC;
KSP ksp;
PetscOptionsInsertFile("opciones.txt");
KSPCreate(PETSC_COMM_WORLD,&ksp);
KSPSetOperators(ksp,p_DATOS->G,p_DATOS->G,DIFFERENT_NONZERO_PATTERN);
KSPSetFromOptions(ksp);         
//llamada a funcion que excriba los valores del sistema
PetscViewer escribeDATOStran;
PetscViewerASCIIOpen(PETSC_COMM_WORLD, "Datos_x_trans.data", &escribeDATOStran);
MatView(p_DATOS->G, escribeDATOStran);
VecView(BBCC, escribeDATOStran);                
KSPSolve(ksp,BBCC,p_DATOS->x_trans);//KSPSolve(KSP ksp,Vec b,Vec x)
PetscViewer viewer;
PetscViewerSetFormat(viewer,PETSC_VIEWER_ASCII_COMMON );                
PetscViewerASCIIOpen(PETSC_COMM_WORLD, "x_trans.data", &viewer);
VecView(p_DATOS->x_trans, viewer);
p_DATOS->Wn=p_DATOS->creaWn();
p_DATOS->Wn1=p_DATOS->creaWn1();
return 0;
}