On 10/19/11 5:02 AM, Lorenzo Alessio Botti wrote:
Looking at the implementation of FE::side_map, it seems like it should
be possible to avoid calling inverse_map() in the reinit(elem,side,...), even when
there are user-supplied points, by doing something like:
// Find where the integratio
>
> Looking at the implementation of FE::side_map, it seems like it should
> be possible to avoid calling inverse_map() in the reinit(elem,side,...), even
> when there are user-supplied points, by doing something like:
>
> // Find where the integration points are located on the
> // full eleme
On 10/18/11 3:32 PM, Boyce Griffith wrote:
>
>
> On 10/18/11 1:48 PM, Boyce Griffith wrote:
>>> However, if your points are in the physical frame the best option
>>> (you don't need to modify the reinit behavior) is probably to find the
>>> reference points
>>> in your code with inverse map and t
On 10/18/11 1:48 PM, Boyce Griffith wrote:
>> However, if your points are in the physical frame the best option
>> (you don't need to modify the reinit behavior) is probably to find the
>> reference points
>> in your code with inverse map and then use reinit providing a vector
>> of points
>>
>>
On 10/18/11 3:02 PM, Boyce Griffith wrote:
>
>
> On 10/6/11 3:59 AM, Lorenzo Alessio Botti wrote:
>> However, if your points are in the physical frame the best option
>> (you don't need to modify the reinit behavior) is probably to find the
>> reference points
>> in your code with inverse map and
On 10/6/11 3:59 AM, Lorenzo Alessio Botti wrote:
> However, if your points are in the physical frame the best option
> (you don't need to modify the reinit behavior) is probably to find the
> reference points
> in your code with inverse map and then use reinit providing a vector of points
>
> FE
On 10/6/11 3:59 AM, Lorenzo Alessio Botti wrote:
> Hi Boyce,
>
>>// Find where the integration points are located on the
>>// full element.
>>std::vector qp; //this->inverse_map (elem, xyz, qp, tolerance);
>>this->side_map(elem, side.get(), s, qrule->get_points(), qp);
>>
>> Can
Hi Boyce,
> // Find where the integration points are located on the
> // full element.
> std::vector qp; //this->inverse_map (elem, xyz, qp, tolerance);
> this->side_map(elem, side.get(), s, qrule->get_points(), qp);
>
> Can this revert to using inverse_map if there are user-supplied poin
Is FE::reinit() in fe_boundary.C intended to work when qrule is NULL?
Line 203 of fe_boundary.C calls qrule->get_points() whether or not qrule
has been set to anything:
// Find where the integration points are located on the
// full element.
std::vector qp; //this->inverse_map (elem, xy
