Autor |
Mensaje |
Pilar Lopez
Nivel 3
Registrado: 10 Nov 2005
Mensajes: 30
Ubicación: Ituzaingo
Carrera: Civil
|
|
Hola:
Alguien puede darnos una pista para leer las matrices tridiagonales por bloques del TP con pascal? Nuestra idea es escribir las matrices en un txt. No vemos si la estructura de las matrices se puede desarrollar en forma general (para A6 o A9 o A12).
Gracias, saludos.
Pilar
|
|
|
|
|
|
|
|
|
COKI
Nivel 9
Edad: 45
Registrado: 17 Ene 2006
Mensajes: 2044
Ubicación: Coghlan/Montevideo
Carrera: Mecánica
|
|
|
|
|
facundo.olano
Nivel 8
Edad: 36
Registrado: 10 Ago 2006
Mensajes: 808
Ubicación: encadenado al ánima
Carrera: Informática
|
|
Pilar Lopez escribió:
|
Hola:
Alguien puede darnos una pista para leer las matrices tridiagonales por bloques del TP con pascal? Nuestra idea es escribir las matrices en un txt. No vemos si la estructura de las matrices se puede desarrollar en forma general (para A6 o A9 o A12).
Gracias, saludos.
Pilar
|
Me parece complicarse demasiado usar txts. Creo que te conviene construir "a mano" las B1, B2, B3 y D1, y armarte una función que te haga las An indicándole n (múltiplo de 3), en base a las otras (Fijate que las An son siempre B1 D1 en als primeras 3 filas, D1 B3 en las últimas 3 y D1 B2 D1 en las del medio, tomadas de a tres).
|
|
|
|
_________________
|
|
|
|
|
Gualicho
Nivel 8
Edad: 35
Registrado: 18 Sep 2007
Mensajes: 715
Ubicación: En el templo de Momo...
Carrera: Informática
|
|
facundo.olano escribió:
|
Pilar Lopez escribió:
|
Hola:
Alguien puede darnos una pista para leer las matrices tridiagonales por bloques del TP con pascal? Nuestra idea es escribir las matrices en un txt. No vemos si la estructura de las matrices se puede desarrollar en forma general (para A6 o A9 o A12).
Gracias, saludos.
Pilar
|
Me parece complicarse demasiado usar txts. Creo que te conviene construir "a mano" las B1, B2, B3 y D1, y armarte una función que te haga las An indicándole n (múltiplo de 3), en base a las otras (Fijate que las An son siempre B1 D1 en als primeras 3 filas, D1 B3 en las últimas 3 y D1 B2 D1 en las del medio, tomadas de a tres).
|
Coincido...
Te dejo una que habia usado para el tp en su momento para que uses de ejemplo. Puede ser que hayan cambiado el orden de donde estaba cada matriz y esas cosas...
Código:
|
void cargarMatrizA9(double[][] matb1,double[][] matb2,double[][] matb3,double[][] matd1){
for (int i=0;i<3;i++){
for (int j=0;j<3;j++){
mat[i][j]=matb1[i][j];
mat[i+3][j]=matd1[i][j];
mat[i][j+3]=matd1[i][j];
mat[i+3][j+3]=matb2[i][j];
mat[i+6][j]=0;
mat[i][j+6]=0;
mat[i+6][j+3]=matd1[i][j];
mat[i+3][j+6]=matd1[i][j];
mat[i+6][j+6]=matb3[i][j];
}
}
}
|
Por si no se entiende, voy recorriendo la matriz a llenar (A9) con dos ciclos en los cuales, dependiendo la posicion en la que me encuentro, le cargo las matrices enteras que ya tengo cargadas (B*,D1)...
Suerte!
|
|
|
|
_________________ "Por eso te pido (amigo desconocido), si ves a mi rock perdido, lo traigas por aqui!"
|
|
|
|
|
Pilar Lopez
Nivel 3
Registrado: 10 Nov 2005
Mensajes: 30
Ubicación: Ituzaingo
Carrera: Civil
|
|
Muchisimas gracias por la Info, voy aprobar aver si puedo resolverlo de esa manera.
La entrega del TP es 22 de Abril ( al horno...).
|
|
|
|
|
|
|
|
|
ale_vans
Nivel 7
Edad: 35
Registrado: 22 May 2008
Mensajes: 304
Ubicación: Vte. Lopez
Carrera: No especificada
|
|
no quiero pincharles el globo pero.... en el tp dice explicitamente:
"Es un requisito que el ó los programas
desarrollados tengan en cuenta la estructura de estas matrices y operen solamente
con los elementos no nulos de las mismas, o sea, que dichas matrices An no deberán
estar representadas computacionalmente por un arreglo bidimensional sino que las
fórmulas de los métodos deberán escribirse para estas matrices en particular"
de ahi se sacan dos conclusiones: la primera es que no hace falta crear una funcion que cargue matrices An, ya que las formulas tienen q hacerse para estas matrices en particular.
Y segundo, no puede usarse un array bidimensional, o sea, al horno.
En clase explicaron una forma de hacerlo usando 3 vectores, y solo cargando los elementos no nulos. en el primer vector se carga la posicion de la fila, en el segundo vector la posicion de la columna y en el tercer vector el valor q guarda esa posicion. Por ultimo queda diseñar una funcion de busqueda que devuelva el valor en una posicion determinada, o cero en caso de que no lo encuentre.
Espero q les sirva, saludos
|
|
|
|
|
|
|
|
|
Pilar Lopez
Nivel 3
Registrado: 10 Nov 2005
Mensajes: 30
Ubicación: Ituzaingo
Carrera: Civil
|
|
Te referis a que dijieron que a Ase puede escribir como A= D-L- U????
|
|
|
|
|
|
|
|
|
ale_vans
Nivel 7
Edad: 35
Registrado: 22 May 2008
Mensajes: 304
Ubicación: Vte. Lopez
Carrera: No especificada
|
|
no, no me refiero a eso y te recomiendo q no uses esa notacion porq es solo para fines teoricos, segun el Burden y los profes.
Me refiero a q no podes simplemente crear un vector de dos dimensiones y cargarlo con los elementos que tenes de dato.
Tenes que hacer una funcion como la que dije arriba que "transforme" esos 3 vectores en una matriz. Hay otras alternativas, yo nombre la que use yo.
Saludos
|
|
|
|
|
|
|
|
|
Daniela
Nivel 7
Edad: 37
Registrado: 19 Ene 2006
Mensajes: 420
Ubicación: Un reino Muy Muy Lejano
Carrera: Civil y Sistemas
|
|
Cita:
|
Aprovecho la ocasion para preguntar cuando hay que entregarlo el TP.
Muchas gracias.
|
Es la próxima semana, depende del día de tu práctica la fecha de entrega (el martes/miercoles o jueves)
Cita:
|
Te referis a que dijieron que a Ase puede escribir como A= D-L- U????
|
Yo quería hacer eso, pero si te fijás en G-S hay que invertir una matriz y no es muy lindo hacer eso.
Una pregunta, por que arreglos bidimensionales no?
|
|
|
|
_________________
"Pensá en Rosa" by Edgar
|
|
|
|
|
vir123
Nivel 5
Edad: 34
Registrado: 26 Feb 2008
Mensajes: 178
Ubicación: Almagro
Carrera: Informática
|
|
Daniela escribió:
|
Una pregunta, por que arreglos bidimensionales no?
[/color][/b]
|
En el punto a del tp fijate que dice explícitamente "dichas matrices An no deberán estar represetnadas computacionalmente por un arreglo bidimensional sino que las fórmulas de los métodos deberán escribirse para estas matrices en particular"
|
|
|
|
_________________ Get busy living, or get busy dying
|
|
|
|
|
vir123
Nivel 5
Edad: 34
Registrado: 26 Feb 2008
Mensajes: 178
Ubicación: Almagro
Carrera: Informática
|
|
Tengo un problema: uso el metodo sor, y graficamente encuentro un w optimo de aprox 1,6 para un par de matrices y de 1,3 para las otras dos, pero con la formulita me da 1,2 o 1,1 para todas. Alguien tiene idea de como puedo justificar eso?
|
|
|
|
_________________ Get busy living, or get busy dying
|
|
|
|
|
ale_vans
Nivel 7
Edad: 35
Registrado: 22 May 2008
Mensajes: 304
Ubicación: Vte. Lopez
Carrera: No especificada
|
|
creo q la manera de justificar es diciendo que las formulas para obtener el radio espectral de GS es solo una aproximacion, y que el calculo con la formulita del w optimo tmb, asi q es logico q no te de exactamente lo q deberia dar.
Es mas, no siempre se verifica que ro(J)^2 = ro(gs).
|
|
|
|
|
|
|
|
|
vir123
Nivel 5
Edad: 34
Registrado: 26 Feb 2008
Mensajes: 178
Ubicación: Almagro
Carrera: Informática
|
|
ale_vans escribió:
|
creo q la manera de justificar es diciendo que las formulas para obtener el radio espectral de GS es solo una aproximacion, y que el calculo con la formulita del w optimo tmb, asi q es logico q no te de exactamente lo q deberia dar.
Es mas, no siempre se verifica que ro(J)^2 = ro(gs).
|
en realidad me había equivocado en algo y ahora me dio bastante bien, igual gracias
|
|
|
|
_________________ Get busy living, or get busy dying
|
|
|
|
|
|