Sto lavorando per impostare alcune condizioni al contorno per un modello di falda acquifera, e sono in grado di impostare l'intera prima riga su un valore costante, ma non sull'intera prima colonna. Sto usando np.zeros((11, 1001))
per creare una matrice vuota. Qualcuno sa perché sono riuscito a definire la prima riga, ma non la prima colonna? Ho notato la riga in questione qui sotto.Come impostare la prima colonna su un valore costante di una matrice numPy np.zeros vuota?
import numpy as np
x = range(0, 110, 10)
time = range(0, 5005, 5)
xSize = len(x)
timeSize = len(time)
dx = 10
dt = 5
Sy = 0.1
k = 0.002
head = np.zeros((11, 1001))
head[0:][0] = 16 # sets the first row to 16
head[0][0:] = 16 # DOESN'T set the first column to 16
for t in time:
for i in x[1:len(x)-1]:
head[t+1][i] = head[t][i] + ((dt*k)/(2*Sy)) * (((head[t][i-1]**2) - (2*head[t][i]**2) + (head[t][i+1]**2))/(dx**2))
@Nitish Grazie a entrambi! – AF2k15