Added new template SIAM/ISO (particle-hole symmetric case)

templates/SIAM/ISO/Collect

@@ -0,0 +1,31 @@
+#!/usr/bin/env perl
+
+open (O, ">info") or die;
+
+$x = `cat ops`;
+print O "ops:$x";
+$x = (-e "specs" ? `cat specs` : "\n");
+print O "specs:$x";
+$x = (-e "specd" ? `cat specd` : "\n");
+print O "specd:$x";
+$x = (-e "spect" ? `cat spect` : "\n");
+print O "spect:$x";
+$x = (-e "specq" ? `cat specq` : "\n");
+print O "specq:$x";
+$x = (-e "specot" ? `cat specot` : "\n");
+print O "specot:$x";
+
+open (F, "<param") or die;
+
+foreach (<F>) {
+    last if /\[param\]/;
+    /(.*)=/ or next;
+    $params .= "$1 ";
+}
+
+chop($params); # drop trailing space
+
+print O "params:$params\n";
+
+$x = (-e "polarized" ? `cat polarized` : "\n");
+print O "polarized:$x";

templates/SIAM/ISO/SIAM.m

@@ -0,0 +1,7 @@
+def1ch[nrimp=1];
+Heps = (-U1 / 2) number[ d[] ];
+Hint = U1 hubbard[ d[] ];
+Himp = Heps + Hint;
+Hhyb = gammaPolCh[1] hop[f[0], d[]];
+H = H0 + Himp + Hhyb;
+isoselfopd = ( Chop @ Expand @ komutator[ (2 U1) pow[ isospinz[ d[] ], 2 ], d[#1, #2] ] )&;

templates/SIAM/ISO/chi_struct

@@ -0,0 +1 @@
+imp 1

templates/SIAM/ISO/data.in

@@ -0,0 +1,103 @@
+# Input file for NRG Ljubljana
+# symtype 	ISO
+# Using sneg version 	2.0.7
+#!9
+# Number of channels, chain sites, subspaces: 
+1	0	4
+# SCALE 	1
+# Energies (GS energy subtracted, multiplied by 1/SCALE):
+3	1
+1
+DIAG 	ham_3.1
+2	2
+2
+DIAG 	ham_2.2
+1	3
+1
+DIAG 	ham_1.3
+1	1
+2
+DIAG 	ham_1.1
+# Irreducible matrix elements for Wilson chains:
+f 0 0
+6
+1	1	2	2
+0.	-1.4142135623730951
+1.4142135623730951	0.
+1	3	2	2
+1.4142135623730951	0.
+2	2	1	1
+0.	-0.7071067811865475
+0.7071067811865475	0.
+2	2	1	3
+1.224744871391589
+0.
+2	2	3	1
+0.
+1.224744871391589
+3	1	2	2
+0.	1.4142135623730951
+# GS energy in absolute units:
+e
+1
+# Irreducible matrix elements for other operators:
+sn_d
+4
+3	1	3	1
+op.n_d_3.1_3.1
+2	2	2	2
+op.n_d_2.2_2.2
+1	3	1	3
+op.n_d_1.3_1.3
+1	1	1	1
+op.n_d_1.1_1.1
+sn_d^2
+4
+3	1	3	1
+op.n_d^2_3.1_3.1
+2	2	2	2
+op.n_d^2_2.2_2.2
+1	3	1	3
+op.n_d^2_1.3_1.3
+1	1	1	1
+op.n_d^2_1.1_1.1
+shop0
+4
+3	1	3	1
+op.hop0_3.1_3.1
+2	2	2	2
+op.hop0_2.2_2.2
+1	3	1	3
+op.hop0_1.3_1.3
+1	1	1	1
+op.hop0_1.1_1.1
+dA_d
+6
+1	1	2	2
+op.A_d_1.1_2.2
+1	3	2	2
+op.A_d_1.3_2.2
+2	2	1	1
+op.A_d_2.2_1.1
+2	2	1	3
+op.A_d_2.2_1.3
+2	2	3	1
+op.A_d_2.2_3.1
+3	1	2	2
+op.A_d_3.1_2.2
+tsigma_d
+0
+disoself_d
+6
+1	1	2	2
+op.isoself_d_1.1_2.2
+1	3	2	2
+op.isoself_d_1.3_2.2
+2	2	1	1
+op.isoself_d_2.2_1.1
+2	2	1	3
+op.isoself_d_2.2_1.3
+2	2	3	1
+op.isoself_d_2.2_3.1
+3	1	2	2
+op.isoself_d_3.1_2.2

templates/SIAM/ISO/discretize

@@ -0,0 +1,4 @@
+#!/bin/sh
+echo dos=Gamma_imp_00.dat >>param
+adapt P
+adapt N

templates/SIAM/ISO/gf_struct

@@ -0,0 +1 @@
+imp 1

templates/SIAM/ISO/ham_1.1

@@ -0,0 +1 @@
+{{0, 2*gammaPolCh[1]}, {2*gammaPolCh[1], -1/2*U1}}

templates/SIAM/ISO/ham_1.3

@@ -0,0 +1 @@
+{{-1/2*U1}}

templates/SIAM/ISO/ham_2.2

@@ -0,0 +1 @@
+{{-1/2*U1 + coefzeta[1, 0], -gammaPolCh[1]}, {-gammaPolCh[1], 0}}

templates/SIAM/ISO/ham_3.1

@@ -0,0 +1 @@
+{{coefzeta[1, 0]}}

templates/SIAM/ISO/info

@@ -0,0 +1,8 @@
+ops:n_d n_d^2 hop0 A_d isoself_d sigma_d
+specs:n_d-n_d
+specd:A_d-A_d isoself_d-A_d A_d-isoself_d isoself_d-isoself_d
+spect:sigma_d-sigma_d
+specq:
+specot:
+params:U1
+polarized: