Skip to main content
se eliminaron 178 caracteres en el cuerpo
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez. Se puede hacer uso de numpy.ndarray.resize() que hace uso de realloc (C), pero nada nos asegura que aún así no sea necesario reasignar y copiar todo el bloque de memoria si no es posible aumentar su tamaño.

Dado que sabes de antemano el tamaño de X e Y es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. Puedes hacer esto con numpy.empty es una buena opción en este caso, simplemente reserva la memoria necesaria para el(no inicializa los items del array sin inicializar ningún valor), dado que vamos a iterar sobre él al completo inicializar (np.zeros /o np.ones) no es necesario y es menos eficiente:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.zeros(50000)
y = np.zeros(50000)

for n in range(50000):
    r = random.uniform(0, 100)

    if r < 1.0:
        x[n] = 0
        y[n] = 0.16 * y[n - 1]

    elif r < 86.0:
        x[n] = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y[n] = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6

    elif r < 93.0:
        x[n] = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y[n] = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6

    else:
        x[n] = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y[n] = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44


'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez. Se puede hacer uso de numpy.ndarray.resize() que hace uso de realloc (C), pero nada nos asegura que aún así no sea necesario reasignar y copiar todo el bloque de memoria si no es posible aumentar su tamaño.

Dado que sabes de antemano el tamaño de X e Y es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. numpy.empty es una buena opción en este caso, simplemente reserva la memoria necesaria para el array sin inicializar ningún valor, dado que vamos a iterar sobre él al completo inicializar (np.zeros / np.ones) no es necesario y es menos eficiente:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.zeros(50000)
y = np.zeros(50000)

for n in range(50000):
    r = random.uniform(0, 100)

    if r < 1.0:
        x[n] = 0
        y[n] = 0.16 * y[n - 1]

    elif r < 86.0:
        x[n] = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y[n] = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6

    elif r < 93.0:
        x[n] = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y[n] = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6

    else:
        x[n] = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y[n] = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44


'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez. Se puede hacer uso de numpy.ndarray.resize() que hace uso de realloc (C), pero nada nos asegura que aún así no sea necesario reasignar y copiar todo el bloque de memoria si no es posible aumentar su tamaño.

Dado que sabes de antemano el tamaño de X e Y es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. Puedes hacer esto con numpy.empty (no inicializa los items del array), np.zeros o np.ones:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.zeros(50000)
y = np.zeros(50000)

for n in range(50000):
    r = random.uniform(0, 100)

    if r < 1.0:
        y[n] = 0.16 * y[n - 1]

    elif r < 86.0:
        x[n] = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y[n] = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6

    elif r < 93.0:
        x[n] = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y[n] = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6

    else:
        x[n] = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y[n] = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44


'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()
se añadieron 335 caracteres en el cuerpo
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez. Se puede hacer uso de numpy.ndarray.resize() que hace uso de realloc (C), pero nada nos asegura que aún así no sea necesario reasignar y copiar todo el bloque de memoria si no es posible aumentar su tamaño.

Dado que sabes de antemano el tamaño de XX e YY es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. numpy.empty es una buena opción en este caso, simplemente reserva la memoria necesaria para el array sin inicializar ningún valor, dado que vamos a iterar sobre él al completo inicializar (np.zeros / np.ones) no es necesario y es menos eficiente:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.emptyzeros(50000)
y = np.emptyzeros(50000)

for n in range(50000):
    r = random.uniform(0, 100) 

    if r < 1.0:
        x_x[n] = 0
        y_y[n] = 0.16 * y[n - 1]

    elif r < 86.0:
        x_x[n] = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y_y[n] = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6 

    elif r < 93.0:
        x_x[n] = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y_y[n] = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6 

    else:
        x_x[n] = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y_y[n] = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44
 
    x[n] = x_
    y[n] = y_

'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez.

Dado que sabes de antemano el tamaño de X e Y es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. numpy.empty es una buena opción en este caso, simplemente reserva la memoria necesaria para el array sin inicializar ningún valor:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.empty(50000)
y = np.empty(50000)

for n range(50000):
    r = random.uniform(0, 100)
    if r < 1.0:
        x_ = 0
        y_ = 0.16 * y[n - 1]

    elif r < 86.0:
        x_ = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y_ = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6
    elif r < 93.0:
        x_ = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y_ = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6
    else:
        x_ = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y_ = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44
 
    x[n] = x_
    y[n] = y_

'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez. Se puede hacer uso de numpy.ndarray.resize() que hace uso de realloc (C), pero nada nos asegura que aún así no sea necesario reasignar y copiar todo el bloque de memoria si no es posible aumentar su tamaño.

Dado que sabes de antemano el tamaño de X e Y es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. numpy.empty es una buena opción en este caso, simplemente reserva la memoria necesaria para el array sin inicializar ningún valor, dado que vamos a iterar sobre él al completo inicializar (np.zeros / np.ones) no es necesario y es menos eficiente:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.zeros(50000)
y = np.zeros(50000)

for n in range(50000):
    r = random.uniform(0, 100) 

    if r < 1.0:
        x[n] = 0
        y[n] = 0.16 * y[n - 1]

    elif r < 86.0:
        x[n] = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y[n] = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6 

    elif r < 93.0:
        x[n] = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y[n] = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6 

    else:
        x[n] = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y[n] = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44


'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

NumPy tiene su propio método append, no obstante modificar el tamaño de un array NumPy es muy ineficiente ya que no se puede hacer in-place (se crea un nuevo objeto) copiando un bloque de memoria entero cada vez.

Dado que sabes de antemano el tamaño de X e Y es mucho más simple y eficiente inicializar los arrays con ese tamaño y luego asignar los valores mediante indizado. numpy.empty es una buena opción en este caso, simplemente reserva la memoria necesaria para el array sin inicializar ningún valor:

import random
import numpy as np
import matplotlib.pyplot as plt



x = np.empty(50000)
y = np.empty(50000)

for n range(50000):
    r = random.uniform(0, 100)
    if r < 1.0:
        x_ = 0
        y_ = 0.16 * y[n - 1]

    elif r < 86.0:
        x_ = 0.85 * x[n - 1] + 0.04 * y[n - 1]
        y_ = -0.04 * x[n - 1] + 0.85 * y[n - 1] + 1.6
    elif r < 93.0:
        x_ = 0.2 * x[n - 1] - 0.26 * y[n - 1]
        y_ = 0.23 * x[n - 1] + 0.22 * y[n - 1] + 1.6
    else:
        x_ = -0.1 * x[n - 1] + 0.28 * y[n - 1]
        y_ = 0.26 * x[n - 1] + 0.24 * y[n - 1] + 0.44

    x[n] = x_
    y[n] = y_

'''Make a plot'''
plt.figure(figsize=[50, 50])
plt.scatter(x, y, color='g', marker='.')
plt.show()