| 56 | 
  | 
  Amat = new double[nElements*9]; | 
| 57 | 
  | 
  mu   = new double[nElements]; | 
| 58 | 
  | 
  ul   = new double[nElements*3]; | 
| 59 | 
+ | 
  rc = new double[nElements*3]; | 
| 60 | 
+ | 
  massRatio = new double[nElements]; | 
| 61 | 
  | 
   | 
| 62 | 
  | 
  // init directional values to zero | 
| 63 | 
  | 
   | 
| 100 | 
  | 
  if(Amat != NULL) delete[] Amat; | 
| 101 | 
  | 
  if(mu   != NULL) delete[] mu; | 
| 102 | 
  | 
  if(ul   != NULL) delete[] ul;   | 
| 103 | 
+ | 
  if(rc != NULL) delete[] rc; | 
| 104 | 
+ | 
  if(massRatio != NULL) delete[] massRatio; | 
| 105 | 
  | 
 | 
| 106 | 
  | 
  pos  = NULL; | 
| 107 | 
  | 
  vel  = NULL; | 
| 110 | 
  | 
  Amat = NULL; | 
| 111 | 
  | 
  mu   = NULL; | 
| 112 | 
  | 
  ul   = NULL; | 
| 113 | 
+ | 
  rc = NULL; | 
| 114 | 
+ | 
  massRatio = NULL; | 
| 115 | 
  | 
 | 
| 116 | 
  | 
  arraysAllocated = false; | 
| 117 | 
  | 
  nElements = 0; | 
| 124 | 
  | 
                                double** the_trq, | 
| 125 | 
  | 
                                double** the_Amat, | 
| 126 | 
  | 
                                double** the_mu,  | 
| 127 | 
< | 
                                double** the_ul ){ | 
| 127 | 
> | 
                                double** the_ul,  | 
| 128 | 
> | 
                                double** the_rc, | 
| 129 | 
> | 
                                double** the_massRatio){ | 
| 130 | 
  | 
  int index3, index9; | 
| 131 | 
  | 
 | 
| 132 | 
  | 
  if( arraysAllocated ){ | 
| 141 | 
  | 
    *the_Amat = &(Amat[index9]); | 
| 142 | 
  | 
    *the_mu   = &(mu[index]);  | 
| 143 | 
  | 
    *the_ul   = &(ul[index3]); | 
| 144 | 
+ | 
    *the_rc = &(rc[index3]); | 
| 145 | 
+ | 
    *the_massRatio = &(massRatio[index]); | 
| 146 | 
  | 
  } | 
| 147 | 
  | 
  else{ | 
| 148 | 
  | 
 |