def simulate(self):
starttime = time.time()
timings = np.zeros(self.nspinup)
bgcTime = np.zeros(self.ntimestep)
#bgcStepTime = np.zeros(self.ntimestep)
multTime = np.zeros(self.ntimestep)
interpolationTime = np.zeros(self.ntimestep)
for spin in range(self.nspinup):
for step in range(self.ntimestep):
if((step == 2879) & (step > 0)):
N = self.U_PODN.dot(self.y_redN)
DOP = self.U_PODDOP.dot(self.y_redDOP)
#y_hig = io.read_PETSc_vec(self.monitor_path % (spin,step))
io.write_PETSc_vec(N, self.out_pathN % (spin,step))
io.write_PETSc_vec(DOP, self.out_pathDOP % (spin,step))
print("time: ", time.time() - starttime ,"spin: ", spin,"step: ", step,"t:", self.t, "spinup norm: ", np.linalg.norm(N-self.yN) ,np.linalg.norm(DOP-self.yDOP))
timings[spin] = time.time() - starttime
starttime = time.time()
self.yN = N
self.yDOP = DOP
bgcTime[step],multTime[step],interpolationTime[step] = self.__TimeStep(step)
#self.__TimeStep(step)
np.save(self.out_pathN %(spin,step) + "_bgc_timeings.npy",bgcTime)
np.save(self.out_pathN %(spin,step) + "_timeings.npy",timings)
np.save(self.out_pathN %(spin,step) + "_mult_timeings.npy",multTime)
np.save(self.out_pathN %(spin,step) + "_interpolation_timeings.npy",interpolationTime)
np.save(self.out_pathN %(spin,step) + "_timeings.npy",timings)
def simulate(self):
starttime = time.time()
timings = np.zeros(self.nspinup)
bgcTime = np.zeros(self.ntimestep)
#bgcStepTime = np.zeros(self.ntimestep)
multTime = np.zeros(self.ntimestep)
interpolationTime = np.zeros(self.ntimestep)
for spin in range(self.nspinup):
for step in range(self.ntimestep):
if ((step % 240 == 239) & (step > 0)):
matrixIndex = np.int(np.floor(step / 240))
y = np.dot(self.U_POD[(matrixIndex - 1) % 12, :],
self.y_red)
y_hig = io.read_PETSc_vec(self.monitor_path % (spin, step))
io.write_PETSc_vec(y, self.out_path % (spin, step))
print("time: ",
time.time() - starttime, "spin: ", spin, "step: ",
step, "t:", self.t, "error norm: ",
np.linalg.norm(y - y_hig), "spinup norm: ",
np.linalg.norm(y - self.y))
timings[spin] = time.time() - starttime
starttime = time.time()
self.y = y
#bgcTime[step],multTime[step],interpolationTime[step] = self.__TimeStep(step)
self.__TimeStep(step)
#np.save(self.out_path %(spin,step) + "_bgc_timeings.npy",bgcTime)
#np.save(self.out_path %(spin,step) + "_bgcStep_timeings.npy",bgcStepTime)
#np.save(self.out_path %(spin,step) + "_mult_timeings.npy",multTime)
#np.save(self.out_path %(spin,step) + "_interpolation_timeings.npy",interpolationTime)
np.save(self.out_path % (spin, step) + "_timeings.npy", timings)
def simulate(self):
starttime = time.time()
timings = np.zeros(self.nspinup)
bgcTime = np.zeros(self.ntimestep)
#bgcStepTime = np.zeros(self.ntimestep)
multTime = np.zeros(self.ntimestep)
interpolationTime = np.zeros(self.ntimestep)
for spin in range(self.nspinup):
for step in range(self.ntimestep):
if((step == 2879) & (step > 0)):
y = self.U_POD.dot(self.y_red)
y_hig = io.read_PETSc_vec(self.monitor_path % (spin,step))
io.write_PETSc_vec(y, self.out_path % (spin,step))
print("time: ", time.time() - starttime ,"spin: ", spin,"step: ", step,"t:", self.t, "spinup norm: ", np.linalg.norm(y-self.y), "norm: ", np.linalg.norm(y_hig -y))
timings[spin] = time.time() - starttime
starttime = time.time()
self.y = y
bgcTime[step],multTime[step],interpolationTime[step] = self.__TimeStep(step)
#self.__TimeStep(step)
np.save(self.out_path %(spin,step) + "_bgc_timeings.npy",bgcTime)
np.save(self.out_path %(spin,step) + "_timeings.npy",timings)
np.save(self.out_path %(spin,step) + "_mult_timeings.npy",multTime)
np.save(self.out_path %(spin,step) + "_interpolation_timeings.npy",interpolationTime)
np.save(self.out_path %(spin,step) + "_timeings.npy",timings)
def simulate(self):
starttime = time.time()
timings = np.zeros(self.nspinup)
bgcTime = np.zeros(self.ntimestep)
#bgcStepTime = np.zeros(self.ntimestep)
multTime = np.zeros(self.ntimestep)
interpolationTime = np.zeros(self.ntimestep)
for spin in range(self.nspinup):
for step in range(self.ntimestep):
if((step % 240 == 239) & (step > 0)):
matrixIndex = np.int(np.floor(step/240))
y = np.dot(self.U_POD[(matrixIndex-1) % 12,:],self.y_red)
y_hig = io.read_PETSc_vec(self.monitor_path % (spin,step))
io.write_PETSc_vec(y, self.out_path % (spin,step))
print("time: ", time.time() - starttime ,"spin: ", spin,"step: ", step,"t:", self.t,"error norm: ", np.linalg.norm(y-y_hig), "spinup norm: ", np.linalg.norm(y-self.y))
timings[spin] = time.time() - starttime
starttime = time.time()
self.y = y
#bgcTime[step],multTime[step],interpolationTime[step] = self.__TimeStep(step)
self.__TimeStep(step)
#np.save(self.out_path %(spin,step) + "_bgc_timeings.npy",bgcTime)
#np.save(self.out_path %(spin,step) + "_bgcStep_timeings.npy",bgcStepTime)
#np.save(self.out_path %(spin,step) + "_mult_timeings.npy",multTime)
#np.save(self.out_path %(spin,step) + "_interpolation_timeings.npy",interpolationTime)
np.save(self.out_path %(spin,step) + "_timeings.npy",timings)
def test(self,nspinup,ntimestep):
y = np.ones(52749,dtype=np.float_) * 2.17
#load high dim matrices
Ae = []
Ai = []
print("load matrices")
for i in range(12):
print(i)
Ai.append(io.read_PETSc_mat('data/TMM/2.8/Transport/Matrix5_4/1dt/Ai_'+str(i).zfill(2)+'.petsc'))
Ae.append(io.read_PETSc_mat('data/TMM/2.8/Transport/Matrix5_4/1dt/Ae_'+str(i).zfill(2)+'.petsc'))
#check if q is zero in fortran routine
q = np.zeros(52749,dtype=np.float_)
t = 0
#q_select = np.zeros(p.shape[0],dtype=np.float_)
starttime =time.time()
for spin in range(nspinup):
for step in range(ntimestep):
t = np.fmod(0 + step*self.dt, 1.0);
counter = 0
for i in range(4448):
self.bc[0] = self.lat[i]
self.bc[1] = self.interpolation_a[step]*self.fice[i,self.interpolation_j[step]] + self.interpolation_b[step]*self.fice[i,self.interpolation_k[step]]
q[self.J[i]:self.J[i+1]] = modeln.metos3dbgc(self.dt,t,y[self.J[i]:self.J[i+1]],self.u,self.bc,self.dc[self.J[i]:self.J[i+1],:])[:,0]
#print("q:", q[self.J[i]:self.J[i+1]])
Aiint = self.interpolation_a[step]*Ai[self.interpolation_j[step]] + self.interpolation_b[step]*Ai[self.interpolation_k[step]]
Aeint = self.interpolation_a[step]*Ae[self.interpolation_j[step]] + self.interpolation_b[step]*Ae[self.interpolation_k[step]]
#Aey = io.read_PETSc_vec("simulation/compare/Aey_sp%.4dts%.4dN.petsc" % (spin,step))
#Aeq = io.read_PETSc_vec("simulation/compare/Ae+q_sp%.4dts%.4dN.petsc" % (spin,step))
#q_v = io.read_PETSc_vec("simulation/compare/q_sp%.4dts%.4dN.petsc" % (spin,step))
#Aiint_metos = io.read_PETSc_mat("simulation/compare/A%.4d.petsc" % (step))
#print("norm A interplaton: ", (Aiint-Aiint_metos))
ye = Aeint.dot(y)
yeq = ye +q
#io.write_PETSc_vec(yeq,"yeqts%.4dN.petsc" % step)
# A_saved = io.read_PETSc_mat("Ai_interpolatedts%.4d.petsc" % step)
y_j = Aiint.dot(yeq)
#print("q:", np.linalg.norm(q_v-q))
#print("before Ai:", np.linalg.norm(Aeq-yeq))
#print(step,spin)
if(step == 2879):
v = io.read_PETSc_vec("simulation/POD_DEIM/sp%.4dts%.4dN.petsc" % (spin,step))
print("time: ", time.time() - starttime ,spin,step,np.linalg.norm(y-v))
io.write_PETSc_vec(y_j,"simulation/compare/exp01/sp%.4dts%.4dN.petsc" % (spin,step))
starttime = time.time()
#io.write_PETSc_mat(Aiint,"Ai%.4dN.petsc" % step)
# print(Aiint-A_saved.T)
y = y_j
def simulate(self):
starttime = time.time()
timings = np.zeros(self.nspinup)
bgcTime = np.zeros(self.ntimestep)
#bgcStepTime = np.zeros(self.ntimestep)
multTime = np.zeros(self.ntimestep)
interpolationTime = np.zeros(self.ntimestep)
for spin in range(self.nspinup):
for step in range(self.ntimestep):
if ((step == 2879) & (step > 0)):
N = self.U_PODN.dot(self.y_redN)
DOP = self.U_PODDOP.dot(self.y_redDOP)
#y_hig = io.read_PETSc_vec(self.monitor_path % (spin,step))
io.write_PETSc_vec(N, self.out_pathN % (spin, step))
io.write_PETSc_vec(DOP, self.out_pathDOP % (spin, step))
print("time: ",
time.time() - starttime, "spin: ", spin, "step: ",
step, "t:", self.t, "spinup norm: ",
np.linalg.norm(N - self.yN),
np.linalg.norm(DOP - self.yDOP))
timings[spin] = time.time() - starttime
starttime = time.time()
self.yN = N
self.yDOP = DOP
bgcTime[step], multTime[step], interpolationTime[
step] = self.__TimeStep(step)
#self.__TimeStep(step)
np.save(self.out_pathN % (spin, step) + "_bgc_timeings.npy",
bgcTime)
np.save(self.out_pathN % (spin, step) + "_timeings.npy", timings)
np.save(self.out_pathN % (spin, step) + "_mult_timeings.npy",
multTime)
np.save(
self.out_pathN % (spin, step) + "_interpolation_timeings.npy",
interpolationTime)
np.save(self.out_pathN % (spin, step) + "_timeings.npy", timings)