# How to print values on gauss points

(Difference between revisions)
 Revision as of 08:09, 3 May 2012 (view source)Kazem (Talk | contribs) (Created page with "Following steps are necessary to be abale to print values on '''Gauss points''': gid_io.PrintOnGaussPoints(CONTACT_STICK, self.model_part, time);") Revision as of 08:25, 3 May 2012 (view source)Kazem (Talk | contribs) Newer edit → Line 1: Line 1: Following steps are necessary to be abale to print values on '''Gauss points''': Following steps are necessary to be abale to print values on '''Gauss points''': − gid_io.PrintOnGaussPoints(CONTACT_STICK, self.model_part, time); + 1. Call '''PrintOnGaussPoints''' function of the '''gid_io''' with the variable that you want to see, i.e '''CONTACT_STICK''', to + the ''python'' file + + gid_io.PrintOnGaussPoints(CONTACT_STICK, self.model_part, time); + + Note that the variable '''must''' be a double, a vector of doubles or a Matrix of doubles. + + 2. This function in his turn calls a function inside your '''element''' that is called '''GetValueOnIntegrationPoints'''. + + This is an standard function of class of '''element''' and has 3 representatives for double, vector of doubles and Matrix of doubles. + + One example of it can be found as , + + void PfemContactElement3D::GetValueOnIntegrationPoints(const Variable& rVariable, std::vector& rValues, const ProcessInfo& rCurrentProcessInfo) { + if (rVariable == CONTACT_STICK) { + for (unsigned int PointNumber = 0; + PointNumber < 1; PointNumber++) { + rValues[PointNumber] = double(this->GetValue(IS_CONTACT_MASTER)); + } + } + + }

## Revision as of 08:25, 3 May 2012

Following steps are necessary to be abale to print values on Gauss points:

1. Call PrintOnGaussPoints function of the gid_io with the variable that you want to see, i.e CONTACT_STICK, to the python file

``` gid_io.PrintOnGaussPoints(CONTACT_STICK, self.model_part, time);
```

Note that the variable must be a double, a vector of doubles or a Matrix of doubles.

2. This function in his turn calls a function inside your element that is called GetValueOnIntegrationPoints.

This is an standard function of class of element and has 3 representatives for double, vector of doubles and Matrix of doubles.

One example of it can be found as ,

```	void PfemContactElement3D::GetValueOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rValues, const ProcessInfo& rCurrentProcessInfo) {
```

if (rVariable == CONTACT_STICK) { for (unsigned int PointNumber = 0; PointNumber < 1; PointNumber++) { rValues[PointNumber] = double(this->GetValue(IS_CONTACT_MASTER)); } }

}