Grupo de Bosónicos
Códigos de aprendizaje
- Script que calcula el factor de convergencia de Q, para factor entre resoluciones = 2, es decir, el
Q de 2^Q. El formato del archivo de datos es el siguiente: i) seis columnas, 1,3,5 con la coordenada
espacial, 2,4,6 con la res baja, media y alta respectivamente. Cada bloque de datos esta separado por
tres espacios, o pueden ser dos espacios y un hash (#) desde el primer bloque; ii) Los datos deben
tener el downsampling de la resolución mas baja, es decir, deben estar salvados ya
solamente en las posiciones donde se pueden comparar. Yo llamé phi.xl a mi archivo de datos y el
script se puede modificar para que le llames como quieras. Aquí el script y un archivo de parámetros. Sobra decir que lo programé en fortran 90.
Compilar y ejecutar "ejecutable < parametros". También sobra decir que espero que los estudiantes
generalicen este programa para otros casos y todo eso.
- Resuelve la ecuación de onda con una discretización
simple
simple.tgz. Este código es parte del trabajo enviado a la
Rev. Mex. Fis. E ``Solución de la ecuación de onda como un
problema de valores iniciales usando diferencias finitas''.
- Resuelve la ecuación de onda usando variables de primer orden
1st.tgz. Este código es parte del trabajo enviado a la Rev.
Mex. Fis. E ``Solución de la ecuación de onda como un
problema de valores iniciales usando diferencias finitas''.
- Resuelve las ecuaciones de Einstein para un agujero de gusano
soportado por un campo escalar ghost y perturbado. Se trata de uno
de los dos codigos independientes utilizados para mostrar que estos agujeros de gusano
son inestables, resultado que fue publicado en:
``Instability of wormholes supported by a ghost scalar field.
II. Nonlinear evolution.''
J. A. Gonzalez, F. S. Guzman, O. Sarbach.
Class. Quantum Grav. {\bf 26}
(2009) 015011 (20pp)
El codigo esta escrito en Fortran 90 wh_cd.tgz
,
y se recomienda leer
el
README. Funciona de manera optima en plataformas UNIX.
- Resuelve la ecuación de onda en 1+1 usando un método
implícito
Implicit_wave1p1.f90.
Un archivo de parámetros es
Implicit_wave1p1_input.par.
- Resuelve la ecuación de onda en 1+1 usando un método simple (la
discretizacion mas sencilla, el primer código)
Simple_wave1p1.tgz.
- Resuelve la ecuación de onda en 1+1 usando un método
explícito Runge-Kutta de tercer orden
Explicit_wave1p1.f90.
Un archivo de parámetros es
Explicit_wave1p1_input.par.
Un manuscrito que contiene información muy general respecto a
este caso es Guzman.pdf.
- Resuelve la ecuación de onda en 1+1 usando un método
explícito Runge-Kutta de tercer orden
pero tiene una estructura modular, es decir, contiene un Makefile
que hay que modificar modular.tgz.
Un manuscrito que contiene información muy general respecto a
este caso es Guzman.pdf,
y de hecho este código fue usado para lograr la primera parte de
ese trabajo.
Aquí hay un README.
- Resuelve la ecuación de onda en 3+1 con simetría
esférica usando un método explícito Runge-Kutta de
tercer orden
Explicit_wave3p1.f90.
Un archivo de parámetros es
Explicit_wave3p1_input.par.
- Resuelve la ecuación de Schroedinger en 1+1
usando el metodo implícito. El sistema es el de una partícula en
una caja 1-dimensional. El código fuente es Implicit_Schroedinger.f90. Un
archivo de parámetros es implicit_schroedinger_input.par
.
- Resuelve la ecuación de Schroedinger en 1+1
usando el metodo explícito RK3. El sistema es el de una partícula
en una caja 1-dimensional. El código fuente es Explicit_Schroedinger.f90. Un
archivo de parámetros es explicit_schroedinger_input.par
.
- Resuelve la ecuación de Schroedinger en 2+1
usando el metodo explícito RK3. El sistema es el de una partícula
en una caja 2-dimensional. El código fuente es Explicit_Schroedinger2D.f90. Un
archivo de parámetros es explicit_schroedinger_input2d.par
.
- Resuelve la ecuación de Schroedinger en 1+1
para el oscilador armónico usando el metodo explícito RK3. El
código fuente es OsciladorArmonico1D.f90. Un
archivo de parámetros es OsciladorArmonico1D_input.par
.
- El código ShockTube.f90
presenta la solución del problema de Shock Tube de Sod. El
archivo de parámetros es ShockTube_input.par. Para
entender el problema ve a esta
página y para
comparar los resultados obtenidos de ejecutar este archivo de
parámetros ve a esta
página. A este código hay que ponerle condiciones de
frontera reflejantes, para hacer la prueba final. Este código usa
disipación, o sea que no es lo máximo, hay que programar
un HRSC.