-- Make AWS output into m2 input

fnbase = "BR7valuesByTerm";
fn = openOut "BR7valuesByTerm.m2";
fn << "BR7valuesByTerm={" << endl;
for i from 0 to 807 do (
    stri = toString(i);
    if #stri==1 then stri = concatenate("00",toString(i));
    if #stri==2 then stri = concatenate("0",toString(i)); 
    fni = concatenate(fnbase,".",stri,".out");
    s=get fni;
    fn << s;
);
fn << "};" << endl;
close fn;

-- Delete the comma after item 807
-- Then:


-- Load v1, v2, and the AWS output

needsPackage("LieTypes");
load "../LieAlgebraRepresentations.v2.3.m2";
load "../SpinRepresentations.v1.7.m2";


n=5;
k=7;
lambda1={5,0,0,0,0};


EWB = evenWedgeBasis(n);
snk = sort subsets(apply(2^(n-1), i -> i),k);
K=QQ;
ringSubscripts = apply(snk, s -> apply(s, i -> EWB_i));
R=K[apply(ringSubscripts, i -> yy_i),MonomialOrder=>Lex];


load "../MyBasisWords.m2";
load "../Sym2 V1 highest weight vector for 5w1.m2";
load "../Sym2 V2 highest weight vector for 5w1.m2";
T = join(terms(v1),terms(v2));
T = apply(T, t -> leadMonomial(t));
T = sort unique T;
load "BR7valuesByTerm.m2";



v1BR7= last(sum apply(#T, i -> coefficient(T_i,v1)*(BR7valuesByTerm_i)));
v2BR7= last(sum apply(#T, i -> coefficient(T_i,v2)*(BR7valuesByTerm_i)));

fn = openOut "v1andv2OnBR7.m2";
fn << toString(v1BR7) << endl;
fn << toString(v2BR7) << endl;
close fn;