Python plotSingular Beispiele

Beispiel #1

import numpy as np
import util
import petsc_io as io
import matplotlib
import matplotlib.pyplot as plt

matplotlib.rcParams.update({'font.size': 22})
matplotlib.rcParams.update({'font.weight': 'bold'})
matplotlib.rcParams.update({'lines.linewidth': 2})
matplotlib.rcParams.update({'figure.figsize': (8, 6)})
axis_font = {'fontname': 'Bitstream Vera Sans', 'weight': 'bold', 'size': '22'}

s_deim = io.read_PETSc_vec(
    "reduced_basis/full_trajectory/3_snapshots_per_month/1000/s_DEIM.petsc")

s_pod = io.read_PETSc_vec(
    "reduced_basis/full_trajectory/3_snapshots_per_month/1000/s_POD.petsc")

pod_k = util.computePODError(np.square(s_pod), 1e-5)
deim_m = util.computePODError(np.square(s_deim), 1e-5)
print(pod_k, deim_m)
util.plotSingular(np.square(s_pod), np.square(s_deim), 1e-6, bounds=False)
plt.xlabel(r'\textbf{Sigma }  $i$', **axis_font)
plt.ylabel(r'\textbf{Value}', **axis_font)
plt.legend(["Singular values FS", "Singular values of q"], loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 36\_3000}',**axis_font)
plt.savefig('singularvalues_36_1000.png', bbox_inches='tight')
plt.show()

Beispiel #2

Datei: plot_singular_values_3_3000.py Projekt: neeljp/pod_deim

import numpy as np
import util
import petsc_io as io
import matplotlib
import matplotlib.pyplot as plt

matplotlib.rcParams.update({'font.size': 22})
matplotlib.rcParams.update({'font.weight': 'bold'})
matplotlib.rcParams.update({'lines.linewidth': 2})
matplotlib.rcParams.update({'figure.figsize': (8, 6)})
axis_font = {'fontname':'Bitstream Vera Sans', 'weight':'bold' ,'size':'22'}



s_deim = io.read_PETSc_vec("reduced_basis/full_trajectory/one_month/s_DEIM.petsc")


s_pod = io.read_PETSc_vec("reduced_basis/full_trajectory/one_month/s_POD.petsc")


pod_k  = util.computePODError(np.square(s_pod),1e-5)
deim_m = util.computePODError(np.square(s_deim),1e-5)
print(pod_k,deim_m)
util.plotSingular(np.square(s_pod),np.square(s_deim),1e-6,bounds=False)
plt.xlabel(r'\textbf{Sigma }  $i$',**axis_font)
plt.ylabel(r'\textbf{Value}',**axis_font)
plt.legend(["Singular values FS","Singular values of q"],loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 3\_3000}',**axis_font)
plt.savefig('singularvalues_3_3000.png', bbox_inches='tight')
plt.show()

Beispiel #3

matplotlib.rcParams.update({'font.weight': 'bold'})
matplotlib.rcParams.update({'lines.linewidth': 2})
matplotlib.rcParams.update({'figure.figsize': (8, 6)})
axis_font = {'fontname': 'Bitstream Vera Sans', 'weight': 'bold', 'size': '22'}

s_deimN = io.read_PETSc_vec("reduced_basis/MITgcm/s_DEIMPO4.petsc")

s_podN = io.read_PETSc_vec("reduced_basis/MITgcm/s_PODPO4.petsc")
s_deimDOP = io.read_PETSc_vec("reduced_basis/MITgcm/s_DEIMDOP.petsc")

s_podDOP = io.read_PETSc_vec("reduced_basis/MITgcm/s_PODDOP.petsc")

#pod_k  = util.computePODError(np.square(s_pod),1e-5)
#deim_m = util.computePODError(np.square(s_deim),1e-5)
#print(pod_k,deim_m)
util.plotSingular(np.square(s_podN), np.square(s_deimN), 1e-6, bounds=False)
plt.xlabel(r'\textbf{Sigma }  $i$', **axis_font)
plt.ylabel(r'\textbf{Value}', **axis_font)
plt.legend(["Singular values FS", "Singular values of q"], loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 36\_1000 of tracer N}',**axis_font)
plt.savefig('singularvalues_MITgcm_PO4.png', bbox_inches='tight')
plt.show()

util.plotSingular(np.square(s_podDOP),
                  np.square(s_deimDOP),
                  1e-6,
                  bounds=False)
plt.xlabel(r'\textbf{Sigma }  $i$', **axis_font)
plt.ylabel(r'\textbf{Value}', **axis_font)
plt.legend(["Singular values FS", "Singular values of q"], loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 36\_1000 of tracer DOP}',**axis_font)

