| 56 |
|
Amat = new double[nElements*9]; |
| 57 |
|
mu = new double[nElements]; |
| 58 |
|
ul = new double[nElements*3]; |
| 59 |
+ |
quat = new double[nElements*4]; |
| 60 |
|
|
| 61 |
|
// init directional values to zero |
| 62 |
|
|
| 85 |
|
ul[index3+0] = 1.0; |
| 86 |
|
ul[index3+1] = 0.0; |
| 87 |
|
ul[index3+2] = 0.0; |
| 88 |
+ |
|
| 89 |
+ |
quat[i*4] = 1.0; |
| 90 |
+ |
quat[i*4+1] = 0.0; |
| 91 |
+ |
quat[i*4+2] = 0.0; |
| 92 |
+ |
quat[i*4+3] = 0.0; |
| 93 |
+ |
|
| 94 |
|
} |
| 95 |
|
|
| 96 |
|
arraysAllocated = true; |
| 113 |
|
Amat = NULL; |
| 114 |
|
mu = NULL; |
| 115 |
|
ul = NULL; |
| 116 |
+ |
quat = NULL; |
| 117 |
|
|
| 118 |
|
|
| 119 |
|
arraysAllocated = false; |
| 127 |
|
double** the_trq, |
| 128 |
|
double** the_Amat, |
| 129 |
|
double** the_mu, |
| 130 |
< |
double** the_ul){ |
| 130 |
> |
double** the_ul, |
| 131 |
> |
double** the_quat){ |
| 132 |
|
int index3, index9; |
| 133 |
|
|
| 134 |
|
if( arraysAllocated ){ |
| 143 |
|
*the_Amat = &(Amat[index9]); |
| 144 |
|
*the_mu = &(mu[index]); |
| 145 |
|
*the_ul = &(ul[index3]); |
| 146 |
+ |
*the_quat = &(quat[index*4]); |
| 147 |
|
} |
| 148 |
|
else{ |
| 149 |
|
|