FHGetSelfgl

       FHGetSelfgl  -  calculate  the  renormalized  Higgs self-energies in the
       gaugeless limit

       integer error
       double complex k2
       integer key, dkey, ren
       double complex sig(13), dsig(13)

       integer h0h0, HHHH, A0A0, HmHp
       integer h0HH, h0A0, HHA0
       integer G0G0, h0G0, HHG0, A0G0
       integer GmGp, HmGp
       parameter (h0h0 = 1, HHHH = 2, A0A0 = 3, HmHp = 4)
       parameter (h0HH = 5, h0A0 = 6, HHA0 = 7)
       parameter (G0G0 = 8, h0G0 = 9, HHG0 = 10, A0G0 = 11)
       parameter (GmGp = 12, HmGp = 13)
       #define Key(se) 2**(se-1)

       subroutine FHGetSelfgl(error, k2, key, sig, dkey, dsig, ren)

       FHGetSelfgl[k2, key, dkey, ren]

       FHGetSelfgl returns the renormalized Higgs self-energies in  the  gauge-
       less  limit  at momentum-squared k2.  The flags and parameters must have
       been set before with FHSetFlags(1) and FHSetPara(1).

       error (output)
              zero if successful, otherwise the line number in GetSelfgl.F from
              which the error message was emitted

       k2 (input)
              the momentum-squared k^2 at which the  self-energies  are  evalu-
              ated.

       key (input)
              a flag determining which of the self-energies are actually evalu-
              ated,  e.g.  to  evaluate the h0-h0 self-energy, add Key(h0h0) to
              key.

       dkey (input)
              a flag determining which of the derivatives of the  self-energies
              are  actually  evaluated,  e.g. to evaluate the derivative of the
              h0-h0 self-energy, add Key(h0h0) to dkey.

       ren (input)
              whether the unrenormalized (0) or renormalized (1)  self-energies
              are output.

       sig(h0h0), sig(HHHH), sig(A0A0), sig(HmHp) (output)
              the h0, HH, A0, and Hp self-energies at k^2 = k2.

       sig(h0HH), sig(h0A0), sig(HHA0) (output)
              the h0-HH, h0-A0, and HH-A0 mixing self-energies at k^2 = k2.

       sig(G0G0), sig(h0G0), sig(HHG0), sig(A0G0) (output)
              the neutral Goldstone self-energies at k^2 = k2.

       sig(GmGp), sig(HmGp) (output)
              the charged Goldstone self-energies at k^2 = k2.

       dsig(i) (output)
              the derivatives of the self-energies with respect to k^2 at k^2 =
              k2, where the index i runs as for the sig(i).

       libFH(1)

                                  18-Jul-2018                    FHGETSELFGL(1)