Toru Shiozaki
Toru Shiozaki
I know this is not a bug but has to be reported. The MP2 assembly step becomes slower by a lot if I use BTAS' native operator(). (I previously wrote...
I only see the following. I assume that we need to explicitly define this (have you guys discussed before?). Important because one may want to implement a custom storage class,...
Recent changes change the behavior of storage and there is a bug in BAGEL (I am wasting a lot of time now; we are sure that code mentioned in #63...
double gemm with alpha = 1.0 and beta = 2 would not be accepted. Also types are handled in a very error-prone way because everything is passed around using const...
The following code returns 1.0, which is wrong. ``` #include #include #include using namespace std; using namespace btas; int main() { Range b(5); Tensor a(b); a(0) = 1.0; a(1) =...
I hit a bug that most likely originates from a bug in btas::group for TensorView. Here is my code. @emstoudenmire Do you have any ideas? ``` // slowest index of...
As long as it is specified by initialier_list of constexpr int, there must be a way in theory (In practice I do not know). If it would be very powerful...
I think it is _very_ useful for optimization of contract if one could obtain information on rank() and contiguous() at compile time. Otherwise there will be a bunch of if...
What was the consensus on the complex conjugate? In other words do we allow the users to specify complex conjugation in the interface? As far as I see there is...
btas::contract(1.0, *this, {0,1,3}, *cmat, {3,2}, 0.0, *out, {0,1,2}); with the column major should not permute indices, but it does. Not sure if this is a problem in the index notation,...