Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
From: Sergei Isakov isakov@itp.phys.ethz.ch Date: October 28, 2010 18:59:30 GMT+02:00 To: Matthias Troyer troyer@phys.ethz.ch Subject: Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
From: Urs Gerber ITP gerberu@itp.unibe.ch Date: October 28, 2010 15:43:38 GMT+02:00 To: comp-phys-alps-users@lists.phys.ethz.ch Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 Reply-To: comp-phys-alps-users@lists.phys.ethz.ch
Here are the files I use and the output I get:
bilayer-test:
LATTICE_LIBRARY="./../lattices-ug.xml"; LATTICE="latticegraph square bilayer"; MODEL_LIBRARY="./../models.xml"; MODEL="spin"; local_S=1/2; L = 12; Jxy0=1; Jz0=3.3; Jxy1=3.45; Jz1=3.45; THERMALIZATION=1000; SWEEPS=1000; beta=48; {h0=7.5;}
bilayer-test.in.xml:
<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd"> <OUTPUT file="bilayer-test.out.xml"/> <TASK status="new"> <INPUT file="bilayer-test.task1.in.xml"/> <OUTPUT file="bilayer-test.task1.out.xml"/> </TASK> </JOB>
output in the terminal:
Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov isakov@itp.phys.ethz.ch
using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer troyer@comp-phys.org. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).
based on the ALPS libraries version 2.0.0b4 available from http://alps.comp-phys.org/ copyright (c) 1994-2010 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).
parsing task files ... Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.001 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Created run 1 locally Starting task 1.
After this it stops creating output in the terminal and no output files are produced.
Hello Matthias Here is the lattices-ug.xml file. With directed loop algorithm v. 1.1 and ALPS libraries version 1.3.3 the same input files as below work perfectly fine. Urs
Matthias Troyer wrote:
Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
*From:* Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch> *Date:* October 28, 2010 18:59:30 GMT+02:00 *To:* Matthias Troyer <troyer@phys.ethz.ch mailto:troyer@phys.ethz.ch> *Subject:* *Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3*
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
From: Urs Gerber ITP <gerberu@itp.unibe.ch mailto:gerberu@itp.unibe.ch> Date: October 28, 2010 15:43:38 GMT+02:00 To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 Reply-To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch
Here are the files I use and the output I get:
bilayer-test:
LATTICE_LIBRARY="./../lattices-ug.xml"; LATTICE="latticegraph square bilayer"; MODEL_LIBRARY="./../models.xml"; MODEL="spin"; local_S=1/2; L = 12; Jxy0=1; Jz0=3.3; Jxy1=3.45; Jz1=3.45; THERMALIZATION=1000; SWEEPS=1000; beta=48; {h0=7.5;}
bilayer-test.in.xml:
<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd">
<OUTPUT file="bilayer-test.out.xml"/> <TASK status="new"> <INPUT file="bilayer-test.task1.in.xml"/> <OUTPUT file="bilayer-test.task1.out.xml"/> </TASK> </JOB>
output in the terminal:
Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch>
using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer@comp-phys.org mailto:troyer@comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).
based on the ALPS libraries version 2.0.0b4 available from http://alps.comp-phys.org/ copyright (c) 1994-2010 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).
parsing task files ... Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.001 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Created run 1 locally Starting task 1.
After this it stops creating output in the terminal and no output files are produced.
Dear Urs,
I've tested your input file and It seems to work fine with my installation. Which version of lp_solve do you use ?
Best Fabien
Le 1 nov. 10 à 15:46, Urs Gerber ITP a écrit :
Hello Matthias Here is the lattices-ug.xml file. With directed loop algorithm v. 1.1 and ALPS libraries version 1.3.3 the same input files as below work perfectly fine. Urs
Matthias Troyer wrote:
Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
*From:* Sergei Isakov <isakov@itp.phys.ethz.ch <mailto:isakov@itp.phys.ethz.ch
*Date:* October 28, 2010 18:59:30 GMT+02:00 *To:* Matthias Troyer <troyer@phys.ethz.ch <mailto:troyer@phys.ethz.ch
*Subject:* *Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3*
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
From: Urs Gerber ITP <gerberu@itp.unibe.ch <mailto:gerberu@itp.unibe.ch
>
Date: October 28, 2010 15:43:38 GMT+02:00 To: comp-phys-alps-users@lists.phys.ethz.ch <mailto:comp-phys-alps-users@lists.phys.ethz.ch
Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 Reply-To: comp-phys-alps-users@lists.phys.ethz.ch <mailto:comp-phys-alps-users@lists.phys.ethz.ch
Here are the files I use and the output I get:
bilayer-test:
LATTICE_LIBRARY="./../lattices-ug.xml"; LATTICE="latticegraph square bilayer"; MODEL_LIBRARY="./../models.xml"; MODEL="spin"; local_S=1/2; L = 12; Jxy0=1; Jz0=3.3; Jxy1=3.45; Jz1=3.45; THERMALIZATION=1000; SWEEPS=1000; beta=48; {h0=7.5;}
bilayer-test.in.xml:
<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd ">
<OUTPUT file="bilayer-test.out.xml"/> <TASK status="new"> <INPUT file="bilayer-test.task1.in.xml"/> <OUTPUT file="bilayer-test.task1.out.xml"/> </TASK> </JOB>
output in the terminal:
Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch>
using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer@comp-phys.org mailto:troyer@comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).
based on the ALPS libraries version 2.0.0b4 available from http://alps.comp-phys.org/ copyright (c) 1994-2010 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).
parsing task files ... Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.001 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Created run 1 locally Starting task 1.
After this it stops creating output in the terminal and no output files are produced.
<LATTICES>
<LATTICE name="square bilayer" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="0d lattice" dimension="0"/>
<LATTICE name="chain lattice" dimension="1"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="square lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/a</ VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="a"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</ VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="a/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</ VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="4"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</ VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</ VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey cubic lattice" dimension="3"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="triangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>a/2 a*sqrt(3)/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/sqrt(3)</VECTOR><VECTOR>0 4*pi/a/sqrt(3)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="centered rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></ BASIS> <RECIPROCALBASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</ VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="oblique lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="phi" default="Pi/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>b*sin(phi) b*cos(phi)</VECTOR></ BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/tan(phi)</VECTOR><VECTOR>0 2*pi/b/cos(phi)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="simple cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 a</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="face-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a/2 a/2 0 </VECTOR> <VECTOR>a/2 0 a/2</VECTOR> <VECTOR>0 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 2*pi/a -2*pi/a </VECTOR> <VECTOR>2*pi/a -2*pi/a 2*pi/a </VECTOR> <VECTOR>-2*pi/a 2*pi/a 2*pi/a </VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="body-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR> a/2 a/2 -a/2</VECTOR> <VECTOR> a/2 -a/2 a/2</VECTOR> <VECTOR>-a/2 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR> 2*pi/a 2*pi/a 0</VECTOR> <VECTOR> 2*pi/a 0 2*pi/a</VECTOR> <VECTOR> 0 2*pi/a 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="tetragonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="orthorhombic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="hexagonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>a/2 a*sqrt(3)/2 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a -2*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 4*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<UNITCELL name="isolated" dimension="1"> <VERTEX/> </UNITCELL>
<UNITCELL name="simple1d" dimension="1"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/ ></EDGE> </UNITCELL>
<UNITCELL name="simple2d" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> </UNITCELL>
<UNITCELL name="triangular" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="-1 1"/></EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d" dimension="2"> <VERTEX/> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_inner_anisotropy_layers" dimension="3"> <VERTEX type="0"><COORDINATE>0 0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="2" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="3" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="1" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0 0"/> <TARGET vertex="2" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="3" offset="0 0 0"/> <TARGET vertex="3" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="4" offset="0 0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_staggered" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_isotropy" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder_1" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.75</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.75</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="6"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="7"/> <TARGET vertex="8"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="5" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="nnn2d" dimension="2"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 1"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 -1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d" dimension="3" vertices="1"> <VERTEX/> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d_1" dimension="3" vertices="1"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="2band1d" dimension="1"> <VERTEX type="0"/> <VERTEX type="1"/> <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="2"><SOURCE vertex="2" offset="0"/><TARGET vertex="2" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="Kagome" dimension="2"> <VERTEX><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0 0.5</COORDINATE></VERTEX>
<EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="2" offset="-1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3" offset="1 -1"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3" offset="0 -1"/></EDGE>
</UNITCELL>
<UNITCELL name="coupled ladders" dimension="2"> <VERTEX type="0"/> <VERTEX type="0"/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0"/ ></EDGE> <EDGE type="0"><SOURCE vertex="2"/><TARGET vertex="2" offset="1 0"/ ></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2" offset="0 1"/ ></EDGE> </UNITCELL>
<UNITCELL name="unitcell square bilayer" dimension="3"> <VERTEX id="1" type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 0 1"/> </EDGE> </UNITCELL>
<LATTICEGRAPH name = "square lattice 3x3"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="3"/> <EXTENT dimension="2" size="3"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice 4x4"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="4"/> <EXTENT dimension="2" size="4"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "dimer"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="2"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "site"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="1"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="isolated"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice 1"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "coupled ladders"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="coupled ladders"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "triangular lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="triangular"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "frustrated square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="nnn2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="ladder rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey staggered anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_staggered"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey isotropic lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_isotropy"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey inner anisotropic layers lattice"> <FINITELATTICE> <LATTICE ref=" honey cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="Z" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="Z"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_inner_anisotropy_layers"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder_1 anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" ladder honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="open"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "depleted square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <DEPLETION> <VERTEX probability="DEPLETION" seed="DEPLETION_SEED"/> </DEPLETION> </LATTICEGRAPH>
<LATTICEGRAPH name = "Kagome lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="Kagome"/> </LATTICEGRAPH>
<LATTICEGRAPH name="latticegraph square bilayer" dimension="3"> <FINITELATTICE> <LATTICE ref="square bilayer"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="periodic"/> <BOUNDARY dimension="3" type="open"/> </FINITELATTICE> <UNITCELL ref="unitcell square bilayer"/> </LATTICEGRAPH>
<GRAPH name="5-site dimerized" vertices="5"> <EDGE type="0" source="1" target="2"/> <EDGE type="1" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="1" source="4" target="5"/> </GRAPH>
<GRAPH name="4-site mixed" vertices="4"> <VERTEX id="1" type="0"/> <VERTEX id="2" type="1"/> <VERTEX id="3" type="0"/> <VERTEX id="4" type="1"/> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="0" source="4" target="1"/> <EDGE type="1" source="1" target="3"/> <EDGE type="1" source="2" target="4"/> </GRAPH>
<GRAPH name="triangle" vertices="3"> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="1"/> </GRAPH>
</LATTICES>
Urs, can you try with lpsolve version 4 instead of 5.1?
On 1 Nov 2010, at 15:46, Urs Gerber ITP wrote:
Hello Matthias Here is the lattices-ug.xml file. With directed loop algorithm v. 1.1 and ALPS libraries version 1.3.3 the same input files as below work perfectly fine. Urs
Matthias Troyer wrote:
Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
*From:* Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch> *Date:* October 28, 2010 18:59:30 GMT+02:00 *To:* Matthias Troyer <troyer@phys.ethz.ch mailto:troyer@phys.ethz.ch> *Subject:* *Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3*
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
From: Urs Gerber ITP <gerberu@itp.unibe.ch mailto:gerberu@itp.unibe.ch> Date: October 28, 2010 15:43:38 GMT+02:00 To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 Reply-To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch
Here are the files I use and the output I get:
bilayer-test:
LATTICE_LIBRARY="./../lattices-ug.xml"; LATTICE="latticegraph square bilayer"; MODEL_LIBRARY="./../models.xml"; MODEL="spin"; local_S=1/2; L = 12; Jxy0=1; Jz0=3.3; Jxy1=3.45; Jz1=3.45; THERMALIZATION=1000; SWEEPS=1000; beta=48; {h0=7.5;}
bilayer-test.in.xml:
<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd"> <OUTPUT file="bilayer-test.out.xml"/> <TASK status="new"> <INPUT file="bilayer-test.task1.in.xml"/> <OUTPUT file="bilayer-test.task1.out.xml"/> </TASK> </JOB>
output in the terminal:
Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch>
using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer@comp-phys.org mailto:troyer@comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).
based on the ALPS libraries version 2.0.0b4 available from http://alps.comp-phys.org/ copyright (c) 1994-2010 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).
parsing task files ... Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.001 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Created run 1 locally Starting task 1.
After this it stops creating output in the terminal and no output files are produced.
<LATTICES>
<LATTICE name="square bilayer" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="0d lattice" dimension="0"/>
<LATTICE name="chain lattice" dimension="1"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="square lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="a"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="a/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="4"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey cubic lattice" dimension="3"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="triangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>a/2 a*sqrt(3)/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/sqrt(3)</VECTOR><VECTOR>0 4*pi/a/sqrt(3)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="centered rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="oblique lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="phi" default="Pi/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>b*sin(phi) b*cos(phi)</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/tan(phi)</VECTOR><VECTOR>0 2*pi/b/cos(phi)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="simple cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 a</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="face-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a/2 a/2 0 </VECTOR> <VECTOR>a/2 0 a/2</VECTOR> <VECTOR>0 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 2*pi/a -2*pi/a </VECTOR> <VECTOR>2*pi/a -2*pi/a 2*pi/a </VECTOR> <VECTOR>-2*pi/a 2*pi/a 2*pi/a </VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="body-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR> a/2 a/2 -a/2</VECTOR> <VECTOR> a/2 -a/2 a/2</VECTOR> <VECTOR>-a/2 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR> 2*pi/a 2*pi/a 0</VECTOR> <VECTOR> 2*pi/a 0 2*pi/a</VECTOR> <VECTOR> 0 2*pi/a 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="tetragonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="orthorhombic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="hexagonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>a/2 a*sqrt(3)/2 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a -2*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 4*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<UNITCELL name="isolated" dimension="1"> <VERTEX/> </UNITCELL>
<UNITCELL name="simple1d" dimension="1"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="simple2d" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> </UNITCELL>
<UNITCELL name="triangular" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="-1 1"/></EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d" dimension="2"> <VERTEX/> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_inner_anisotropy_layers" dimension="3"> <VERTEX type="0"><COORDINATE>0 0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="2" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="3" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="1" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0 0"/> <TARGET vertex="2" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="3" offset="0 0 0"/> <TARGET vertex="3" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="4" offset="0 0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_staggered" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_isotropy" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder_1" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.75</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.75</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="6"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="7"/> <TARGET vertex="8"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="5" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="nnn2d" dimension="2"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 1"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 -1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d" dimension="3" vertices="1"> <VERTEX/> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d_1" dimension="3" vertices="1"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="2band1d" dimension="1"> <VERTEX type="0"/> <VERTEX type="1"/> <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="2"><SOURCE vertex="2" offset="0"/><TARGET vertex="2" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="Kagome" dimension="2"> <VERTEX><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0 0.5</COORDINATE></VERTEX>
<EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="2" offset="-1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3" offset="1 -1"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3" offset="0 -1"/></EDGE>
</UNITCELL>
<UNITCELL name="coupled ladders" dimension="2"> <VERTEX type="0"/> <VERTEX type="0"/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="0"><SOURCE vertex="2"/><TARGET vertex="2" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2" offset="0 1"/></EDGE> </UNITCELL>
<UNITCELL name="unitcell square bilayer" dimension="3"> <VERTEX id="1" type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 0 1"/> </EDGE> </UNITCELL>
<LATTICEGRAPH name = "square lattice 3x3"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="3"/> <EXTENT dimension="2" size="3"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice 4x4"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="4"/> <EXTENT dimension="2" size="4"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "dimer"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="2"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "site"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="1"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="isolated"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice 1"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "coupled ladders"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="coupled ladders"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "triangular lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="triangular"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "frustrated square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="nnn2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="ladder rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey staggered anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_staggered"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey isotropic lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_isotropy"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey inner anisotropic layers lattice"> <FINITELATTICE> <LATTICE ref=" honey cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="Z" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="Z"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_inner_anisotropy_layers"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder_1 anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" ladder honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="open"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "depleted square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <DEPLETION> <VERTEX probability="DEPLETION" seed="DEPLETION_SEED"/> </DEPLETION> </LATTICEGRAPH>
<LATTICEGRAPH name = "Kagome lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="Kagome"/> </LATTICEGRAPH>
<LATTICEGRAPH name="latticegraph square bilayer" dimension="3"> <FINITELATTICE> <LATTICE ref="square bilayer"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="periodic"/> <BOUNDARY dimension="3" type="open"/> </FINITELATTICE> <UNITCELL ref="unitcell square bilayer"/> </LATTICEGRAPH>
<GRAPH name="5-site dimerized" vertices="5"> <EDGE type="0" source="1" target="2"/> <EDGE type="1" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="1" source="4" target="5"/> </GRAPH>
<GRAPH name="4-site mixed" vertices="4"> <VERTEX id="1" type="0"/> <VERTEX id="2" type="1"/> <VERTEX id="3" type="0"/> <VERTEX id="4" type="1"/> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="0" source="4" target="1"/> <EDGE type="1" source="1" target="3"/> <EDGE type="1" source="2" target="4"/> </GRAPH>
<GRAPH name="triangle" vertices="3"> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="1"/> </GRAPH>
</LATTICES>
I use lp_solve_5.1. Actually I'm running ALPS on a supercomputer. During the CECAM workshop Ryo Igarashi and Bela Bauer made the installation for me. When I remember correctly there was a problem with installing lpsolve version 4 and therefore we took 5.1 . I will try to reinstall version 4 by myself and tell you what the result of this test is.
Fabien, which installation did you use to test ?
Urs
Matthias Troyer wrote:
Urs, can you try with lpsolve version 4 instead of 5.1?
On 1 Nov 2010, at 15:46, Urs Gerber ITP wrote:
Hello Matthias Here is the lattices-ug.xml file. With directed loop algorithm v. 1.1 and ALPS libraries version 1.3.3 the same input files as below work perfectly fine. Urs
Matthias Troyer wrote:
Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
*From:* Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch> *Date:* October 28, 2010 18:59:30 GMT+02:00 *To:* Matthias Troyer <troyer@phys.ethz.ch mailto:troyer@phys.ethz.ch> *Subject:* *Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3*
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
From: Urs Gerber ITP <gerberu@itp.unibe.ch mailto:gerberu@itp.unibe.ch> Date: October 28, 2010 15:43:38 GMT+02:00 To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 Reply-To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch
Here are the files I use and the output I get:
bilayer-test:
LATTICE_LIBRARY="./../lattices-ug.xml"; LATTICE="latticegraph square bilayer"; MODEL_LIBRARY="./../models.xml"; MODEL="spin"; local_S=1/2; L = 12; Jxy0=1; Jz0=3.3; Jxy1=3.45; Jz1=3.45; THERMALIZATION=1000; SWEEPS=1000; beta=48; {h0=7.5;}
bilayer-test.in.xml:
<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd"> <OUTPUT file="bilayer-test.out.xml"/> <TASK status="new"> <INPUT file="bilayer-test.task1.in.xml"/> <OUTPUT file="bilayer-test.task1.out.xml"/> </TASK> </JOB>
output in the terminal:
Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch>
using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer@comp-phys.org mailto:troyer@comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).
based on the ALPS libraries version 2.0.0b4 available from http://alps.comp-phys.org/ copyright (c) 1994-2010 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).
parsing task files ... Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.001 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Created run 1 locally Starting task 1.
After this it stops creating output in the terminal and no output files are produced.
<LATTICES>
<LATTICE name="square bilayer" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="0d lattice" dimension="0"/>
<LATTICE name="chain lattice" dimension="1"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="square lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="a"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="a/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="4"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey cubic lattice" dimension="3"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="triangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>a/2 a*sqrt(3)/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/sqrt(3)</VECTOR><VECTOR>0 4*pi/a/sqrt(3)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="centered rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="oblique lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="phi" default="Pi/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>b*sin(phi) b*cos(phi)</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/tan(phi)</VECTOR><VECTOR>0 2*pi/b/cos(phi)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="simple cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 a</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="face-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a/2 a/2 0 </VECTOR> <VECTOR>a/2 0 a/2</VECTOR> <VECTOR>0 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 2*pi/a -2*pi/a </VECTOR> <VECTOR>2*pi/a -2*pi/a 2*pi/a </VECTOR> <VECTOR>-2*pi/a 2*pi/a 2*pi/a </VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="body-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR> a/2 a/2 -a/2</VECTOR> <VECTOR> a/2 -a/2 a/2</VECTOR> <VECTOR>-a/2 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR> 2*pi/a 2*pi/a 0</VECTOR> <VECTOR> 2*pi/a 0 2*pi/a</VECTOR> <VECTOR> 0 2*pi/a 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="tetragonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="orthorhombic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="hexagonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>a/2 a*sqrt(3)/2 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a -2*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 4*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<UNITCELL name="isolated" dimension="1"> <VERTEX/> </UNITCELL>
<UNITCELL name="simple1d" dimension="1"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="simple2d" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> </UNITCELL>
<UNITCELL name="triangular" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="-1 1"/></EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d" dimension="2"> <VERTEX/> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_inner_anisotropy_layers" dimension="3"> <VERTEX type="0"><COORDINATE>0 0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="2" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="3" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="1" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0 0"/> <TARGET vertex="2" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="3" offset="0 0 0"/> <TARGET vertex="3" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="4" offset="0 0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_staggered" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_isotropy" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder_1" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.75</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.75</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="6"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="7"/> <TARGET vertex="8"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="5" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="nnn2d" dimension="2"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 1"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 -1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d" dimension="3" vertices="1"> <VERTEX/> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d_1" dimension="3" vertices="1"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="2band1d" dimension="1"> <VERTEX type="0"/> <VERTEX type="1"/> <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="2"><SOURCE vertex="2" offset="0"/><TARGET vertex="2" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="Kagome" dimension="2"> <VERTEX><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0 0.5</COORDINATE></VERTEX>
<EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="2" offset="-1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3" offset="1 -1"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3" offset="0 -1"/></EDGE>
</UNITCELL>
<UNITCELL name="coupled ladders" dimension="2"> <VERTEX type="0"/> <VERTEX type="0"/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="0"><SOURCE vertex="2"/><TARGET vertex="2" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2" offset="0 1"/></EDGE> </UNITCELL>
<UNITCELL name="unitcell square bilayer" dimension="3"> <VERTEX id="1" type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 0 1"/> </EDGE> </UNITCELL>
<LATTICEGRAPH name = "square lattice 3x3"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="3"/> <EXTENT dimension="2" size="3"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice 4x4"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="4"/> <EXTENT dimension="2" size="4"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "dimer"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="2"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "site"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="1"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="isolated"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice 1"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "coupled ladders"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="coupled ladders"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "triangular lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="triangular"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "frustrated square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="nnn2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="ladder rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey staggered anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_staggered"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey isotropic lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_isotropy"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey inner anisotropic layers lattice"> <FINITELATTICE> <LATTICE ref=" honey cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="Z" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="Z"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_inner_anisotropy_layers"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder_1 anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" ladder honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="open"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "depleted square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <DEPLETION> <VERTEX probability="DEPLETION" seed="DEPLETION_SEED"/> </DEPLETION> </LATTICEGRAPH>
<LATTICEGRAPH name = "Kagome lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="Kagome"/> </LATTICEGRAPH>
<LATTICEGRAPH name="latticegraph square bilayer" dimension="3"> <FINITELATTICE> <LATTICE ref="square bilayer"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="periodic"/> <BOUNDARY dimension="3" type="open"/> </FINITELATTICE> <UNITCELL ref="unitcell square bilayer"/> </LATTICEGRAPH>
<GRAPH name="5-site dimerized" vertices="5"> <EDGE type="0" source="1" target="2"/> <EDGE type="1" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="1" source="4" target="5"/> </GRAPH>
<GRAPH name="4-site mixed" vertices="4"> <VERTEX id="1" type="0"/> <VERTEX id="2" type="1"/> <VERTEX id="3" type="0"/> <VERTEX id="4" type="1"/> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="0" source="4" target="1"/> <EDGE type="1" source="1" target="3"/> <EDGE type="1" source="2" target="4"/> </GRAPH>
<GRAPH name="triangle" vertices="3"> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="1"/> </GRAPH>
</LATTICES>
I use lpsolve 4, on my Mac, and found no problem with your input file.
Fabien
On 11/02/2010 04:39 PM, Urs Gerber ITP wrote:
I use lp_solve_5.1. Actually I'm running ALPS on a supercomputer. During the CECAM workshop Ryo Igarashi and Bela Bauer made the installation for me. When I remember correctly there was a problem with installing lpsolve version 4 and therefore we took 5.1 . I will try to reinstall version 4 by myself and tell you what the result of this test is.
Fabien, which installation did you use to test ?
Urs
Urs,
can you try with the latest nightly snapshot?
Matthias
On Nov 2, 2010, at 4:39 PM, Urs Gerber ITP wrote:
I use lp_solve_5.1. Actually I'm running ALPS on a supercomputer. During the CECAM workshop Ryo Igarashi and Bela Bauer made the installation for me. When I remember correctly there was a problem with installing lpsolve version 4 and therefore we took 5.1 . I will try to reinstall version 4 by myself and tell you what the result of this test is.
Fabien, which installation did you use to test ?
Urs
Matthias Troyer wrote:
Urs, can you try with lpsolve version 4 instead of 5.1?
On 1 Nov 2010, at 15:46, Urs Gerber ITP wrote:
Hello Matthias Here is the lattices-ug.xml file. With directed loop algorithm v. 1.1 and ALPS libraries version 1.3.3 the same input files as below work perfectly fine. Urs
Matthias Troyer wrote:
Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
*From:* Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch> *Date:* October 28, 2010 18:59:30 GMT+02:00 *To:* Matthias Troyer <troyer@phys.ethz.ch mailto:troyer@phys.ethz.ch> *Subject:* *Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3*
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
> From: Urs Gerber ITP <gerberu@itp.unibe.ch mailto:gerberu@itp.unibe.ch> > Date: October 28, 2010 15:43:38 GMT+02:00 > To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch > Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 > Reply-To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch > > Here are the files I use and the output I get: > > bilayer-test: > > LATTICE_LIBRARY="./../lattices-ug.xml"; > LATTICE="latticegraph square bilayer"; > MODEL_LIBRARY="./../models.xml"; > MODEL="spin"; > local_S=1/2; > L = 12; > Jxy0=1; > Jz0=3.3; > Jxy1=3.45; > Jz1=3.45; > THERMALIZATION=1000; > SWEEPS=1000; > beta=48; > {h0=7.5;} > > > bilayer-test.in.xml: > > <?xml version="1.0" encoding="UTF-8"?> > <?xml-stylesheet type="text/xsl" href="ALPS.xsl"?> > <JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd"> > <OUTPUT file="bilayer-test.out.xml"/> > <TASK status="new"> > <INPUT file="bilayer-test.task1.in.xml"/> > <OUTPUT file="bilayer-test.task1.out.xml"/> > </TASK> > </JOB> > > output in the terminal: > > Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 > available from http://alps.comp-phys.org/ > copyright (c) 2003-2010 by Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch> > > using the ALPS parallelizing scheduler > copyright (c) 1994-2006 by Matthias Troyer <troyer@comp-phys.org mailto:troyer@comp-phys.org>. > see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998). > > based on the ALPS libraries version 2.0.0b4 > available from http://alps.comp-phys.org/ > copyright (c) 1994-2010 by the ALPS collaboration. > Consult the web page for license details. > For details see the publication: > A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007). > > parsing task files ... > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0.84127 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0.84127 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 5 > > Final solution 0.84127 at iteration 5. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 5, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 5.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 5 > > Final solution 0.84127 at iteration 5. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 5, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 5.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0.465116 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0.465116 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 5 > > Final solution 0.465116 at iteration 5. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 5, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 5.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 42 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 5 > > Final solution 0.465116 at iteration 5. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 5, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 5.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.001 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Model name: '' - run #1 > Objective: Minimize(R0) > > SUBMITTED > Model size: 18 constraints, 9 variables, 48 non-zeros. > Constraints: 18 equality, 0 GUB, 0 SOS. > Variables: 0 integer, 0 semi-cont., 0 SOS. > > Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2. > > Optimal solution with dual simplex at iteration 6 > > Final solution 0 at iteration 6. > > Excellent numeric accuracy ||*|| = 0 > > Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. > In the total iteration count 6, 0 (0.0%) were minor/bound swaps. > There were 0 refactorizations, 0 triggered by time and 0 by density. > ... on average 6.0 major pivots per refactorization. > Total solver time was 0.000 seconds. > Created run 1 locally > Starting task 1. > > > After this it stops creating output in the terminal and no output files are produced. > >
<LATTICES>
<LATTICE name="square bilayer" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="0d lattice" dimension="0"/>
<LATTICE name="chain lattice" dimension="1"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="square lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="a"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="a/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="4"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey cubic lattice" dimension="3"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="triangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>a/2 a*sqrt(3)/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/sqrt(3)</VECTOR><VECTOR>0 4*pi/a/sqrt(3)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="centered rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="oblique lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="phi" default="Pi/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>b*sin(phi) b*cos(phi)</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/tan(phi)</VECTOR><VECTOR>0 2*pi/b/cos(phi)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="simple cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 a</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="face-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a/2 a/2 0 </VECTOR> <VECTOR>a/2 0 a/2</VECTOR> <VECTOR>0 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 2*pi/a -2*pi/a </VECTOR> <VECTOR>2*pi/a -2*pi/a 2*pi/a </VECTOR> <VECTOR>-2*pi/a 2*pi/a 2*pi/a </VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="body-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR> a/2 a/2 -a/2</VECTOR> <VECTOR> a/2 -a/2 a/2</VECTOR> <VECTOR>-a/2 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR> 2*pi/a 2*pi/a 0</VECTOR> <VECTOR> 2*pi/a 0 2*pi/a</VECTOR> <VECTOR> 0 2*pi/a 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="tetragonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="orthorhombic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="hexagonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>a/2 a*sqrt(3)/2 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a -2*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 4*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<UNITCELL name="isolated" dimension="1"> <VERTEX/> </UNITCELL>
<UNITCELL name="simple1d" dimension="1"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="simple2d" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> </UNITCELL>
<UNITCELL name="triangular" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="-1 1"/></EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d" dimension="2"> <VERTEX/> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_inner_anisotropy_layers" dimension="3"> <VERTEX type="0"><COORDINATE>0 0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="2" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="3" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="1" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0 0"/> <TARGET vertex="2" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="3" offset="0 0 0"/> <TARGET vertex="3" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="4" offset="0 0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_staggered" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_isotropy" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder_1" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.75</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.75</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="6"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="7"/> <TARGET vertex="8"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="5" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="nnn2d" dimension="2"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 1"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 -1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d" dimension="3" vertices="1"> <VERTEX/> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d_1" dimension="3" vertices="1"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="2band1d" dimension="1"> <VERTEX type="0"/> <VERTEX type="1"/> <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="2"><SOURCE vertex="2" offset="0"/><TARGET vertex="2" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="Kagome" dimension="2"> <VERTEX><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0 0.5</COORDINATE></VERTEX>
<EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="2" offset="-1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3" offset="1 -1"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3" offset="0 -1"/></EDGE>
</UNITCELL>
<UNITCELL name="coupled ladders" dimension="2"> <VERTEX type="0"/> <VERTEX type="0"/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="0"><SOURCE vertex="2"/><TARGET vertex="2" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2" offset="0 1"/></EDGE> </UNITCELL>
<UNITCELL name="unitcell square bilayer" dimension="3"> <VERTEX id="1" type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 0 1"/> </EDGE> </UNITCELL>
<LATTICEGRAPH name = "square lattice 3x3"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="3"/> <EXTENT dimension="2" size="3"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice 4x4"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="4"/> <EXTENT dimension="2" size="4"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "dimer"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="2"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "site"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="1"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="isolated"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice 1"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "coupled ladders"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="coupled ladders"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "triangular lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="triangular"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "frustrated square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="nnn2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="ladder rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey staggered anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_staggered"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey isotropic lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_isotropy"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey inner anisotropic layers lattice"> <FINITELATTICE> <LATTICE ref=" honey cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="Z" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="Z"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_inner_anisotropy_layers"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder_1 anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" ladder honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="open"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "depleted square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <DEPLETION> <VERTEX probability="DEPLETION" seed="DEPLETION_SEED"/> </DEPLETION> </LATTICEGRAPH>
<LATTICEGRAPH name = "Kagome lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="Kagome"/> </LATTICEGRAPH>
<LATTICEGRAPH name="latticegraph square bilayer" dimension="3"> <FINITELATTICE> <LATTICE ref="square bilayer"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="periodic"/> <BOUNDARY dimension="3" type="open"/> </FINITELATTICE> <UNITCELL ref="unitcell square bilayer"/> </LATTICEGRAPH>
<GRAPH name="5-site dimerized" vertices="5"> <EDGE type="0" source="1" target="2"/> <EDGE type="1" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="1" source="4" target="5"/> </GRAPH>
<GRAPH name="4-site mixed" vertices="4"> <VERTEX id="1" type="0"/> <VERTEX id="2" type="1"/> <VERTEX id="3" type="0"/> <VERTEX id="4" type="1"/> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="0" source="4" target="1"/> <EDGE type="1" source="1" target="3"/> <EDGE type="1" source="2" target="4"/> </GRAPH>
<GRAPH name="triangle" vertices="3"> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="1"/> </GRAPH>
</LATTICES>
Dear Urs,
Nobody can reproduce your problems? Could you please try again with the latest nightly snapshot?
Matthias
On Nov 1, 2010, at 3:46 PM, Urs Gerber ITP wrote:
Hello Matthias Here is the lattices-ug.xml file. With directed loop algorithm v. 1.1 and ALPS libraries version 1.3.3 the same input files as below work perfectly fine. Urs
Matthias Troyer wrote:
Urs, can you please send your input files, especially the lattice file. Could you please also check whether the dirloop_sse_v1 code, which is the old code, works on your input?
Matthias
Begin forwarded message:
*From:* Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch> *Date:* October 28, 2010 18:59:30 GMT+02:00 *To:* Matthias Troyer <troyer@phys.ethz.ch mailto:troyer@phys.ethz.ch> *Subject:* *Re: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3*
I don't really understand this output. I probably need that lattices-ug.xml file to check. Does dirloop_sse_v1 work?
Best, Sergei
On 28 Oct 2010, at 18:47, Matthias Troyer wrote:
Begin forwarded message:
From: Urs Gerber ITP <gerberu@itp.unibe.ch mailto:gerberu@itp.unibe.ch> Date: October 28, 2010 15:43:38 GMT+02:00 To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch Subject: [ALPS-users] Problem in ALPS 2.0b4 which does not occur in ALPS 1.3.3 Reply-To: comp-phys-alps-users@lists.phys.ethz.ch mailto:comp-phys-alps-users@lists.phys.ethz.ch
Here are the files I use and the output I get:
bilayer-test:
LATTICE_LIBRARY="./../lattices-ug.xml"; LATTICE="latticegraph square bilayer"; MODEL_LIBRARY="./../models.xml"; MODEL="spin"; local_S=1/2; L = 12; Jxy0=1; Jz0=3.3; Jxy1=3.45; Jz1=3.45; THERMALIZATION=1000; SWEEPS=1000; beta=48; {h0=7.5;}
bilayer-test.in.xml:
<?xml version="1.0" encoding="UTF-8"?>
<?xml-stylesheet type="text/xsl" href="ALPS.xsl"?>
<JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd"> <OUTPUT file="bilayer-test.out.xml"/> <TASK status="new"> <INPUT file="bilayer-test.task1.in.xml"/> <OUTPUT file="bilayer-test.task1.out.xml"/> </TASK> </JOB>
output in the terminal:
Quantum Monte Carlo simulations using the SSE algorithm v. 4.0 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov <isakov@itp.phys.ethz.ch mailto:isakov@itp.phys.ethz.ch>
using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer@comp-phys.org mailto:troyer@comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).
based on the ALPS libraries version 2.0.0b4 available from http://alps.comp-phys.org/ copyright (c) 1994-2010 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).
parsing task files ... Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.84127 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.84127 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0.465116 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 42 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 5
Final solution 0.465116 at iteration 5.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 5, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 5.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.001 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Model name: '' - run #1 Objective: Minimize(R0)
SUBMITTED Model size: 18 constraints, 9 variables, 48 non-zeros. Constraints: 18 equality, 0 GUB, 0 SOS. Variables: 0 integer, 0 semi-cont., 0 SOS.
Using DUAL simplex for phase 1 and PRIMAL simplex for phase 2.
Optimal solution with dual simplex at iteration 6
Final solution 0 at iteration 6.
Excellent numeric accuracy ||*|| = 0
Memo: Largest [etaPFI v1.0] inv(B) had 0 NZ entries, 0.0x largest basis. In the total iteration count 6, 0 (0.0%) were minor/bound swaps. There were 0 refactorizations, 0 triggered by time and 0 by density. ... on average 6.0 major pivots per refactorization. Total solver time was 0.000 seconds. Created run 1 locally Starting task 1.
After this it stops creating output in the terminal and no output files are produced.
<LATTICES>
<LATTICE name="square bilayer" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="0d lattice" dimension="0"/>
<LATTICE name="chain lattice" dimension="1"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="square lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="a"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="a/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="ladder honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="4"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey rectangular lattice" dimension="2"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 b</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/b</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="honey cubic lattice" dimension="3"> <PARAMETER name="a" default="2"/> <PARAMETER name="b" default="2"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="triangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>a/2 a*sqrt(3)/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/sqrt(3)</VECTOR><VECTOR>0 4*pi/a/sqrt(3)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="centered rectangular lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>a 0</VECTOR><VECTOR>sqrt(b*b-a*a/4) a/2</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="oblique lattice" dimension="2"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="phi" default="Pi/2"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>b*sin(phi) b*cos(phi)</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a -2*pi/a/tan(phi)</VECTOR><VECTOR>0 2*pi/b/cos(phi)</VECTOR></RECIPROCALBASIS> </LATTICE>
<LATTICE name="simple cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 a</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="face-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR>a/2 a/2 0 </VECTOR> <VECTOR>a/2 0 a/2</VECTOR> <VECTOR>0 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 2*pi/a -2*pi/a </VECTOR> <VECTOR>2*pi/a -2*pi/a 2*pi/a </VECTOR> <VECTOR>-2*pi/a 2*pi/a 2*pi/a </VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="body-centered cubic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <BASIS> <VECTOR> a/2 a/2 -a/2</VECTOR> <VECTOR> a/2 -a/2 a/2</VECTOR> <VECTOR>-a/2 a/2 a/2</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR> 2*pi/a 2*pi/a 0</VECTOR> <VECTOR> 2*pi/a 0 2*pi/a</VECTOR> <VECTOR> 0 2*pi/a 2*pi/a</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="tetragonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 a 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/a 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="orthorhombic lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="b" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>0 b 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a 0 0</VECTOR> <VECTOR>0 2*pi/b 0</VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<LATTICE name="hexagonal lattice" dimension="3"> <PARAMETER name="a" default="1"/> <PARAMETER name="c" default="1"/> <BASIS> <VECTOR>a 0 0</VECTOR> <VECTOR>a/2 a*sqrt(3)/2 0</VECTOR> <VECTOR>0 0 c</VECTOR> </BASIS> <RECIPROCALBASIS> <VECTOR>2*pi/a -2*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 4*pi/a/sqrt(3) 0 </VECTOR> <VECTOR>0 0 2*pi/c</VECTOR> </RECIPROCALBASIS> </LATTICE>
<UNITCELL name="isolated" dimension="1"> <VERTEX/> </UNITCELL>
<UNITCELL name="simple1d" dimension="1"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="simple2d" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> </UNITCELL>
<UNITCELL name="triangular" dimension="2"> <VERTEX/> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="-1 1"/></EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d" dimension="2"> <VERTEX/> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="anisotropic2d_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="1" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_inner_anisotropy_layers" dimension="3"> <VERTEX type="0"><COORDINATE>0 0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="2" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="3" offset="0 1 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/> <TARGET vertex="1" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2" offset="0 0 0"/> <TARGET vertex="2" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="3" offset="0 0 0"/> <TARGET vertex="3" offset="0 0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0 0"/> <TARGET vertex="4" offset="0 0 1"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_staggered" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_isotropy" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="1"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="honeycomb_ladder_1" dimension="2"> <VERTEX type="0"><COORDINATE>0 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.25</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.75</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.5</COORDINATE></VERTEX> <VERTEX type="0"><COORDINATE>0.5 0.75</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="0"> <SOURCE vertex="3"/> <TARGET vertex="4"/> </EDGE> <EDGE type="0"> <SOURCE vertex="4"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5"/> <TARGET vertex="6"/> </EDGE> <EDGE type="0"> <SOURCE vertex="6"/> <TARGET vertex="7"/> </EDGE> <EDGE type="1"> <SOURCE vertex="7"/> <TARGET vertex="8"/> </EDGE> <EDGE type="1"> <SOURCE vertex="5" offset="0 0"/> <TARGET vertex="2" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="3" offset="0 1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="4" offset="0 0"/> <TARGET vertex="1" offset="1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="8" offset="0 0"/> <TARGET vertex="5" offset="1 0"/> </EDGE> </UNITCELL>
<UNITCELL name="nnn2d" dimension="2"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 1"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 -1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d" dimension="3" vertices="1"> <VERTEX/> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="simple3d_1" dimension="3" vertices="1"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0 0"/></EDGE> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="1" offset="0 0 1"/></EDGE> </UNITCELL>
<UNITCELL name="2band1d" dimension="1"> <VERTEX type="0"/> <VERTEX type="1"/> <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="2"><SOURCE vertex="2" offset="0"/><TARGET vertex="2" offset="1"/></EDGE> </UNITCELL>
<UNITCELL name="Kagome" dimension="2"> <VERTEX><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0.5 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>0 0.5</COORDINATE></VERTEX>
<EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="2" offset="-1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="3" offset="1 -1"/></EDGE> <EDGE><SOURCE vertex="1"/><TARGET vertex="3" offset="0 -1"/></EDGE>
</UNITCELL>
<UNITCELL name="coupled ladders" dimension="2"> <VERTEX type="0"/> <VERTEX type="0"/> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="0"><SOURCE vertex="2"/><TARGET vertex="2" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2" offset="0 1"/></EDGE> </UNITCELL>
<UNITCELL name="unitcell square bilayer" dimension="3"> <VERTEX id="1" type="0"><COORDINATE>0 0 0</COORDINATE></VERTEX> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="1 0 0"/> </EDGE> <EDGE type="0"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 1 0"/> </EDGE> <EDGE type="1"> <SOURCE vertex="1" offset="0 0 0"/><TARGET vertex="1" offset="0 0 1"/> </EDGE> </UNITCELL>
<LATTICEGRAPH name = "square lattice 3x3"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="3"/> <EXTENT dimension="2" size="3"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice 4x4"> <FINITELATTICE> <LATTICE ref="square lattice"/> <EXTENT dimension="1" size="4"/> <EXTENT dimension="2" size="4"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "dimer"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="2"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "site"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size="1"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="isolated"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "simple cubic lattice 1"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple3d_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "coupled ladders"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="coupled ladders"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "triangular lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="triangular"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "frustrated square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="nnn2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "2 band open chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <EXTENT dimension="1" size ="L"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="2band1d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref="ladder rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey staggered anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_staggered"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey isotropic lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_isotropy"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey inner anisotropic layers lattice"> <FINITELATTICE> <LATTICE ref=" honey cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="Z" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="Z"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_inner_anisotropy_layers"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "honey ladder_1 anisotropic square lattice"> <FINITELATTICE> <LATTICE ref=" ladder honey rectangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb_ladder_1"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "open ladder"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="open"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "inhomogeneous simple cubic lattice"> <FINITELATTICE> <LATTICE ref="simple cubic lattice"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="W"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY type="open"/> </FINITELATTICE> <UNITCELL ref="simple3d"/> <INHOMOGENEOUS><VERTEX/></INHOMOGENEOUS> </LATTICEGRAPH>
<LATTICEGRAPH name = "depleted square lattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="simple2d"/> <DEPLETION> <VERTEX probability="DEPLETION" seed="DEPLETION_SEED"/> </DEPLETION> </LATTICEGRAPH>
<LATTICEGRAPH name = "Kagome lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="Kagome"/> </LATTICEGRAPH>
<LATTICEGRAPH name="latticegraph square bilayer" dimension="3"> <FINITELATTICE> <LATTICE ref="square bilayer"/> <PARAMETER name="W" default="L"/> <PARAMETER name="H" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <EXTENT dimension="3" size="H"/> <BOUNDARY dimension="1" type="periodic"/> <BOUNDARY dimension="2" type="periodic"/> <BOUNDARY dimension="3" type="open"/> </FINITELATTICE> <UNITCELL ref="unitcell square bilayer"/> </LATTICEGRAPH>
<GRAPH name="5-site dimerized" vertices="5"> <EDGE type="0" source="1" target="2"/> <EDGE type="1" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="1" source="4" target="5"/> </GRAPH>
<GRAPH name="4-site mixed" vertices="4"> <VERTEX id="1" type="0"/> <VERTEX id="2" type="1"/> <VERTEX id="3" type="0"/> <VERTEX id="4" type="1"/> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="4"/> <EDGE type="0" source="4" target="1"/> <EDGE type="1" source="1" target="3"/> <EDGE type="1" source="2" target="4"/> </GRAPH>
<GRAPH name="triangle" vertices="3"> <EDGE type="0" source="1" target="2"/> <EDGE type="0" source="2" target="3"/> <EDGE type="0" source="3" target="1"/> </GRAPH>
</LATTICES>
comp-phys-alps-users@lists.phys.ethz.ch