3j symbol in progress. Very inelegant, but good enough for now.
Merge branch 'master' of github.com:shiozaki/BAGEL into ecp

# By Toru Shiozaki (12) and Shane Parker (6)
# Via Toru Shiozaki
* 'master' of github.com:shiozaki/BAGEL:
  seems synchronization is a real issue. Also not sure why I declared matrices in Davidson to be localized. At least it is better as diagonalization is always sync'ed in this way
  to compile with -O3 -Wunused-variable
  -Wunused-variable for parallel builds
  More strict warning checks with -Wunused -Wnounsed-function
  Resolves #30
  davidson related cleanup
  tentatively disables SuperCI tests
  removing orthog in the Davidson context. normalization is still necessary, though.
  removing orthog
  davidson issue is most likely fixed. Convergence is virtually the same even if the space is reduced during the iteration.
  remove debug printing
  Revamped davidson that does not assume orthogonality anywhere. Test cases pass, but I’m still pretty sure there is a bug somewhere.
  zquartev was returning wrong eigenvectors. fixed.
  fixed the code to elminate warnings by clang (one was clearly a bug).
  cleanup
  Avoiding sigma formation in MEH iterations.
  SparseMatrix::print_table for debugging purposes
  bug fixes
Wigner 3j symbols up to j = 6.
Merge branch 'master' of github.com:shiozaki/BAGEL into ecp

# By Shane Parker
# Via Shane Parker
* 'master' of github.com:shiozaki/BAGEL:
  circumventing buggy code
  Move spin_decontaminate out to the whole Dvec (so it can be parallelized).
  When building Gamma matrices in MEH-RAS, tasks now look ahead to check whether terms will contribute BEFORE computing them. Also reorganizing code base so MEHRAS and MEHDistRAS use the same “next block” function.
  bug fixes and cleanup in MEH
  fix for using Muffle in parallel
compute the modified bessel function of the first kind for the angular integrals
compute the spherical harmonics.
still working on the angular integrals
expand a real spherical harmonics in terms of xyz, working for |00> and |10>
one bug fixed but still not done yet
an attempt to tabulate the carsph coefficients.
changed my mind. I think I can just use carsphlist and all the _carsph files.
Implemented the angular integral. To be checked…
spherical harmonics to cartesian, checked.
Merge branch 'integral' of github.com:shiozaki/BAGEL into ecp

# By Toru Shiozaki (99) and others
# Via Ryan Dave Reynolds (8) and Toru Shiozaki (6)
* 'integral' of github.com:shiozaki/BAGEL: (157 commits)
  some cleanup
  bug fix in testing interface
  Hessian works with kramers pairs.
  attempt to do G(1,2), but totally wrong
  minor fixes
  intermediate checkin - still wrong
  gradients fine with kramers on
  support functions for active indices
  reproduced the closed-virtual diagonal hessian
  numerical derivatives of gradient and hessian. Rotation angle can be controlled from the input (using "debugrot" keyword).
  moving things from Jeb's branch. Still under way
  easier construction of StaticDist with alignments
  allowing custom alignments in DFDistT via an optionally passing StaticDist.
  minor improvements
  zfci restartable
  (1) boost::serialization::make_array does not bit-wise serialize std::complex<double>. Made a wrapper for that; (2) dvec needed some work
  More tools for Resolution of Identity testing
  Interface London integral codes to MixedBasis class
  preparing zfci for serialization
  hz restartable.
  ...
Merge branch 'master' of github.com:shiozaki/BAGEL into ecp

# By Ryan Dave Reynolds
# Via Ryan Dave Reynolds (1) and Toru Shiozaki (1)
* 'master' of github.com:shiozaki/BAGEL:
  Quick fixes (after pull request feedback)
Merge branch 'master' of github.com:shiozaki/BAGEL into ecp

# By Shane Parker (2) and Toru Shiozaki (2)
# Via Shane Parker
* 'master' of github.com:shiozaki/BAGEL:
  missed one
  Restructuring localization classes so that arbitrary subspaces can be localized using the input keyword "subspaces." This can be helpful when trying to choose active orbitals among virtuals. Additionally, changing definition of metric in Pipek Mezey localization to account for number of orbitals. Now it is more like an RMS.
  cosmetic changes in configure.ac
  Hopefully resolves #32
Evaluate the angular integral (2SHs and 1USP) by transforming each SH into a USP.
Bugs fixed
update the projection integral
part of testing the projection integral