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 #define MPICH_SKIP_MPICXX
00036 #include <mpi.h>
00037 #include <util/ref/ref.h>
00038 #include <util/group/memory.h>
00039 #include <chemistry/qc/mbptr12/r12ia.h>
00040
00041 namespace sc {
00042
00044
00045
00046
00047
00048
00049
00050
00051
00052 class R12IntsAcc_MPIIOFile: public R12IntsAcc {
00053
00054 protected:
00055 Ref<MemoryGrp> mem_;
00056 int nproc_;
00057 size_t nints_per_block_;
00058 int icounter_;
00059 char *filename_;
00060 MPI_File datafile_;
00061
00062 struct PairBlkInfo {
00063 double* ints_[max_num_te_types_];
00064 int refcount_[max_num_te_types_];
00065 MPI_Offset offset_;
00066 } *pairblk_;
00067
00069 void check_error_code_(int errcod) const;
00070
00071 public:
00072 R12IntsAcc_MPIIOFile(Ref<MemoryGrp>& mem, const char *filename, int num_te_types, int nbasis1, int nbasis2,
00073 int nocc, int nfzc, bool restart);
00074 ~R12IntsAcc_MPIIOFile();
00075
00077 void store_pair_block(int i, int j, double *ints);
00080 void commit();
00083 void deactivate();
00085 void release_pair_block(int i, int j, tbint_type oper_type);
00087 bool is_local(int i, int j) const { return true;};
00089 bool is_avail(int i, int j) const { return true;};
00091 bool has_access(int proc) const { return true;};
00093 bool can_restart() const { return true; };
00094
00095
00096 int ij_index(int i, int j) const { return i*nocc_act_ + j; };
00097 };
00098
00100
00101
00102
00103
00104
00105
00106
00107 class R12IntsAcc_MPIIOFile_Ind: public R12IntsAcc_MPIIOFile {
00108
00109 public:
00110 R12IntsAcc_MPIIOFile_Ind(Ref<MemoryGrp>& mem, const char *filename, int num_te_types, int nbasis1, int nbasis2,
00111 int nocc, int nfzc, bool restart);
00112 ~R12IntsAcc_MPIIOFile_Ind();
00113
00119 void store_memorygrp(Ref<MemoryGrp>& mem, int ni);
00121 double* retrieve_pair_block(int i, int j, tbint_type oper_type);
00122 };
00123
00124 }
00125
00126 #endif
00127
00128
00129
00130
00131
Generated at Sat Dec 18 15:14:22 2004 for MPQC
2.2.3 using the documentation package Doxygen
1.3.7-20040617.