nutbar/NuBasis.hh at master · ragnarstroberg/nutbar · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
#ifndef NuBasis_h

#include <iostream>
#include <fstream>
#include <stdlib.h>
#include <stdint.h>
#include <iomanip>
#include <vector>
#include <string>
#include <cstdint>
#include <bitset>

#define KEY_BITS 64

using namespace std;

typedef int16_t part_type;
typedef int64_t mvec_type;
//typedef vector<mvec_type> key_type;
typedef bitset<KEY_BITS> key_type;


struct MschemeOrbit;

class NuBasis
{
 public:

  int32_t no_spart,sum_nJT;
  vector<int32_t> pindx,nJT,j,t,ibf,max_cM2,max_cTz2;
  vector<vector<part_type>> partition;
  vector<vector<key_type>> vec;
  static const int gwords;
  vector<MschemeOrbit> m_orbits;

  NuBasis();
  void ReadFile(string fname);
  void ReadSPS(string fname);
  void PrintBasis();
  void PrintSPS();
  void Clear();

  //TODO: These should be applied to a J state
//  void LoweringOperator( key_type& mvec_in, vector<key_type>& mvec_out, vector<float>& coefs);
//  void RaisingOperator( key_type& mvec_in, vector<key_type>& mvec_out, vector<float>& coefs);
//  void LoweringOperator( key_type& mvec_in, vector<key_type>& mvec_out, vector<float>& coefs);
//  void RaisingOperator( key_type& mvec_in, vector<key_type>& mvec_out, vector<float>& coefs);
//  key_type TimeReverse( key_type& mvec_in);
//  key_type TimeReverse( key_type& mvec_in);

};


struct MschemeOrbit
{
  int n,l2,j2,mj2,tz2;

  MschemeOrbit(int n, int l2, int j2, int mj2, int tz2): n(n),l2(l2),j2(j2),mj2(mj2),tz2(tz2){};
};


#define NuBasis_h
#endif