38 lines
1 KiB
C
38 lines
1 KiB
C
|
|
/*************************************************************************/
|
|
/*************************************************************************/
|
|
/*************************************************************************/
|
|
|
|
void DEF_DEN(double x_gas[][3], double m_gas[], double h_gas[], int sosed_gas[][NB], double DEN_gas[])
|
|
{
|
|
|
|
double Xij, Yij, Zij, Rij, tmp;
|
|
|
|
/* SPH style */
|
|
|
|
for(i=0;i<N_GAS;i++)
|
|
{
|
|
tmp = 0.0;
|
|
|
|
for(k=0;k<NB;k++)
|
|
{
|
|
j = sosed_gas[i][k];
|
|
|
|
Xij = x_gas[i][0]-x_gas[j][0];
|
|
Yij = x_gas[i][1]-x_gas[j][1];
|
|
Zij = x_gas[i][2]-x_gas[j][2];
|
|
|
|
Rij = sqrt( Xij*Xij + Yij*Yij + Zij*Zij );
|
|
|
|
tmp = tmp + 0.5 * m_gas[j] * ( W(Rij,h_gas[i]) + W(Rij,h_gas[j]) );
|
|
} /* j */
|
|
|
|
DEN_gas[i] = tmp;
|
|
|
|
} /* i */
|
|
|
|
}
|
|
|
|
/*************************************************************************/
|
|
/*************************************************************************/
|
|
/*************************************************************************/
|