> For the complete documentation index, see [llms.txt](https://applied-physics.gitbook.io/warp-factory/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://applied-physics.gitbook.io/warp-factory/modules/analyzer-module/getenergyconditions.md). # getEnergyConditions ## Description The energy conditions are evaluated point-wise for all spacetime gridpoints of the stress-energy tensor. The function only evaluates the Null Energy Condition (NEC), Weak Energy Condition (WEC), Strong Energy Condition (SEC), and Dominant Energy Condition (DEC).
Null Energy Condition (NEC) The NEC expresses the observed energy density for null (light-ray) observers. For this to be physical, energy density must be non-negative as viewed by any light rays. The NEC is determined by: $$\Xi\_{NEC}(X) = T\_{\hat{\mu}\hat{\nu}}(X) k^{\hat{\mu}} k^{\hat{\nu}} \ge 0$$ where $$T\_{\hat{\mu}\hat{\nu}}$$ is the Eulerian stress-energy tensor and $$k^{\hat{\mu}}$$ is a null vector field.
Weak Energy Condition (WEC) The WEC expresses the observed energy density for timelike (matter) observers. For this to be physical, energy density must be non-negative as viewed by any timelike observer. The WEC is determined by: $$\Xi\_{WEC}(X) = T\_{\hat{\mu}\hat{\nu}}(X) V^{\hat{\mu}} V^{\hat{\nu}} \ge 0$$ where $$T\_{\hat{\mu}\hat{\nu}}$$ is the Eulerian stress-energy tensor and $$V^{\hat{\mu}}$$ is a timelike vector field.
Strong Energy Condition (SEC) The SEC expresses the tidal effect of gravity acting on observers. For this to be physical, the gravitational tidal effect should be non-negative, i.e., matter should always gravitate toward other matter. The SEC is determined by: $$\Xi\_{SEC}(X) = \left(T\_{\hat{\mu}\hat{\nu}}(X) - \frac{1}{2} T(X) \eta\_{\hat{\mu}\hat{\nu}} \right)V^{\hat{\mu}} V^{\hat{\nu}} \ge 0$$ where $$T\_{\hat{\mu}\hat{\nu}}$$ is the Eulerian stress-energy tensor, $$\eta\_{\hat{\mu}\hat{\nu}}$$ is the Minkowski metric, and $$V^{\hat{\mu}}$$ is a timelike vector field.
Dominant Energy Condition (DEC) The DEC expresses the velocity of observed matter flow. For this to be physical, the matter should not be seen flowing faster than light by any timelike observer. $$\Upsilon^{\hat{\mu}}(X) = -T^{\hat{\mu}}\_{\ \ \hat{\nu}}(X) V^{\hat{\nu}}$$ where $$T^{\hat{\mu}}\_{\ \ \hat{\nu}}$$ is the mixed form of the Eulerian stress-energy tensor and $$V^{\hat{\mu}}$$ is a timelike vector field. The energy condition is for this to be future pointing, so this is further evaluated as: $$\xi\_D(X) = \eta\_{\hat{\mu}\hat{\nu}} \Upsilon^{\hat{\mu}}(X) \Upsilon^{\hat{\nu}}(X) \le 0$$ which is the returned energy condition. **The value returned is additionally flipped in sign such that violations are negative values, so as to match the returns of other energy conditions.**
For more general background on the energy conditions, please read: {% embed url="" %} ## Method The evaluation of the point-wise energy conditions samples a set of observer vector fields. The null vector is defined within a Cartesian locally-Minkowskian space. The direction of the spatial velocity is sampled from a set of vectors that map to evenly distributed points on a sphere, which, in a spherical coordinate representation, result in the set of four-velocities that map to different observers. The resulting vectors for each observer are then normalized. For timelike vectors, the four-velocity is defined in the same manner but is additionally scaled from 0 to 1 for a specified number of timelike samples for each of the spatial velocity direction samples. All vectors are then normalized. ## Syntax `[``map``,`` ``vec``,`` ``vectorFieldOut``] = getEnergyConditions(``energyTensor``,`` ``metric``,`` ``condition``,`` ``numAngularVec``,`` ``numTimeVec``,`` ``returnVec``,`` ``tryGPU``)` ### Input Arguments {% hint style="info" %} blue are required inputs. orange are optional inputs with native default values. {% endhint %}
InputsFormatTypeDescription
energyTensor structobjectInput stress-energy tensor to determine energy conditions.
metricstructobjectInput metric which determined the stress-energy tensor.
condition1x1 arraystringThe selected energy condition. Can be either "Null", "Weak", "Strong", or "Dominant".
numAngularVec1x1 arrayintegerThe set number of even-spaced angular samples of four-velocity orientations. This sets the number of directions to sample.
numTimeVec1x1 arrayintegerThe set number of time samples. This evenly reduces the magnitude of the four-velocity spatial part in the given sample number. The default value is 10.
returnVec1x1 arrayintegerReturn the vectors and each of their respective energy condition evaluation results. The default value is 0.
tryGPU1x1 arrayintegerUse the GPU for evaluations, input either 1 for true or 0 for false. The default value is 0.
### Output Arguments
OutputsFormatTypeDescription
map4D arraydouble

Returned energy condition at each spacetime grid point. The array format is specified as:

[x_0, x_1, x_2, x_3]

where x_\muis the spacetime dimension.

vec6D arraydouble

The energy condition evaluation for each of the sampled vectors at each spacetime grid point. The array format is specified as:

[x_0, x_1, x_2, x_3, N_{ang}, N_{time}]

where x_\mu is the spacetime dimensions, N_{ang} is the angular samples, and N_{time} is the time samples.

This is only returned if returnVec is not 0.

vectorFieldOut3D or 4D arraydouble

The set of observer 4-vectors that are used to sample the energy conditions at each point in spacetime. If the energy condition requires a null vector field then the return is a 3D array specified as:

[n_{obs}, N_{ang}]

where n_{obs} is the observer four-velocity (length of 4) and N_{ang} is the length of the angular samples. If the energy condition requires a timelike vector field, then the return is a 4D array specified as:

[n_{obs}, N_{ang}, N_{time}]

where N_{time} is the length of the time samples.

This is only returned if returnVec is not 0.