Beispiel #4

Datei: plot_singular_values_one_month.py Projekt: neeljp/pod_deim

import numpy as np
import util
import petsc_io as io
import matplotlib
import matplotlib.pyplot as plt

matplotlib.rcParams.update({'font.size': 22})
matplotlib.rcParams.update({'font.weight': 'bold'})
matplotlib.rcParams.update({'lines.linewidth': 2})
matplotlib.rcParams.update({'figure.figsize': (12, 10)})
axis_font = {'fontname':'Bitstream Vera Sans', 'weight':'bold' ,'size':'22'}



s_deim = io.read_PETSc_vec("reduced_basis/full_trajectory/one_month/s_DEIM.petsc")


s_pod = io.read_PETSc_vec("reduced_basis/full_trajectory/one_month/s_POD.petsc")


pod_k  = util.computePODError(np.square(s_pod),1e-5)
deim_m = util.computePODError(np.square(s_deim),1e-5)
print(pod_k,deim_m)
util.plotSingular(np.square(s_pod),np.square(s_deim),1e-6)
plt.xlabel(r'\textbf{Sigma }  $i$',**axis_font)
plt.ylabel(r'\textbf{Value}',**axis_font)
plt.legend(["Singular values FS","Singular values of q"],loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 3\_3000}',**axis_font)
plt.savefig('singularvalues_one_month.png', bbox_inches='tight')
plt.show()

Beispiel #5

Datei: plot_singular_values_MITgcm_3000.py Projekt: neeljp/pod_deim

s_deimN = io.read_PETSc_vec("reduced_basis/MITgcm/s_DEIMPO4.petsc")


s_podN = io.read_PETSc_vec("reduced_basis/MITgcm/s_PODPO4.petsc")
s_deimDOP = io.read_PETSc_vec("reduced_basis/MITgcm/s_DEIMDOP.petsc")


s_podDOP = io.read_PETSc_vec("reduced_basis/MITgcm/s_PODDOP.petsc")


#pod_k  = util.computePODError(np.square(s_pod),1e-5)
#deim_m = util.computePODError(np.square(s_deim),1e-5)
#print(pod_k,deim_m)
util.plotSingular(np.square(s_podN),np.square(s_deimN),1e-6, bounds=False)
plt.xlabel(r'\textbf{Sigma }  $i$',**axis_font)
plt.ylabel(r'\textbf{Value}',**axis_font)
plt.legend(["Singular values FS","Singular values of q"],loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 36\_1000 of tracer N}',**axis_font)
plt.savefig('singularvalues_MITgcm_PO4.png', bbox_inches='tight')
plt.show()

util.plotSingular(np.square(s_podDOP),np.square(s_deimDOP),1e-6, bounds=False)
plt.xlabel(r'\textbf{Sigma }  $i$',**axis_font)
plt.ylabel(r'\textbf{Value}',**axis_font)
plt.legend(["Singular values FS","Singular values of q"],loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 36\_1000 of tracer DOP}',**axis_font)
plt.savefig('singularvalues_MITgcm_DOP.png', bbox_inches='tight')
plt.show()

Beispiel #6

import numpy as np
import util
import petsc_io as io
import matplotlib
import matplotlib.pyplot as plt

matplotlib.rcParams.update({'font.size': 22})
matplotlib.rcParams.update({'font.weight': 'bold'})
matplotlib.rcParams.update({'lines.linewidth': 2})
matplotlib.rcParams.update({'figure.figsize': (12, 10)})
axis_font = {'fontname': 'Bitstream Vera Sans', 'weight': 'bold', 'size': '22'}

s_deim = io.read_PETSc_vec(
    "reduced_basis/full_trajectory/one_month/s_DEIM.petsc")

s_pod = io.read_PETSc_vec(
    "reduced_basis/full_trajectory/one_month/s_POD.petsc")

pod_k = util.computePODError(np.square(s_pod), 1e-5)
deim_m = util.computePODError(np.square(s_deim), 1e-5)
print(pod_k, deim_m)
util.plotSingular(np.square(s_pod), np.square(s_deim), 1e-6)
plt.xlabel(r'\textbf{Sigma }  $i$', **axis_font)
plt.ylabel(r'\textbf{Value}', **axis_font)
plt.legend(["Singular values FS", "Singular values of q"], loc='upper right')
#plt.title(r'\textbf{Singularvalues of SVD 3\_3000}',**axis_font)
plt.savefig('singularvalues_one_month.png', bbox_inches='tight')
plt.show()