OpenCores
URL https://opencores.org/ocsvn/openrisc/openrisc/trunk

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-old/] [gcc-4.2.2/] [gcc/] [testsuite/] [gcc.dg/] [ultrasp6.c] - Blame information for rev 149

Go to most recent revision | Details | Compare with Previous | View Log

Line No. Rev Author Line
1 149 jeremybenn
/* PR target/7784 */
2
/* Originator: Peter van Hoof <p.van-hoof@qub.ac.uk> */
3
/* { dg-do compile { target sparc-*-* } } */
4
/* { dg-options "-O2 -mcpu=ultrasparc" } */
5
 
6
typedef struct
7
{
8
  float EnergyErg;
9
  float ots;
10
} EmLine;
11
 
12
extern const int ipH_LIKE ;
13
extern const int ipHYDROGEN ;
14
extern const int ipH1s;
15
extern const int ipH2s;
16
extern const int ipH2p;
17
 
18
extern EmLine ****EmisLines;
19
 
20
typedef struct
21
{
22
  long n;
23
  long s;
24
  long l;
25
} Elevels;
26
 
27
extern struct t_iso
28
{
29
  float ***Pop2Ion;
30
  long int numLevels[2][30L];
31
} iso;
32
 
33
extern struct t_LineSave
34
{
35
  long int nsum;
36
  long int ndsum;
37
  long int nComment;
38
  long int npxdd;
39
  long int ipass;
40
  char chHoldComments[10][200];
41
} LineSave;
42
 
43
extern struct t_hydro
44
{
45
  int lgHydEmiss;
46
  float **pestrk ;
47
} hydro;
48
 
49
extern struct t_dense
50
{
51
  double DensityLaw[10];
52
  float frad[500];
53
  float fhden[500];
54
  float den0;
55
  double eden;
56
} dense;
57
 
58
extern struct t_abund
59
{
60
  float xIonFracs[30L +3][30L +1];
61
} abund;
62
 
63
extern struct t_CaseBHS
64
{
65
  long int nDensity[2][8] , ntemp[2][8] , ncut[2][8] ;
66
  int lgHCaseBOK[2][8];
67
} CaseBHS ;
68
 
69
extern struct t_smbeta
70
{
71
  float SimHBeta,
72
    cn4861,
73
    cn1216,
74
    sv4861,
75
    sv1216;
76
} smbeta;
77
 
78
extern struct t_phycon
79
{
80
  float te;
81
} phycon;
82
 
83
 
84
extern struct t_sphere
85
{
86
  int lgSphere;
87
  float covgeo;
88
} sphere;
89
 
90
void linadd(double xInten, float wavelength, char *chLab, char chInfo);
91
 
92
extern struct t_radiusVar
93
{
94
  int lgDrNeg;
95
  double dVeff;
96
} radius;
97
 
98
void lines_hydro(void)
99
{
100
  long int i, nelem, ipHi, ipLo;
101
  double hbetab, em , EmisFac, pump;
102
  char chLabel[5];
103
 
104
  linadd(abund.xIonFracs[ipHYDROGEN][1]*iso.Pop2Ion[ipH_LIKE][ipHYDROGEN][3]*hydro.pestrk[3][2]*3.025e-12, 6563,"Strk",'i');
105
 
106
  linadd(abund.xIonFracs[ipHYDROGEN][1]*iso.Pop2Ion[ipH_LIKE][ipHYDROGEN][4]*hydro.pestrk[4][2]*4.084e-12, 4861,"Strk",'i');
107
 
108
  linadd(abund.xIonFracs[ipHYDROGEN][1]*iso.Pop2Ion[ipH_LIKE][ipHYDROGEN][4]*hydro.pestrk[4][3]*1.059e-12, 18751,"Strk",'i');
109
 
110
  linadd(abund.xIonFracs[ipHYDROGEN][1]*iso.Pop2Ion[ipH_LIKE][ipHYDROGEN][5]*hydro.pestrk[5][4]*4.900e-13, 40512,"Strk",'i');
111
 
112
  ((void)((LineSave.ipass <1 || EmisLines[ipH_LIKE][ipHYDROGEN][ipH2p][ipH1s].ots>= 0.) || (__assert("LineSave.ipass <1 || EmisLines[ipH_LIKE][ipHYDROGEN][ipH2p][ipH1s].ots>= 0.", "lines_hydro.c", 118), 0)));
113
 
114
  linadd(EmisLines[ipH_LIKE][ipHYDROGEN][3][ipH2s].ots*EmisLines[ipH_LIKE][ipHYDROGEN][3][ipH2s].EnergyErg, 6563,"Dest",'i');
115
 
116
  linadd(EmisLines[ipH_LIKE][ipHYDROGEN][5][4].ots*EmisLines[ipH_LIKE][ipHYDROGEN][5][4].EnergyErg,40516, "Dest",'i');
117
 
118
  smbeta.SimHBeta = smbeta.SimHBeta/(float)radius.dVeff*sphere.covgeo;
119
 
120
  linadd(smbeta.SimHBeta,4861,"Q(H)",'i');
121
 
122
  smbeta.SimHBeta = smbeta.SimHBeta*(float)radius.dVeff/sphere.covgeo;
123
 
124
  for( nelem=0; nelem < 30L; nelem++ )
125
  {
126
    int iCase;
127
    for( iCase=0; iCase<2; ++iCase )
128
    {
129
      char chAB[2]={'A','B'};
130
      char chLab[5]="Ca  ";
131
 
132
      for( ipLo=1+iCase; ipLo<(((6)<(iso.numLevels[ipH_LIKE][nelem])) ? (6) : (5)); ++ipLo )
133
      {
134
        for( ipHi=ipLo+1; ipHi< (((ipLo+5)<(iso.numLevels[ipH_LIKE][nelem])) ? (ipLo+5) : (iso.numLevels[ipH_LIKE][nelem])); ++ipHi )
135
        {
136
          float wl;
137
 
138
          hbetab = HSRate( ipHi,ipLo , nelem+1, phycon.te , dense.eden, chAB[iCase] );
139
          if( hbetab<=0. )
140
            CaseBHS.lgHCaseBOK[iCase][nelem] = 0;
141
 
142
          if( !hydro.lgHydEmiss )
143
            hbetab *= abund.xIonFracs[nelem][nelem+1]*dense.eden;
144
 
145
         linadd(hbetab,wl,chLab,'i' );
146
        }
147
      }
148
    }
149
  }
150
}

powered by: WebSVN 2.1.0

© copyright 1999-2024 OpenCores.org, equivalent to Oliscience, all rights reserved. OpenCores®, registered trademark.