Fecha y hora actual: Lunes 23 Sep 2019 08:52
Índice del Foro

Foros de programación informática, diseño gráfico y Web

En esta comunidad intentaremos dar soporte de programación a todos los niveles, desde principiantes a profesionales de la informática, desarrollo de programas, programación web y mucho más.

Paralelizacion en python

Responder al Tema

Índice del Foro > Python > Paralelizacion en python

Autor Mensaje
ssy31



Registrado: 13 Sep 2017
Mensajes: 1

Mensaje Publicado: Miércoles 13 Sep 2017 23:50

Título del mensaje: Paralelizacion en python

Responder citando

Hola a todos me gustaría saber si alguno sabría como paralelizar este programa en python sin usar OpenMP, yo la verdad no tengo ni idea de como poder hacerlo, aquí les dejo el codigo en serial y muchas gracias
Código:
  1. # coding=utf-8
  2. import numpy as np
  3. import matplotlib.pyplot as plt
  4.  
  5. #constantes
  6. Vxo = 0
  7. Vyo = 29780.
  8. xo = 1.5E11
  9. yo = 0
  10. M = 1.98855E30
  11. G = 6.674E-11
  12.  
  13. Xl = np.array([1.5E11])
  14. Yl = np.array([0])
  15. Vxl = np.array([0])
  16. Vyl = np.array([29780.])
  17. El = np.array([])
  18. tl = np.array([])
  19. for i in range(0, 1825 + 1): # numero de dias en 5 años
  20.  
  21.  
  22. h = 86400 # segundos por dia
  23.  
  24. rm = (Xl[-1] ** 2. + Yl[-1] ** 2.) ** 0.5 # radio para Vmed
  25.  
  26. dVxm = -1 * G * M * Xl[-1] / (rm ** 3.) # derivada respecto al tiempo de la Vel.x para Vmed
  27. dVym = -1 * G * M * Yl[-1] / (rm ** 3.) # derivada respecto al tiempo de la Vel.y para Vmed
  28.  
  29. Vximed = Vxl[-1] + (h / 2.) * dVxm # hayando vel intermedia
  30. Vyimed = Vyl[-1] + (h / 2.) * dVym
  31.  
  32. xi = Xl[-1] + h * Vximed # aplicando Leapfrog para x
  33. yi = Yl[-1] + h * Vyimed # aplicando Leapfrog para y
  34.  
  35. r = (xi ** 2. + yi ** 2.) ** 0.5 # radio para Vi+1
  36.  
  37. dVxo = -1 * G * M * xi / (r ** 3.) # derivada respecto al tiempo de la Vel.x para Vi+1
  38. dVyo = -1 * G * M * yi / (r ** 3.) # derivada respecto al tiempo de la Vel.y para Vi+1
  39.  
  40. Vxi = Vximed + (h / 2.) * dVxo # Obteniendo Vx
  41. Vyi = Vyimed + (h / 2.) * dVyo # Obteniendo Vy
  42.  
  43. Ei = (0.5) * (Vxi ** 2 + Vyi ** 2) - G * M / ((xi ** 2. + yi ** 2.) ** 0.5) # Calculando la E total (LeapF)
  44.  
  45. # añadiendo valores a los arreglos
  46. Xl = np.append(Xl, xi)
  47. Yl = np.append(Yl, yi)
  48. Vxl = np.append(Vxl, Vxi)
  49. Vyl = np.append(Vyl, Vyi)
  50. El = np.append(El, Ei)
  51. tl = np.append(tl, i)
  52.  
  53. plt.plot(Xl, Yl)
  54. plt.title("Grafico de posicion en Leap Frog")
  55. plt.xlabel("Posicion en x")
  56. plt.ylabel("Posicion en y")
  57. plt.xlim(-2e11, 2e11)
  58. plt.show()
  59.  
  60. plt.plot(tl, El, color="green", label="Leap Frog Energy")
  61. plt.legend(loc='upper right', shadow=True, fontsize='small')
  62. plt.title("Grafico de energias contra tiempo")
  63. plt.xlabel("Tiempo [dias]")
  64. plt.ylabel("Energia")
  65. plt.show()
  66.  

Volver arriba
Ver perfil del usuario Enviar mensaje privado
Pedrolo
Usuario Inquieto


Registrado: 07 May 2017
Mensajes: 50
Ubicación: Posiblemente lejos de ti

Mensaje Publicado: Jueves 21 Sep 2017 01:24

Título del mensaje: Paralelizacion en python

Responder citando

https://docs.python.org/2/library/multiprocessing.html#module-multiprocessing

Volver arriba
Ver perfil del usuario Enviar mensaje privado Visitar sitio web del autor
Responder al Tema
Mostrar mensajes anteriores:   
Ir a:  
Todas las horas están en GMT + 2 Horas

Temas relacionados

Tema Autor Foros Respuestas Publicado
El foro no contiene ningún mensaje nuevo

Python abrir sin darle al start

Sergi99 Python 0 Lunes 20 May 2019 10:54 Ver último mensaje
El foro no contiene ningún mensaje nuevo

necesito ayuda con el algoritmo de dijkstra en ...

mery Python 0 Jueves 13 Dic 2018 18:07 Ver último mensaje
El foro no contiene ningún mensaje nuevo

Tutorial Python Rápido para Principiantes

Errodringer Python 0 Viernes 14 Sep 2018 19:38 Ver último mensaje
El foro no contiene ningún mensaje nuevo

Duda codigo en Python - Red Neuronal con numpy ...

Daniel_wUNTG Python 1 Sábado 21 Oct 2017 22:27 Ver último mensaje
El foro no contiene ningún mensaje nuevo

Problema para extracción de datos utilizando py...

liz Python 1 Viernes 06 Oct 2017 11:05 Ver último mensaje
Panel de Control
No puede crear mensajes, No puede responder temas, No puede editar sus mensajes, No puede borrar sus mensajes, No puede votar en encuestas,