| 230 | } |
| 231 | |
| 232 | realnum *multilevel_susceptibility::cinternal_notowned_ptr(int inotowned, component c, int cmp, |
| 233 | int n, void *P_internal_data) const { |
| 234 | multilevel_data *d = (multilevel_data *)P_internal_data; |
| 235 | if (!d || !d->P[c][cmp] || inotowned < 0 || inotowned >= T) // never true |
| 236 | return NULL; |
| 237 | return d->P[c][cmp][inotowned] + n; |
| 238 | } |
| 239 | |
| 240 | void multilevel_susceptibility::update_P(realnum *W[NUM_FIELD_COMPONENTS][2], |
| 241 | realnum *W_prev[NUM_FIELD_COMPONENTS][2], realnum dt, |