Este lunes a partir de las 11:00 hs estaremos entregando los 1eros recuperatorios y haciendo consultas para el 2do, en el aula de siempre o en el bar si estuviera ocupada.

Les recordamos que para firmar libretas ademas es obligatorio tener completa la encuesta docente.

Este viernes (01/12), yo (Quimey) voy a estar disponible para consultas de 16:30 a 18:30. Pueden pasar

a buscarme por el Laboratorio de Procesado de Imagenes (LPI). En caso de no encontrarme ahí, fijense en el bar.

Este viernes 24/11 daremos consultas extras para el parcial, en distintos horarios.

• Quimey estará disponible en el LPI  desde las 11:00 hs hasta las 12:30 hs.
• Pablo estará en la bar del I de 15:00 a 17:00 hs.
• Martin estará en el bar del I de 16:30 a 18:30 hs.

Los esperamos.

con este generan los datos:ising_quf

Aqui les dejo una corrida de 32×32 sitios, con temperaturas de 1 a 5 con paso 0.1 y npre=2000 npasos=8000:ising_dat.tar

con este grafican: print_ising

IMPORTANTE… a veces surge un error al generar los vectores "uposta" y "cvposta" (solucion de Onsager). Si ese es el caso, simplemente usen los archivos homonimos en el zip para "loadear". (esta accion esta comentada en el codigo ising_print, solo hace falta descomentar, con ctrl+5)

La solucion de Onsager esta en terminos de integrales elipticas… no logré una manera simple de implementar estas, pero si un paper (https://arxiv.org/abs/1411.2495) donde explican como ponerlas en funcion de hipergeométricas, las cuales si pude implementar.

Finalmente aqui dejo un programita para visualizar la convergencia …convergencia_cumsum

Aqui la energia:

aqui la diferencia entre onsager y metropolis:

vemos que es maxima cerca de Tc. Esto se debe es donde mas le cuesta converger al algoritmo.

Aqui el cv:

y la magnetizacion:

Finalemente vemos como difieren para un sistema bien chico (4×4) donde claramente metropolis difiere muchisimo del cv analitico (correspondiente al limite termodinámico)

Les dejamos, a modo de ejemplo, el 2do parcial que tomamos el año pasado. Resuelvanlo y traigan sus consultas para la próxima clase.  Aca el parcial 2doParcial .

Les dejo los programas de Ising en Matlab (ising2D ising2Dpaso En ) y en Python ( ising2D ). 

Les dejamos un programa escrito en Matlab/Octave p2Bose para calcular las propiedades del gas de Bose (\(N_0, U,C_v\),etc) del problema 2 de la guia en función de número de partículas y temperatura. 

A modo de ejemplo les dejamos unos gráficos producidos por el programa p2Bose para el gas.

Continue reading →

Este lunes tendremos primero la práctica, y luego la teórica. Empezamos la práctica a las 10:00 y la teórica será a las 13:00. 

Nos vemos.

Este viernes de 29/09 de 16:00 a 19:00 hs daremos consultas para el parcial de este lunes. Los esperamos en el bar del 1.

Este viernes 22/09 daremos consultas extras para el 1er parcial en distintos horarios:

• Quimey estará disponible en el LPI (Laboratorio de Procesamiento de Imágenes) de 11:00 a 13:00hs.
• Martin estará en el Bar de Pab. I de 16:00 a 18:00 hs.
• Pablo estará en el Bar del Pab. I de 13:30 a 15:30hs.

Los esperamos.