ROM object of type DMD (see libROM) More...

#include <rom_object_dmd.h>

Inheritance diagram for DMDROMObject:

Public Member Functions
	DMDROMObject (const int, const double, const int, const int, const int, const int, const int a_var_idx=-1)

virtual	~DMDROMObject ()

virtual void	projectInitialSolution (CAROM::Vector &a_U)

virtual void	takeSample (const CAROM::Vector &, const double)

virtual void	train ()

virtual const CAROM::Vector *const	predict (const double a_t) const

virtual void	save (const std::string &a_fname_root) const

virtual void	load (const std::string &a_fname_root)

Public Member Functions inherited from ROMObject
virtual	~ROMObject ()

Protected Attributes
std::vector< CAROM::DMD * >	m_dmd

std::vector< bool >	m_dmd_is_trained

std::vector< Interval >	m_intervals

int	m_rank

int	m_nproc

int	m_vec_size

double	m_dt

double	m_t_final

int	m_rdim

int	m_num_window_samples

int	m_tic

int	m_curr_win

int	m_sim_idx

int	m_var_idx

bool	m_write_snapshot_mat

std::string	m_dirname

Detailed Description

ROM object of type DMD (see libROM)

ROM object of type DMD (see libROM) DMD-type ROM (see libROM for details)

DMD-type ROM (see libROM for details). DMD stands for "Dynamic Mode Decomposition".

Definition at line 37 of file rom_object_dmd.h.

Constructor & Destructor Documentation

◆ DMDROMObject()

DMDROMObject	(	const int	a_vec_size,
		const double	a_dt,
		const int	a_rdim,
		const int	a_rank,
		const int	a_nproc,
		const int	a_sim_idx,
		const int	a_var_idx = `-1`
	)

Constructor

Constructor This function will also look into the file librom.inp for DMD-specific options. Rank 0 reads in the inputs and broadcasts them to all the processors.

The format of librom.inp is as follows:

begin
    <keyword>   <value>
    ...
    <keyword>   <value>
end

where the list of keywords and their type (that this function will look for) are:

Keyword name	Type	Variable	Default value --—
dmd_num_win_samples	int	DMDROMObject::m_num_window_samples	INT_MAX
dmd_dirname	string	DMDROMObject::m_dirname	"DMD"
dmd_write_snapshot_mat	bool	DMDROMObject::m_write_snapshot_mat	false

Note: other keywords in this file may be read by other functions.

Parameters

