[Paraview] streamline tracer does not work on cut filter results

[Dominik Szczerba]

Done, id = 0007197
-- Dominik

Utkarsh Ayachit wrote:
> Dominik, can you please report a bug for this issue and assign it to 
> Berk. I will make sure he takes a look at it when he gets back.
> 
> Thanks,
> Utkarsh
> 
> Dominik Szczerba wrote:
>> vtkStreamTracer.cxx:
>>
>> void vtkStreamTracer::SetMaximumPropagation( double max)
>> {
>>   if ( max == this->MaximumPropagation.Interval )
>> ...
>>
>> void vtkStreamTracer::SetMinimumIntegrationStep(double step)
>> {
>>   if ( step == this->MinimumIntegrationStep.Interval )
>> ...
>>
>> void vtkStreamTracer::SetMaximumPropagation( double max)
>> {
>>   if ( max == this->MaximumPropagation.Interval )
>> ...
>>
>> void vtkStreamTracer::SetMinimumIntegrationStep(double step)
>> {
>>   if ( step == this->MinimumIntegrationStep.Interval )
>> ...
>>
>>
>> (!!!)
>>
>> -- Dominik
>>
>> Dominik Szczerba wrote:
>>> OK I traced the problem down to using the Intel compiler for my 
>>> build. The same code compiled with GNU works.
>>>
>>> To the developer of the stream tracer class: are you using things like
>>>
>>> if(x==y)
>>>
>>> for x,y float/double?
>>>
>>> I have already noticed GNU and Intel respond differently, with GNU 
>>> more forgiving. The proper (= correct, robust, meaningful) comparison 
>>> should be
>>>
>>> if(fabs(x-y)<eps)
>>>
>>> I do not know if that is the problem, but it is one of the 
>>> differences between the compilers that I noticed.
>>>
>>> -- Dominik
>>>
>>> Dominik Szczerba wrote:
>>>> Ha it even works with official 3.3.0, only not my own compilation 
>>>> from CVS 05 June. So either it is my own compilation, or something 
>>>> still got broken in CVS slightly later. Attached is my CMakeCache. 
>>>> Any ideas?
>>>> - Dominik
>>>>
>>>>
>>>> Dominik Szczerba wrote:
>>>>> Wooo it works with PV 3.2.1 official. So something got broken in 
>>>>> the CVS.
>>>>> -- Dominik
>>>>>
>>>>> Dominik Szczerba wrote:
>>>>>> I already did it, as mentioned in a previous post. No, it does not 
>>>>>> work for me - does it work on my dataset for you?
>>>>>>
>>>>>> In any case, that is a trivialized example, my real case is not 
>>>>>> symmetric and a cut plane is not easily aligned with the axes.
>>>>>>
>>>>>> I suspect this might have to do with:
>>>>>>
>>>>>> 1) the funny range like 1.52e-18 - 1.22e-17 for point1/point2 in 
>>>>>> line source sub-window that I get hitting 'X/Y/Z Axis' button.
>>>>>>
>>>>>> 2) my mesh being quadratic and streamlining not working in such 
>>>>>> cases.
>>>>>>
>>>>>> Thanks for any hints.
>>>>>> -- Dominik
>>>>>>
>>>>>> Michael Ahlmann wrote:
>>>>>>> I had a similar problem when trying to visualize an axi-symmetric 
>>>>>>> data set.  The problem was that the streamlines were not confined 
>>>>>>> to the the axi-symmetric plane, and thus they were cut short at 
>>>>>>> the point where the flow went in the circumferential direction.  
>>>>>>> The fix in my case was to use the calculator to develop a new, 
>>>>>>> simplified data set.
>>>>>>>
>>>>>>> V = V_x * ihat + V_y * jhat + 0 * khat
>>>>>>>
>>>>>>> This way, the axi-symmetric component is zero, and doesn't affect 
>>>>>>> the in plane streamlines.  I wonder if you could do something 
>>>>>>> similar in your case?  Although I presume it would be more 
>>>>>>> complex given that your cutting plane may not lay perfectly 
>>>>>>> within one of the Cartesian planes...
>>>>>>>
>>>>>>> -Michael
>>>>>>> _______________________________________________
>>>>>>> ParaView mailing list
>>>>>>> ParaView at paraview.org
>>>>>>> http://www.paraview.org/mailman/listinfo/paraview
>>>>>>
>>>>>
>>>>
>>>>
>>>> ------------------------------------------------------------------------ 
>>>>
>>>>
>>>> _______________________________________________
>>>> ParaView mailing list
>>>> ParaView at paraview.org
>>>> http://www.paraview.org/mailman/listinfo/paraview
>>>
>>

-- 
Dominik Szczerba, Ph.D.
Biomedical Simulation Group
Computer Vision Lab CH-8092 Zurich
http://www.vision.ee.ethz.ch/~domi