> but instead of of getting back the full vector, I want to pull out > each component separately > > double phi_i_s_l = fe_values[saturation].gradient (i, > q)[l]; > > by looping over > > for (unsigned int l=0; l<dim;++l) > > I've tried the above, in various guises, and am getting variations of > the following error: > > error: cannot convert ‘dealii::Tensor<1, 2>’ to ‘double’ in > initialization
That really doesn't make any sense. Are you sure 'saturation' is really a *scalar* extractor? If it is, can you send me your code so I can take a look? Cheers W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
