That is the wrong mental model for using NumPy efficiently. NumPy arrays are stored in contiguous blocks of memory. To append rows or columns to an existing array, the entire array needs to be copied to a new block of memory, creating gaps for the new elements to be stored. This is very inefficient if done repeatedly.
Instead of appending rows, allocate a suitably sized array, and then assign to it row-by-row:
>>> import numpy as np
>>> a = np.zeros(shape=(3, 2))
>>> a
array([[ 0., 0.],
[ 0., 0.],
[ 0., 0.]])
>>> a[0] = [1, 2]
>>> a[1] = [3, 4]
>>> a[2] = [5, 6]
>>> a
array([[ 1., 2.],
[ 3., 4.],
[ 5., 6.]])
That is the wrong mental model for using NumPy efficiently. NumPy arrays are stored in contiguous blocks of memory. To append rows or columns to an existing array, the entire array needs to be copied to a new block of memory, creating gaps for the new elements to be stored. This is very inefficient if done repeatedly.
Instead of appending rows, allocate a suitably sized array, and then assign to it row-by-row:
>>> import numpy as np
>>> a = np.zeros(shape=(3, 2))
>>> a
array([[ 0., 0.],
[ 0., 0.],
[ 0., 0.]])
>>> a[0] = [1, 2]
>>> a[1] = [3, 4]
>>> a[2] = [5, 6]
>>> a
array([[ 1., 2.],
[ 3., 4.],
[ 5., 6.]])
A NumPy array is a very different data structure from a list and is designed to be used in different ways. Your use of hstack is potentially very inefficient... every time you call it, all the data in the existing array is copied into a new one. (The append function will have the same issue.) If you want to build up your matrix one column at a time, you might be best off to keep it in a list until it is finished, and only then convert it into an array.
e.g.
mylist = []
for item in data:
mylist.append(item)
mat = numpy.array(mylist)
item can be a list, an array or any iterable, as long
as each item has the same number of elements.
In this particular case (data is some iterable holding the matrix columns) you can simply use
mat = numpy.array(data)
(Also note that using list as a variable name is probably not good practice since it masks the built-in type by that name, which can lead to bugs.)
EDIT:
If for some reason you really do want to create an empty array, you can just use numpy.array([]), but this is rarely useful!
How can I create a truly empty numpy array which can be merged onto (by a recursive function)?
In Python in what ways can u make an empty NumPy array? - Stack Overflow
Creating empty nxn square matrix of 0
numpy.empty equivalent in julia?
I'm kind of stuck conceptually on how to make this happen. I have a recursive method that builds a binary tree, and stores the tree as an instance variable. However, the function is not allowed to return anything, so each recursive call should (according to me) modify in-place the tree instance variable. However, I'm not sure how to set up my instance variable such that all said and done it holds a multidimensional array that represents the tree.
Say I set initialize it as a 1x1 array with element zero as a placeholder. Then as I go about recursing through my tree I can merge to it... but at the end I'm left with a spare [0] element that I don't need. In this case, I'd need some kind of final stop condition and function to remove that unnecessary placeholder stump. I don't think this is possible?
Otherwise, say I initialize the instance variable as None. Then when the first series of recursive calls, it would have to reassign the tree variable to change from None to an ndarray object, but all future calls would have to merge to the array. I don't think this is what the function should be asked to do?
Is there a way to make a truly empty array that I can merge onto? (e.g. np.empty doesn't reallly give an empty array, it gives an array with placeholder values so I'm still left with a useless stump at the end).
If by empty array, you mean an array with 0 dimensions, you can simply convert to a NumPy array an empty list:
import numpy as np
np.array([])
# array([], dtype=float64)
or use any of the NumPy's initialization functions with a (0,) or 0 shape, e.g.:
np.zeros((0,))
# array([], dtype=float64)
np.ones((0,))
# array([], dtype=float64)
np.empty((0,))
# array([], dtype=float64)
np.full((0,), 0.0)
# array([], dtype=float64)
etc.
If by empty you mean an array of given size but not initialized (i.e. you only ask the OS for the required memory), you can use
np.empty() specifying the size you need, e.g.:
np.empty((2, 3))
# array([[1.46643506e-316, 0.00000000e+000, 0.00000000e+000],
# [0.00000000e+000, 0.00000000e+000, 0.00000000e+000]])
There are two ways an empty NumPy array can be created: numpy.zeros and numpy.empty.
The syntax for using numpy.zeros and numpy.empty is shown below:
numpy.zeros(shape, dtype=float, order='C')
numpy.empty(shape, dtype=float, order='C')
# Shape -> Shape of the new array, e.g., (2, 3) or 2.
# dtype -> The desired data-type for the array,e.g., numpy.int8. Default is numpy.float64. This parameter is optional.
# order -> Indicates whether multi-dimensional data should be stored in row-major (C-style) or column-major (Fortran-style) order in memory. This parameter is optional.
how to use numpy.zeros to create an empty NumPy array:
import numpy as np
myArr = np.zeros((2,3))
print(myArr)
how to use numpy.empty to create an empty NumPy array:
import numpy as np
myArr = np.empty((2,3))
print(myArr)
numpy.empty, unlike numpy.zeros, does not set the array values to zero and, therefore, may be marginally faster. On the other hand, it requires the user to manually set all the values in the array and should be used with caution.