def test_xplib(self): # pass # self.assertEqual(funcs.tupleProduct((2, 3)), 6) CTL.setXP(None) with self.assertRaises(AttributeError) as context: simplestExample() CTL.setXP(np) res = simplestExample() print('res = {}'.format(res))
def test_cupy(self): # pass # self.assertEqual(funcs.tupleProduct((2, 3)), 6) try: import cupy as cp except: return CTL.setXP(cp) res = simplestExample() print('res = {}'.format(res))
def main(argv): if len(argv) < 2: print('Usage: ./CORONER [input_file]') sys.exit(1) inpoot_file = argv[1] ast = parse_input(inpoot_file) print("\nCTL ALGORITHM: SERIAL") CTL.execute_ctl_serial(ast) print("\nCTL ALGORITHM: PARALLEL") CTL.execute_ctl_parallel(ast) print("\nMARK SWEEP ALGORITHM") MARKSWEEP.mark_sweep(ast) sys.exit(0)
# -*- coding: utf-8 -*- """ Created on Thu Apr 16 20:52:12 2015 @author: Travis """ from numpy import * import CTL ### # Runtime ### recordList = ["31"] order = 10 #lowPassFilter = [None, [order, 0.1], [order, 0.25], [order, 0.5]] lowPassFilter = [[order, 0.05]] for lpf in lowPassFilter: CTL.fullAnalysis(recordList, 'physionetdata/', lpf, 1)
# -*- coding: utf-8 -*- """ Created on Thu Apr 16 20:52:12 2015 @author: Travis """ from numpy import * import CTL ### # Runtime ### # Had to lower embedding dimension and N due to small datasets embeddingDim = 10 N = 2 recordList = ["SockeyeSalmon", "ElephantSeal", "ChumSalmon", "Beaver"] order = 2 lowPassFilter = [None, [order, 0.1], [order, 0.25], [order, 0.5]] for lpf in lowPassFilter: CTL.fullAnalysis(recordList, 'populationdata/', lpf, 0, embeddingDim, N)
# -*- coding: utf-8 -*- """ Created on Thu Apr 16 20:52:12 2015 @author: Travis """ from numpy import * import CTL ### # Runtime ### # Had to lower embedding dimension and N due to small datasets embeddingDim = 10 N = 2 recordList = ["SockeyeSalmon", "ElephantSeal", "ChumSalmon", "Beaver"] order = 2 lowPassFilter = [None, [order, 0.1], [order, 0.25], [order, 0.5]] for lpf in lowPassFilter: CTL.fullAnalysis(recordList, "populationdata/", lpf, 0, embeddingDim, N)
parentdir = os.path.dirname(currentdir) sys.path.append(os.path.join(parentdir, 'src')) from CTL.tensor.tensor import Tensor from CTL.tensor.contract.link import makeLink from CTL.tensor.contract.optimalContract import contractAndCostWithSequence from CTL.examples.MPS import contractWithMPS from CTL.models.Ising import IsingSiteTensor, IsingEdgeMatrix, IsingTNFromUndirectedGraph, exactZFromGraphIsing from CTL.funcs.graphFuncs import squareLatticeFBC, squareLatticeFBC import numpy as np import cupy as cp import CTL CTL.setXP(cp) def contractHandmadeTN(): print('contractHandmadeTN():') a = Tensor(shape=(3, 5, 7), labels=['a3', 'a5', 'a7']) b = Tensor(shape=(2, 4, 5), labels=['b2', 'b4', 'b5']) c = Tensor(shape=(2, 7, 7, 7), labels=['c2', 'c71', 'c72', 'c73']) d = Tensor(shape=(7, 7, 3, 4), labels=['d71', 'd72', 'd3', 'd4']) e = Tensor(shape=(3, 3, 5), labels=['e31', 'e32', 'e5']) f = Tensor(shape=(2, 2, 5), labels=['f21', 'f22', 'f5']) g = Tensor(shape=(4, 4, 3, 3), labels=['g41', 'g42', 'g31', 'g32']) makeLink('a3', 'e31', a, e) makeLink('a5', 'b5', a, b) makeLink('a7', 'c72', a, c)
# -*- coding: utf-8 -*- """ Created on Thu Apr 16 20:52:12 2015 @author: Travis """ from numpy import * import CTL ### # Runtime ### #recordList = ["DowJonesOpen", "DowJonesClose", "DowJonesHigh", "DowJonesLow"] recordList = ["DowJonesClose"] order = 10 lowPassFilter = [None, [order, 0.1], [order, 0.25], [order, 0.5]] for lpf in lowPassFilter: CTL.fullAnalysis(recordList, 'financialdata/', lpf)
# -*- coding: utf-8 -*- """ Created on Thu Apr 16 20:52:12 2015 @author: Travis """ from numpy import * import CTL ### # Runtime ### recordList = ["AirTemp"] order = 10 lowPassFilter = [None, [order, 0.1], [order, 0.25], [order, 0.5]] for lpf in lowPassFilter: CTL.fullAnalysis(recordList, 'climatedata/', lpf)
def HOTRGImpurityExample(beta=0.5): print('test magnet for Ising model, beta = {}'.format(beta)) # beta = 0.6 symmetryBroken = 1e-5 a = squareIsingTensor(beta=beta, symmetryBroken=symmetryBroken) hotrg = HOTRG(a, chiH=16) for _ in range(20): hotrg.iterate() mTensor = squareIsingTensor(beta=beta, obs="M", symmetryBroken=symmetryBroken) impurityTN = ImpurityTensorNetwork([a, mTensor], 2) impurityTN.setRG(hotrg) for _ in range(20): impurityTN.iterate() M = impurityTN.measureObservables() M = [x[1] for x in M] exactM = infiniteIsingExactM(1.0 / beta) print('magnet = {}'.format(M[-1] * 0.5)) print('exact magnet = {}'.format(exactM)) if __name__ == '__main__': CTL.setXP(None) CTL.setXP(np) simplestExample() HOTRGImpurityExample(beta=0.6) twoTensorsContraction()