mirror of
https://github.com/triqs/dft_tools
synced 2024-12-25 05:43:40 +01:00
Add test for many_body_operators
This commit is contained in:
parent
422cabd52e
commit
89de284cc8
2
test/triqs/operators/CMakeLists.txt
Normal file
2
test/triqs/operators/CMakeLists.txt
Normal file
@ -0,0 +1,2 @@
|
|||||||
|
all_tests()
|
||||||
|
|
100
test/triqs/operators/operator_test.cpp
Normal file
100
test/triqs/operators/operator_test.cpp
Normal file
@ -0,0 +1,100 @@
|
|||||||
|
/*******************************************************************************
|
||||||
|
*
|
||||||
|
* TRIQS: a Toolbox for Research in Interacting Quantum Systems
|
||||||
|
*
|
||||||
|
* Copyright (C) 2013 by I. Krivenko
|
||||||
|
*
|
||||||
|
* TRIQS is free software: you can redistribute it and/or modify it under the
|
||||||
|
* terms of the GNU General Public License as published by the Free Software
|
||||||
|
* Foundation, either version 3 of the License, or (at your option) any later
|
||||||
|
* version.
|
||||||
|
*
|
||||||
|
* TRIQS is distributed in the hope that it will be useful, but WITHOUT ANY
|
||||||
|
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
|
||||||
|
* FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
|
||||||
|
* details.
|
||||||
|
*
|
||||||
|
* You should have received a copy of the GNU General Public License along with
|
||||||
|
* TRIQS. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
*
|
||||||
|
******************************************************************************/
|
||||||
|
#include <triqs/utility/first_include.hpp>
|
||||||
|
#include <triqs/operators/many_body_operator.hpp>
|
||||||
|
#include <iostream>
|
||||||
|
#include <vector>
|
||||||
|
#include <boost/archive/text_oarchive.hpp>
|
||||||
|
#include <boost/archive/text_iarchive.hpp>
|
||||||
|
#include <sstream>
|
||||||
|
|
||||||
|
|
||||||
|
using namespace triqs::utility;
|
||||||
|
|
||||||
|
int main(int argc, char **argv)
|
||||||
|
{
|
||||||
|
// Operators without indices
|
||||||
|
auto op_with_no_indices = c() + c_dag() - n();
|
||||||
|
std::cout << "op_with_no_indices = " << op_with_no_indices << std::endl;
|
||||||
|
|
||||||
|
// Operators with many indices
|
||||||
|
auto op_with_many_indices = c(1,2,"a",true,-2) +
|
||||||
|
c_dag(3,15,"b",false,-5);
|
||||||
|
std::cout << "op_with_many_indices = " << op_with_many_indices << std::endl;
|
||||||
|
|
||||||
|
// Commutation relations
|
||||||
|
std::vector<many_body_operator<double> > C = {c(1), c(2), c(3)};
|
||||||
|
std::vector<many_body_operator<double> > Cd = {c_dag(1), c_dag(2), c_dag(3)};
|
||||||
|
|
||||||
|
std::cout << std::endl << "Anticommutators:" << std::endl;
|
||||||
|
for(auto const& cdi : Cd)
|
||||||
|
for(auto const& ci : C){
|
||||||
|
std::cout << "{" << cdi << ", " << ci << "} = " << cdi*ci + ci*cdi << std::endl;
|
||||||
|
}
|
||||||
|
|
||||||
|
std::cout << std::endl << "Commutators:" << std::endl;
|
||||||
|
for(auto const& cdi : Cd)
|
||||||
|
for(auto const& ci : C){
|
||||||
|
std::cout << "[" << cdi << ", " << ci << "] = " << cdi*ci - ci*cdi << std::endl;
|
||||||
|
}
|
||||||
|
|
||||||
|
// Algebra
|
||||||
|
auto x = c(0);
|
||||||
|
auto y = c_dag(1);
|
||||||
|
|
||||||
|
std::cout << std::endl << "Algebra:" << std::endl;
|
||||||
|
std::cout << "x = " << x << std::endl;
|
||||||
|
std::cout << "y = " << y << std::endl;
|
||||||
|
|
||||||
|
std::cout << "-x = " << -x << std::endl;
|
||||||
|
std::cout << "x + 2.0 = " << x + 2.0 << std::endl;
|
||||||
|
std::cout << "2.0 + x = " << 2.0 + x << std::endl;
|
||||||
|
std::cout << "x - 2.0 = " << x - 2.0 << std::endl;
|
||||||
|
std::cout << "2.0 - x = " << 2.0 - x << std::endl;
|
||||||
|
std::cout << "3.0*y = " << 3.0*y << std::endl;
|
||||||
|
std::cout << "y*3.0 = " << y*3.0 << std::endl;
|
||||||
|
std::cout << "x + y = " << x + y << std::endl;
|
||||||
|
std::cout << "x - y = " << x - y << std::endl;
|
||||||
|
std::cout << "(x + y)*(x - y) = " << (x + y)*(x - y) << std::endl;
|
||||||
|
|
||||||
|
// N^3
|
||||||
|
std::cout << std::endl << "N^3:" << std::endl;
|
||||||
|
auto N = n("up") + n("dn");
|
||||||
|
auto N3 = N*N*N;
|
||||||
|
std::cout << "N = " << N << std::endl;
|
||||||
|
std::cout << "N^3 = " << N3 << std::endl;
|
||||||
|
|
||||||
|
// Serialization
|
||||||
|
std::stringstream ss;
|
||||||
|
boost::archive::text_oarchive oa(ss);
|
||||||
|
oa & N3;
|
||||||
|
|
||||||
|
boost::archive::text_iarchive ia(ss);
|
||||||
|
many_body_operator<double> new_N3;
|
||||||
|
ia & new_N3;
|
||||||
|
|
||||||
|
std::cout << "New N^3 = " << new_N3 << std::endl;
|
||||||
|
|
||||||
|
auto X = c_dag(1) * c_dag(2) * c(3) * c(4);
|
||||||
|
std::cout << "X = "<< X<<std::endl;
|
||||||
|
std::cout << "dagger(X) = "<< dagger(X)<<std::endl;
|
||||||
|
return 0;
|
||||||
|
}
|
45
test/triqs/operators/operator_test.output
Normal file
45
test/triqs/operators/operator_test.output
Normal file
@ -0,0 +1,45 @@
|
|||||||
|
op_with_no_indices = 1*C^+() + 1*C() + -1*C^+()C()
|
||||||
|
op_with_many_indices = 1*C^+(3,15,b,0,-5) + 1*C(1,2,a,1,-2)
|
||||||
|
|
||||||
|
Anticommutators:
|
||||||
|
{1*C^+(1), 1*C(1)} = 1
|
||||||
|
{1*C^+(1), 1*C(2)} = 0
|
||||||
|
{1*C^+(1), 1*C(3)} = 0
|
||||||
|
{1*C^+(2), 1*C(1)} = 0
|
||||||
|
{1*C^+(2), 1*C(2)} = 1
|
||||||
|
{1*C^+(2), 1*C(3)} = 0
|
||||||
|
{1*C^+(3), 1*C(1)} = 0
|
||||||
|
{1*C^+(3), 1*C(2)} = 0
|
||||||
|
{1*C^+(3), 1*C(3)} = 1
|
||||||
|
|
||||||
|
Commutators:
|
||||||
|
[1*C^+(1), 1*C(1)] = -1 + 2*C^+(1)C(1)
|
||||||
|
[1*C^+(1), 1*C(2)] = 2*C^+(1)C(2)
|
||||||
|
[1*C^+(1), 1*C(3)] = 2*C^+(1)C(3)
|
||||||
|
[1*C^+(2), 1*C(1)] = 2*C^+(2)C(1)
|
||||||
|
[1*C^+(2), 1*C(2)] = -1 + 2*C^+(2)C(2)
|
||||||
|
[1*C^+(2), 1*C(3)] = 2*C^+(2)C(3)
|
||||||
|
[1*C^+(3), 1*C(1)] = 2*C^+(3)C(1)
|
||||||
|
[1*C^+(3), 1*C(2)] = 2*C^+(3)C(2)
|
||||||
|
[1*C^+(3), 1*C(3)] = -1 + 2*C^+(3)C(3)
|
||||||
|
|
||||||
|
Algebra:
|
||||||
|
x = 1*C(0)
|
||||||
|
y = 1*C^+(1)
|
||||||
|
-x = -1*C(0)
|
||||||
|
x + 2.0 = 2 + 1*C(0)
|
||||||
|
2.0 + x = 2 + 1*C(0)
|
||||||
|
x - 2.0 = -2 + 1*C(0)
|
||||||
|
2.0 - x = 2 + -1*C(0)
|
||||||
|
3.0*y = 3*C^+(1)
|
||||||
|
y*3.0 = 3*C^+(1)
|
||||||
|
x + y = 1*C^+(1) + 1*C(0)
|
||||||
|
x - y = -1*C^+(1) + 1*C(0)
|
||||||
|
(x + y)*(x - y) = 2*C^+(1)C(0)
|
||||||
|
|
||||||
|
N^3:
|
||||||
|
N = 1*C^+(dn)C(dn) + 1*C^+(up)C(up)
|
||||||
|
N^3 = 1*C^+(dn)C(dn) + 1*C^+(up)C(up) + 6*C^+(dn)C^+(up)C(up)C(dn)
|
||||||
|
New N^3 = 1*C^+(dn)C(dn) + 1*C^+(up)C(up) + 6*C^+(dn)C^+(up)C(up)C(dn)
|
||||||
|
X = -1*C^+(1)C^+(2)C(4)C(3)
|
||||||
|
dagger(X) = -1*C^+(3)C^+(4)C(2)C(1)
|
Loading…
Reference in New Issue
Block a user