546 |
G=U.Gravitational_Constant |
G=U.Gravitational_Constant |
547 |
m_psi_ref=0. |
m_psi_ref=0. |
548 |
x=domain.getX() |
x=domain.getX() |
549 |
for i in range(self.DIM): |
DIM=domain.getDim() |
550 |
m_psi_ref=m_psi_ref + whereZero(x[i]-inf(x[i])) \ |
m_psi_ref=whereZero(x[DIM-1]-sup(x[DIM-1])) # + whereZero(x[DIM-1]-inf(x[DIM-1])) |
|
+ whereZero(x[i]-sup(x[i])) |
|
|
|
|
551 |
if self.getDataType()==DataSource.GRAVITY: |
if self.getDataType()==DataSource.GRAVITY: |
552 |
rho_ref=0. |
rho_ref=0. |
553 |
for f in self._features: |
for f in self._features: |