File : mm1norm_now.adb

with Reps, Messages, Fouriers_Init, Lin, Lin.IO;
use Reps, Messages, Fouriers_Init, Lin, Lin.IO;
with RG_Ops;

with Intervals, Intervals.Ops, Intervals.IO;
use Intervals, Intervals.Ops, Intervals.IO;
pragma Elaborate_All(Intervals.Ops, Intervals.IO);

procedure MM1Norm_noW is

  subtype Scalar is Interval;

  Lmax: constant Positive := 25;
  Wmax: constant Positive := 25;
  Pdeg: constant Positive := 10;
  Dho:  constant Natural  :=  8;

  Sigma: constant Rep := 0.85001;
  Kappa: constant Rep := Sigma/0.4;

  package MM1Norm_noW_Ops is new RG_Ops (Lmax => Lmax, Wmax => Wmax, Pdeg => Pdeg, Dho => Dho,
                                         Sigma => Sigma, Kappa => Kappa, Scalar => Scalar);
  use MM1Norm_noW_Ops;
  use Scal_Fou, Scal_Fou_Ops, Scal_Fou_IO;

  R: constant Weights := (0.85,0.15);      -- for domain of RG

  Iname: constant String := "25-16-8";     -- name stem for input Hamiltonians
  Yname: constant String := "85B04";       -- name stem for input matrices
  HDec:  constant Boolean := False;        -- use decimal IO for Hamiltonians

  Nyy:   constant Positive := 294911;      -- maximum number of entries in sparse matrix W
  Err:   constant Radius := 1.0E-12;       -- precision for NumFixComp3 etc.

  B3: Radius;
  Hq,Hzp,H1,Hdq,Hdp: Hamilton;
  W: Sparse(0 .. Nyy);

begin
  Trace;                                   -- activate TraceEnter/TraceLeave
  ShowParam;
  Show("  R =",R);
  Show("Err =",Err);

  Read("W-" & Yname & ".tdr",W);
  WCheck(W);

  Read("Fix-" & Iname & ".nh",0,(Two,One),H1,HDec);
  Read("Hdq-" & Iname & ".nh",1,(Two,One),Hdq,HDec);
  Read("Hdp-" & Iname & ".nh",0,(Two,One),Hdp,HDec);
  Read("Hzp-" & Iname & ".nh",0,(Two,One),Hzp,HDec);
  B3 := Sup(Norm(R,Hzp));
  Read("Hq-"  & Iname & ".nh",1,(R.Q,B3),Hq,HDec);

  DRG1Norms(Err,R,Hq,Hzp,H1,Hdq,Hdp,W);

end MM1Norm_noW;