a_vec_size	vector size
a_dt	time step size
a_rdim	latent space dimension
a_rank	MPI rank of this process
a_nproc	Number of MPI processes
a_sim_idx	Simulation index (for ensemble simulations
a_var_idx	Vector component index

Definition at line 35 of file DMDROMObject.cpp.

 {
   m_rank = a_rank;
   m_nproc = a_nproc;
 
   m_dmd.clear();
   m_dmd_is_trained.clear();
   m_intervals.clear();
   m_vec_size = a_vec_size;
   m_dt = a_dt;
   m_t_final = -1;
   m_rdim = a_rdim;
 
   m_sim_idx = a_sim_idx;
   m_var_idx = a_var_idx;
 
   m_num_window_samples = INT_MAX;
 
   char dirname_c_str[_MAX_STRING_SIZE_] = "DMD";
   char write_snapshot_mat[_MAX_STRING_SIZE_] = "false";
 
   if (!m_rank) {
 
     FILE *in;
     in = fopen(_LIBROM_INP_FNAME_,"r");
 
     if (in) {
 
       int ferr;
       char word[_MAX_STRING_SIZE_];
       ferr = fscanf(in,"%s", word); if (ferr != 1) return;
 
       if (std::string(word) == "begin") {
         while (std::string(word) != "end") {
           ferr = fscanf(in,"%s",word); if (ferr != 1) return;
           if (std::string(word) == "dmd_num_win_samples") {
             ferr = fscanf(in,"%d", &m_num_window_samples); if (ferr != 1) return;
           } else if (std::string(word) == "dmd_dirname") {
             ferr = fscanf(in,"%s", dirname_c_str); if (ferr != 1) return;
           } else if (std::string(word) == "dmd_write_snapshot_mat") {
             ferr = fscanf(in,"%s", write_snapshot_mat); if (ferr != 1) return;
           }
           if (ferr != 1) return;
         }
       } else {
         fprintf( stderr, "Error: Illegal format in file \"%s\". Word read is: %s\n",
                  _LIBROM_INP_FNAME_, word);
         return;
       }
 
       fclose(in);
 
     }
 
     /* print useful stuff to screen */
     printf("DMDROMObject details:\n");
     printf("  number of samples per window:   %d\n", m_num_window_samples);
     printf("  directory name for DMD onjects: %s\n", dirname_c_str);
     printf("  write snapshot matrix to file:  %s\n", write_snapshot_mat);
     if (m_sim_idx >= 0) {
       printf("  simulation domain:  %d\n", m_sim_idx);
     }
     if (m_var_idx >= 0) {
       printf("  Vector component:  %d\n", m_var_idx);
     }
   }
 
 #ifndef serial
   MPI_Bcast(&m_num_window_samples,1,MPI_INT,0,MPI_COMM_WORLD);
   MPI_Bcast(dirname_c_str,_MAX_STRING_SIZE_,MPI_CHAR,0,MPI_COMM_WORLD);
   MPI_Bcast(write_snapshot_mat,_MAX_STRING_SIZE_,MPI_CHAR,0,MPI_COMM_WORLD);
 #endif
 
   m_dirname = std::string( dirname_c_str );
   m_write_snapshot_mat = (std::string(write_snapshot_mat) == "true");
 
   if (m_num_window_samples <= m_rdim) {
     printf("ERROR:DMDROMObject::DMDROMObject() - m_num_window_samples <= m_rdim!!");
   }
 
   m_tic = 0;
   m_curr_win = 0;
 }

◆ ~DMDROMObject()

virtual ~DMDROMObject ( )

inlinevirtual

Destructor

Definition at line 51 of file rom_object_dmd.h.

     {
       for (int i = 0; i < m_dmd.size(); i++) delete m_dmd[i];
       m_dmd.clear();
       m_intervals.clear();
       m_dmd_is_trained.clear();
     }

Member Function Documentation

◆ projectInitialSolution()

virtual void projectInitialSolution ( CAROM::Vector & a_U )

inlinevirtual

Project initial solution for prediction

Parameters

a_U	solution vector

Implements ROMObject.

Definition at line 60 of file rom_object_dmd.h.

     {
       if (m_dmd.size() == 0) {
         if (!m_rank) {
           printf("ERROR in DMDROMObject::projectInitialSolution() - m_dmd is a vector of size 0.\n");
         }
         return;
       }
 
       m_dmd[0]->projectInitialCondition( &a_U );
       for (int i = 1; i < m_dmd.size(); i++) {
         m_dmd[i]->projectInitialCondition( m_dmd[i-1]->predict(m_intervals[i].first) );
       }
       return;
     }

◆ takeSample()

void takeSample	(	const CAROM::Vector &	a_U,
		const double	a_time
	)

virtual

take a sample (solution snapshot)

Parameters

a_U	solution vector
a_time	sample time

Implements ROMObject.

Definition at line 126 of file DMDROMObject.cpp.

 {
   if (m_tic == 0) {
 
     m_dmd.push_back( new CAROM::DMD(m_vec_size, m_dt) );
     m_dmd_is_trained.push_back(false);
     m_intervals.push_back( Interval(a_time, m_t_final) );
 
     if (!m_rank) {
       printf( "DMDROMObject::takeSample() - creating new DMD object for sim. domain %d, var %d, t=%f (total: %d).\n",
               m_sim_idx, m_var_idx, m_intervals[m_curr_win].first, m_dmd.size());
     }
     m_dmd[m_curr_win]->takeSample( a_U.getData(), a_time );
 
   } else {
 
     m_dmd[m_curr_win]->takeSample( a_U.getData(), a_time );
 
     if (m_tic%m_num_window_samples == 0) {
 
       m_intervals[m_curr_win].second = a_time;
       int ncol = m_dmd[m_curr_win]->getSnapshotMatrix()->numColumns();
       if (!m_rank) {
         printf( "DMDROMObject::train() - training DMD object %d for sim. domain %d, var %d with %d samples.\n", 
                 m_curr_win, m_sim_idx, m_var_idx, ncol );
       }
 
       if (m_write_snapshot_mat) {
         char idx_string[_MAX_STRING_SIZE_];
         sprintf(idx_string, "%04d", m_curr_win);
         std::string fname_root(m_dirname + "/snapshot_mat_"+std::string(idx_string));
         m_dmd[m_curr_win]->getSnapshotMatrix()->write(fname_root);
       }
       m_dmd[m_curr_win]->train(m_rdim);
       m_dmd_is_trained[m_curr_win] = true;
 
       m_curr_win++;
 
       m_dmd.push_back( new CAROM::DMD(m_vec_size, m_dt) );
       m_dmd_is_trained.push_back(false);
       m_intervals.push_back( Interval(a_time, m_t_final) );
       m_dmd[m_curr_win]->takeSample( a_U.getData(), a_time );
       if (!m_rank) {
         printf("DMDROMObject::takeSample() - creating new DMD object for sim. domain %d, var %d, t=%f (total: %d).\n",
                m_sim_idx, m_var_idx, m_intervals[m_curr_win].first, m_dmd.size());
       }
     }
 
   }
 
   m_tic++;
   return;
 }

◆ train()

void train ( )

virtual

train the DMD object

train the DMD objects

Implements ROMObject.

Definition at line 182 of file DMDROMObject.cpp.

 {
   /* make sure the number of columns for the last DMD isn't less than m_rdim */
   {
     int last_win = m_dmd.size() - 1;
     int num_columns( m_dmd[last_win]->getSnapshotMatrix()->numColumns() );
     if (num_columns <= m_rdim) {
       m_dmd.pop_back();
       m_intervals.pop_back();
       m_intervals[m_intervals.size()-1].second = m_t_final;
       if (!m_rank) {
         printf("DMDROMObject::train() - last window DMD for sim. domain %d, var %d has %d sample(s) only; deleted it ",
                m_sim_idx,
                m_var_idx,
                num_columns );
         printf("(total: %d).\n", m_dmd.size());
       }
     }
   }
 
   if (m_dmd.size() > 0) {
     for (int i = 0; i < m_dmd.size(); i++) {
       if (!m_dmd_is_trained[i]) {
         int ncol = m_dmd[i]->getSnapshotMatrix()->numColumns();
         if (!m_rank) {
           printf( "DMDROMObject::train() - training DMD object %d for sim. domain %d, var %d with %d samples.\n", 
                   m_curr_win, m_sim_idx, m_var_idx, ncol );
         }
         if (m_write_snapshot_mat) {
           char idx_string[_MAX_STRING_SIZE_];
           sprintf(idx_string, "%04d", i);
           std::string fname_root(m_dirname + "/snapshot_mat_"+std::string(idx_string));
           m_dmd[i]->getSnapshotMatrix()->write(fname_root);
         }
         m_dmd[i]->train(m_rdim);
       }
     }
   } else {
     printf("ERROR in DMDROMObject::train(): m_dmd is of size zero!");
   }
 
   return;
 }

◆ predict()

virtual const CAROM::Vector* const predict ( const double a_t ) const

inlinevirtual

compute prediction at given time

Parameters

a_t	time at which to predict solution

Implements ROMObject.

Definition at line 84 of file rom_object_dmd.h.

     {
       for (int i = 0; i < m_dmd.size(); i++) {
         if (   (a_t >= m_intervals[i].first) 
             && (  (a_t < m_intervals[i].second) || (m_intervals[i].second < 0)  ) ){
           return m_dmd[i]->predict(a_t);
         }
       }
       printf("ERROR in DMDROMObject::predict(): m_dmd is of size zero or interval not found!");
       return nullptr;
     }

◆ save()

void save ( const std::string & a_fname_root ) const

virtual

save DMD object to file

Save DMD objects to file: the DMD object files will be saved in the subdirectory with the name DMDROMObject::m_dirname. They are in a format that libROM can read from.

Note: If the subdirectory DMDROMObject::m_dirname does not exist, the DMD objects will not be written (even though the screen output will claim they were written)!. The code may not report any error, or one may see HDF5 file writing errors.

Parameters

a_fname_root Filename root

Implements ROMObject.

Definition at line 235 of file DMDROMObject.cpp.

 {
   std::string fname_root = m_dirname + "/";
   std::string summary_fname_root = m_dirname + "/";
   std::string header_fname = m_dirname + "/";
   if (a_fname_root == "") {
     fname_root += "dmdobj_";
     summary_fname_root += "dmd_summary_";
     header_fname += "dmd_header.dat";
   } else {
     fname_root += (a_fname_root+"_dmdobj_");
     summary_fname_root += (a_fname_root+"_dmd_summary_");
     header_fname += (a_fname_root+"_dmd_header.dat");
   }
 
   if (!m_rank) {
     FILE* out;
     out = fopen(header_fname.c_str(), "w");
     fprintf(out, "%d\n", m_dmd.size());
     for (int i = 0; i < m_dmd.size(); i++) {
       fprintf(out, "%1.16e %1.16e\n", m_intervals[i].first, m_intervals[i].second);
     }
     fclose(out);
   }
 
   for (int i = 0; i < m_dmd.size(); i++) {
     char idx_string[_MAX_STRING_SIZE_];
     sprintf(idx_string, "%04d", i);
     std::string fname = fname_root + std::string(idx_string);
     std::string summary_fname = summary_fname_root + std::string(idx_string);
     if (!m_rank) {
       printf( "  Saving DMD object and summary (%s, %s).\n", 
               fname.c_str(), summary_fname.c_str() );
     }
     m_dmd[i]->save(fname);
     //m_dmd[i]->summary(summary_fname);
   }
 
   return;
 }

◆ load()

void load ( const std::string & a_fname_root )

virtual

load DMD object from file

Load DMD objects from file: the DMD object files must be in the subdirectory with the name DMDROMObject::m_dirname. They must be in a format that libROM can read from. The number of objects and their parameters must correspond to the simulation being run for which this object is defined.

Parameters

a_fname_root Filename root

Implements ROMObject.

Definition at line 281 of file DMDROMObject.cpp.

 {
   std::string fname_root = m_dirname + "/";
   std::string header_fname = m_dirname + "/";
   if (a_fname_root == "") {
     fname_root += "dmdobj_";
     header_fname += "dmd_header.dat";
   } else {
     fname_root += (a_fname_root+"_dmdobj_");
     header_fname += (a_fname_root+"_dmd_header.dat");
   }
 
   int num_dmds = 0;
   std::vector<double> intervals_start(0), intervals_end(0);
   if (!m_rank) {
     FILE* in;
     in = fopen(header_fname.c_str(), "r");
     fscanf(in, "%d\n", &num_dmds);
     intervals_start.resize(num_dmds);
     intervals_end.resize(num_dmds);
     for (int i = 0; i < num_dmds; i++) {
       fscanf(in, "%lf", &intervals_start[i]);
       fscanf(in, "%lf", &intervals_end[i]);
     }
     fclose(in);
   }
 #ifndef serial
   MPI_Bcast(&num_dmds,1,MPI_INT,0,MPI_COMM_WORLD);
   if (m_rank) {
     intervals_start.resize(num_dmds);
     intervals_end.resize(num_dmds);
   }
   MPI_Bcast(intervals_start.data(),num_dmds,MPI_DOUBLE,0,MPI_COMM_WORLD);
   MPI_Bcast(intervals_end.data(),num_dmds,MPI_DOUBLE,0,MPI_COMM_WORLD);
 #endif
 
   for (int i = 0; i < num_dmds; i++) {
 
     m_intervals.push_back( Interval(intervals_start[i],intervals_end[i]) );
 
     char idx_string[_MAX_STRING_SIZE_];
     sprintf(idx_string, "%04d", i);
     std::string fname = fname_root + std::string(idx_string);
     if (!m_rank) {
       printf( "  Loading DMD object (%s), time window=[%1.2e,%1.2e].\n", 
               fname.c_str(),
               m_intervals[i].first, m_intervals[i].second );
     }
     m_dmd.push_back( new CAROM::DMD(fname) );
   }
 
   return;
 }

Field Documentation

◆ m_dmd

std::vector<CAROM::DMD*> m_dmd

protected

Vector of DMD objects

Definition at line 104 of file rom_object_dmd.h.

◆ m_dmd_is_trained

std::vector<bool> m_dmd_is_trained

protected

Flag to indicate if DMD is trained

Definition at line 105 of file rom_object_dmd.h.

◆ m_intervals

std::vector<Interval> m_intervals

protected

Time intervals for each DMD object

Definition at line 106 of file rom_object_dmd.h.

◆ m_rank

int m_rank

protected

MPI rank

Definition at line 108 of file rom_object_dmd.h.

◆ m_nproc

int m_nproc

protected

Number of MPI ranks

Definition at line 109 of file rom_object_dmd.h.

◆ m_vec_size

int m_vec_size

protected

Local size of solution vector

Definition at line 111 of file rom_object_dmd.h.

◆ m_dt

double m_dt

protected

Time step size

Definition at line 112 of file rom_object_dmd.h.

◆ m_t_final

double m_t_final

protected

Final time

Definition at line 113 of file rom_object_dmd.h.

◆ m_rdim

int m_rdim

protected

Latent space dimension

Definition at line 114 of file rom_object_dmd.h.

◆ m_num_window_samples

int m_num_window_samples

protected

Number of samples per DMD for time-windowing

Definition at line 116 of file rom_object_dmd.h.

◆ m_tic

int m_tic

protected

private ticker to count number of samples taken

Definition at line 118 of file rom_object_dmd.h.

◆ m_curr_win

int m_curr_win

protected

private index for current window

Definition at line 119 of file rom_object_dmd.h.

◆ m_sim_idx

int m_sim_idx

protected

simulation index of this object for ensemble simulations

Definition at line 121 of file rom_object_dmd.h.

◆ m_var_idx

int m_var_idx

protected

component index of this object if component-wise ROMs are being used

Definition at line 122 of file rom_object_dmd.h.

◆ m_write_snapshot_mat

bool m_write_snapshot_mat

protected

Write snapshot matrix to file or not

Definition at line 124 of file rom_object_dmd.h.

◆ m_dirname

std::string m_dirname

protected

Subdirectory where DMD objects are written to or read from

Definition at line 126 of file rom_object_dmd.h.

The documentation for this class was generated from the following files:

include/rom_object_dmd.h
src/ROM/DMDROMObject.cpp

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DMDROMObject()

◆ ~DMDROMObject()

Member Function Documentation

◆ projectInitialSolution()

◆ takeSample()

◆ train()

◆ predict()

◆ save()

◆ load()

Field Documentation

◆ m_dmd

◆ m_dmd_is_trained

◆ m_intervals

◆ m_rank

◆ m_nproc

◆ m_vec_size

◆ m_dt

◆ m_t_final

◆ m_rdim

◆ m_num_window_samples

◆ m_tic

◆ m_curr_win

◆ m_sim_idx

◆ m_var_idx

◆ m_write_snapshot_mat

◆ m_dirname