r12ia_mpiiofile.h
00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #ifndef _chemistry_qc_mbptr12_r12ia_mpiiofile_h
00029 #define _chemistry_qc_mbptr12_r12ia_mpiiofile_h
00030
00031 #ifdef __GNUC__
00032 #pragma interface
00033 #endif
00034
00035 #include <mpi.h>
00036 #include <util/ref/ref.h>
00037 #include <util/group/memory.h>
00038 #include <chemistry/qc/mbptr12/r12ia.h>
00039
00040 namespace sc {
00041
00043
00044
00045
00046
00047
00048
00049
00050
00051 class R12IntsAcc_MPIIOFile: public R12IntsAcc {
00052
00053 protected:
00054 Ref<MemoryGrp> mem_;
00055 int nproc_;
00056 size_t nints_per_block_;
00057 int icounter_;
00058 char *filename_;
00059 MPI_File datafile_;
00060
00061 struct PairBlkInfo {
00062 double* ints_[max_num_te_types_];
00063 int refcount_[max_num_te_types_];
00064 MPI_Offset offset_;
00065 } *pairblk_;
00066
00068 void check_error_code_(int errcod) const;
00069
00070 public:
00071 R12IntsAcc_MPIIOFile(Ref<MemoryGrp>& mem, const char *filename, int num_te_types, int nbasis1, int nbasis2,
00072 int nocc, int nfzc, bool restart);
00073 ~R12IntsAcc_MPIIOFile();
00074
00076 void store_pair_block(int i, int j, double *ints);
00079 void commit();
00082 void deactivate();
00084 void release_pair_block(int i, int j, tbint_type oper_type);
00086 bool is_local(int i, int j) const { return true;};
00088 bool is_avail(int i, int j) const { return true;};
00090 bool has_access(int proc) const { return true;};
00092 bool can_restart() const { return true; };
00093
00094
00095 int ij_index(int i, int j) const { return i*nocc_act_ + j; };
00096 };
00097
00099
00100
00101
00102
00103
00104
00105
00106 class R12IntsAcc_MPIIOFile_Ind: public R12IntsAcc_MPIIOFile {
00107
00108 public:
00109 R12IntsAcc_MPIIOFile_Ind(Ref<MemoryGrp>& mem, const char *filename, int num_te_types, int nbasis1, int nbasis2,
00110 int nocc, int nfzc, bool restart);
00111 ~R12IntsAcc_MPIIOFile_Ind();
00112
00118 void store_memorygrp(Ref<MemoryGrp>& mem, int ni);
00120 double* retrieve_pair_block(int i, int j, tbint_type oper_type);
00121 };
00122
00123 }
00124
00125 #endif
00126
00127
00128
00129
00130
Generated at Thu Jul 29 19:53:51 2004 for MPQC
2.2.1 using the documentation package Doxygen
1.3.7.