sandialabs / compadre

Compadre (Compatible Particle Discretization and Remap)
Other
22 stars 6 forks source link

Change eigensolver initial iterate (bugfix) #291

Closed kuberry closed 2 years ago

kuberry commented 2 years ago

Previously initial iterate was {1,1,1} for eigensolver in determining tangent plane. If the normal happens to be {1,1,1} this causes undesired behavior (becomes leading tangent direction). Now the initial iterate is chosen as a random vector.

Addresses issue #290