/*! \file valuenode_bline.cpp
** \brief Template File
**
-** $Id: valuenode_bline.cpp,v 1.1.1.1 2005/01/04 01:23:15 darco Exp $
+** $Id$
**
** \legal
** Copyright (c) 2002-2005 Robert B. Quattlebaum Jr., Adrian Bentley
inline Vector
radial_interpolation(const Vector& a, const Vector& b, float c)
{
+ // if either extreme is zero then use linear interpolation instead
+ if (a.is_equal_to(Vector::zero()) || b.is_equal_to(Vector::zero()))
+ return linear_interpolation(a, b, c);
+
affine_combo<Real,float> mag_combo;
affine_combo<Angle,float> ang_combo;
-
+
Real mag(mag_combo(a.mag(),b.mag(),c));
Angle ang(ang_combo(Angle::tan(a[1],a[0]),Angle::tan(b[1],b[0]),c));
-
+
return Point( mag*Angle::cos(ang).get(),mag*Angle::sin(ang).get() );
}
synfig::convert_bline_to_segment_list(const ValueBase& bline)
{
std::vector<Segment> ret;
-
+
// std::vector<BLinePoint> list(bline.operator std::vector<BLinePoint>());
//std::vector<BLinePoint> list(bline);
std::vector<BLinePoint> list(bline.get_list().begin(),bline.get_list().end());
std::vector<BLinePoint>::const_iterator iter;
BLinePoint prev,first;
-
+
//start with prev = first and iter on the second...
-
+
if(list.empty()) return ValueBase(ret,bline.get_loop());
first = prev = list.front();
-
+
for(iter=++list.begin();iter!=list.end();++iter)
{
ret.push_back(
)
);
prev=*iter;
- }
+ }
if(bline.get_loop())
{
ret.push_back(
//std::vector<BLinePoint> list(bline);
std::vector<BLinePoint> list(bline.get_list().begin(),bline.get_list().end());
std::vector<BLinePoint>::const_iterator iter;
-
+
if(bline.empty())
return ValueBase(ValueBase::TYPE_LIST);
-
+
for(iter=list.begin();iter!=list.end();++iter)
ret.push_back(iter->get_width());
{
if(value.get_type()!=ValueBase::TYPE_LIST)
return 0;
-
+
ValueNode_BLine* value_node(new ValueNode_BLine());
if(!value.empty())
std::vector<BLinePoint>::const_iterator iter;
for(iter=bline_points.begin();iter!=bline_points.end();iter++)
- {
+ {
value_node->add(ValueNode::Handle(ValueNode_Composite::create(*iter)));
}
value_node->set_loop(value.get_loop());
// into a list of BLinePoints. We make an assumption
// that the segment list is continuous(sp), but not necessarily
// smooth.
-
+
value_node->set_loop(false);
// std::vector<Segment> segments(value.operator std::vector<Segment>());
// std::vector<Segment> segments(value);
std::vector<Segment>::const_iterator iter,last(segments.end());
--last;
ValueNode_Const::Handle prev,first;
-
+
for(iter=segments.begin();iter!=segments.end();iter++)
- {
+ {
#define PREV_POINT prev->get_value().get(BLinePoint())
#define FIRST_POINT first->get_value().get(BLinePoint())
#define CURR_POINT curr->get_value().get(BLinePoint())
prev=ValueNode_Const::create(ValueBase::TYPE_BLINEPOINT);
{
BLinePoint prev_point(PREV_POINT);
- prev_point.set_vertex(iter->p1);
+ prev_point.set_vertex(iter->p1);
prev_point.set_tangent1(iter->t1);
prev_point.set_width(0.01);
prev_point.set_origin(0.5);
}
first=prev;
value_node->add(ValueNode::Handle(prev));
-
+
}
if(iter==last && iter->p2.is_equal_to(FIRST_POINT.get_vertex()))
{
}
continue;
}
-
+
ValueNode_Const::Handle curr;
curr=ValueNode_Const::create(ValueBase::TYPE_BLINEPOINT);
{
BLinePoint curr_point(CURR_POINT);
- curr_point.set_vertex(iter->p2);
+ curr_point.set_vertex(iter->p2);
curr_point.set_tangent1(iter->t2);
curr_point.set_width(0.01);
curr_point.set_origin(0.5);
value_node->add(ValueNode::Handle(curr));
prev=curr;
}
-
+
}
break;
default:
break;
}
}
-
-
+
+
return value_node;
}
{
ValueNode_BLine::ListEntry ret;
-
+
synfig::BLinePoint prev,next;
int prev_i,next_i;
else
next_i=index;
prev_i=find_prev_valid_entry(index,time);
-
+
synfig::info("index=%d, next_i=%d, prev_i=%d",index,next_i,prev_i);
-
+
next=(*list[next_i].value_node)(time);
prev=(*list[prev_i].value_node)(time);
-
+
etl::hermite<Vector> curve(prev.get_vertex(),next.get_vertex(),prev.get_tangent2(),next.get_tangent1());
etl::derivative< etl::hermite<Vector> > deriv(curve);
bline_point.set_tangent2(bline_point.get_tangent1());
bline_point.set_split_tangent_flag(false);
bline_point.set_origin(origin);
-
+
ret.value_node=ValueNode_Composite::create(bline_point);
-
+
return ret;
}
ValueNode_BLine::operator()(Time t)const
{
std::vector<BLinePoint> ret_list;
-
+
std::vector<ListEntry>::const_iterator iter,first_iter;
bool first_flag(true);
bool rising;
int index(0);
float next_scale(1.0f);
-
+
BLinePoint prev,first;
first.set_origin(100.0f);
-
+
for(iter=list.begin();iter!=list.end();++iter,index++)
{
float amount(iter->amount_at_time(t,&rising));
-
+
assert(amount>=0.0f);
assert(amount<=1.0f);
-
+
if(amount==1.0f)
{
if(first_flag)
ret_list.push_back(first);
continue;
}
-
+
BLinePoint curr;
curr=(*iter->value_node)(t).get(prev);
-
+
if(next_scale!=1.0f)
{
ret_list.back().set_split_tangent_flag(true);
else
{
ret_list.push_back(curr);
-
}
-
+
prev=curr;
}
else
if(amount>0.0f)
{
std::vector<ListEntry>::const_iterator begin_iter,end_iter;
-
+
// This is where the interesting stuff happens
// We need to seek forward in the list to see what the next
// active point is
-
+
BLinePoint curr;
BLinePoint begin; // begin of dynamic group
- BLinePoint end; // End of dynamic group
- Time blend_time;
+ BLinePoint end; // end of dynamic group
int dist_from_begin(0), dist_from_end(0);
BLinePoint ret;
-
- Time begin_time;
- Time end_time;
-
+
+ Time off_time, on_time;
+
if(!rising)
{
- try{ end_time=iter->find_prev(t)->get_time(); }
- catch(...) { end_time=Time::begin(); }
- try{ begin_time=iter->find_next(t)->get_time(); }
- catch(...) { begin_time=Time::end(); }
+ try{ on_time=iter->find_prev(t)->get_time(); }
+ catch(...) { on_time=Time::begin(); }
+ try{ off_time=iter->find_next(t)->get_time(); }
+ catch(...) { off_time=Time::end(); }
}
else
{
- try{ begin_time=iter->find_prev(t)->get_time(); }
- catch(...) { begin_time=Time::begin(); }
- try{ end_time=iter->find_next(t)->get_time(); }
- catch(...) { end_time=Time::end(); }
+ try{ off_time=iter->find_prev(t)->get_time(); }
+ catch(...) { off_time=Time::begin(); }
+ try{ on_time=iter->find_next(t)->get_time(); }
+ catch(...) { on_time=Time::end(); }
}
- blend_time=begin_time;
- curr=(*iter->value_node)(end_time).get(curr);
+ curr=(*iter->value_node)(on_time).get(curr);
// curr=(*iter->value_node)(t).get(curr);
-
+
// Find "end" of dynamic group
end_iter=iter;
// for(++end_iter;begin_iter!=list.end();++end_iter)
for(++end_iter;end_iter!=list.end();++end_iter)
if(end_iter->amount_at_time(t)>amount)
{
- end=(*end_iter->value_node)(blend_time).get(prev);
+ end=(*end_iter->value_node)(off_time).get(prev);
break;
}
-
+
// If we did not find an end of the dynamic group...
if(end_iter==list.end())
{
if(get_loop())
{
end_iter=first_iter;
- end=(*end_iter->value_node)(blend_time).get(prev);
+ end=(*end_iter->value_node)(off_time).get(prev);
// end=first;
}
else
{
// Writeme!
end_iter=first_iter;
- end=(*end_iter->value_node)(blend_time).get(prev);
+ end=(*end_iter->value_node)(off_time).get(prev);
// end=first;
}
}
else
break;
}
-
+
--begin_iter;
dist_from_begin++;
-
+
if(begin_iter==iter)
break;
if(begin_iter->amount_at_time(t)>amount)
{
- begin=(*begin_iter->value_node)(blend_time).get(prev);
+ begin=(*begin_iter->value_node)(off_time).get(prev);
break;
}
}while(begin_iter!=iter);
-
+
// If we did not find a begin
if(begin.get_origin()==100.0f)
{
if(get_loop())
{
begin_iter=first_iter;
- begin=(*begin_iter->value_node)(blend_time).get(prev);
+ begin=(*begin_iter->value_node)(off_time).get(prev);
// begin=first;
}
else
{
// Writeme!
begin_iter=first_iter;
- begin=(*begin_iter->value_node)(blend_time).get(prev);
+ begin=(*begin_iter->value_node)(off_time).get(prev);
// begin=first;
}
}
-
+
etl::hermite<Vector> curve(begin.get_vertex(),end.get_vertex(),begin.get_tangent2(),end.get_tangent1());
etl::derivative< etl::hermite<Vector> > deriv(curve);
-
+
ret.set_vertex(curve(curr.get_origin()));
ret.set_width((end.get_width()-begin.get_width())*curr.get_origin()+begin.get_width());
ret.set_tangent1(deriv(curr.get_origin()));
ret.set_tangent2(deriv(curr.get_origin()));
-
+
float prev_tangent_scalar(1.0f);
float next_tangent_scalar(1.0f);
-
+
//synfig::info("index_%d:dist_from_begin=%d",index,dist_from_begin);
//synfig::info("index_%d:dist_from_end=%d",index,dist_from_end);
-
+
// If we are the next to the begin
if(begin_iter==--std::vector<ListEntry>::const_iterator(iter) || dist_from_begin==1)
{
else
{
float origin=curr.get_origin()-prev.get_origin();
- prev_tangent_scalar=(1.0f-origin)*amount+origin;
+ prev_tangent_scalar=(1.0f-origin)*amount+origin;
}
// If we are the next to the end
BLinePoint next;
next=((*(++std::vector<ListEntry>::const_iterator(iter))->value_node)(t).get(prev));
float origin=next.get_origin()-curr.get_origin();
- next_tangent_scalar=(1.0f-origin)*amount+origin;
+ next_tangent_scalar=(1.0f-origin)*amount+origin;
+ }
+ else
+ {
+ //! \todo this isn't quite right; we should handle looped blines identically no matter where the loop happens
+ //! and we currently don't. this at least makes it a lot better than it was before
+ float origin=end.get_origin()-curr.get_origin();
+ next_tangent_scalar=(1.0f-origin)*amount+origin;
}
next_scale=next_tangent_scalar;
-
+
//ret.set_vertex((curr.get_vertex()-ret.get_vertex())*amount+ret.get_vertex());
if(false)
{
// My first try
Point ref_point_begin(
(
- (*begin_iter->value_node)(begin_time).get(prev).get_vertex() +
- (*end_iter->value_node)(begin_time).get(prev).get_vertex()
+ (*begin_iter->value_node)(off_time).get(prev).get_vertex() +
+ (*end_iter->value_node)(off_time).get(prev).get_vertex()
) * 0.5
);
Point ref_point_end(
(
- (*begin_iter->value_node)(end_time).get(prev).get_vertex() +
- (*end_iter->value_node)(end_time).get(prev).get_vertex()
+ (*begin_iter->value_node)(on_time).get(prev).get_vertex() +
+ (*end_iter->value_node)(on_time).get(prev).get_vertex()
) * 0.5
);
Point ref_point_now(
else
{
// My second try
- Point begin_cord_sys[2], begin_cord_origin;
- Point end_cord_sys[2], end_cord_origin;
- Point curr_cord_sys[2], curr_cord_origin;
-
- {
- const Point a((*end_iter->value_node)(begin_time).get(prev).get_vertex());
- const Point b((*begin_iter->value_node)(begin_time).get(prev).get_vertex());
- begin_cord_origin=(a+b)/2;
- begin_cord_sys[0]=( b - a ).norm();
- begin_cord_sys[1]=begin_cord_sys[0].perp();
- }
- {
- const Point a((*end_iter->value_node)(end_time).get(prev).get_vertex());
- const Point b((*begin_iter->value_node)(end_time).get(prev).get_vertex());
- end_cord_origin=(a+b)/2;
- end_cord_sys[0]=( b - a ).norm();
- end_cord_sys[1]=end_cord_sys[0].perp();
- }
+
+ // define 3 coordinate systems:
+ Point off_coord_sys[2], off_coord_origin; // when the current vertex is completely off
+ Point on_coord_sys[2] , on_coord_origin; // when the current vertex is completely on
+ Point curr_coord_sys[2], curr_coord_origin; // the current state - somewhere in between
+
+ // for each of the 3 systems, the origin is half way between the previous and next active point
+ // and the axes are based on a vector from the next active point to the previous
{
- const Point a((*end_iter->value_node)(t).get(prev).get_vertex());
- const Point b((*begin_iter->value_node)(t).get(prev).get_vertex());
- curr_cord_origin=(a+b)/2;
- curr_cord_sys[0]=( b - a ).norm();
- curr_cord_sys[1]=curr_cord_sys[0].perp();
+ const Point end_pos_at_off_time(( *end_iter->value_node)(off_time).get(prev).get_vertex());
+ const Point begin_pos_at_off_time((*begin_iter->value_node)(off_time).get(prev).get_vertex());
+ off_coord_origin=(begin_pos_at_off_time + end_pos_at_off_time)/2;
+ off_coord_sys[0]=(begin_pos_at_off_time - end_pos_at_off_time).norm();
+ off_coord_sys[1]=off_coord_sys[0].perp();
+
+ const Point end_pos_at_on_time(( *end_iter->value_node)(on_time).get(prev).get_vertex());
+ const Point begin_pos_at_on_time((*begin_iter->value_node)(on_time).get(prev).get_vertex());
+ on_coord_origin=(begin_pos_at_on_time + end_pos_at_on_time)/2;
+ on_coord_sys[0]=(begin_pos_at_on_time - end_pos_at_on_time).norm();
+ on_coord_sys[1]=on_coord_sys[0].perp();
+
+ const Point end_pos_at_current_time(( *end_iter->value_node)(t).get(prev).get_vertex());
+ const Point begin_pos_at_current_time((*begin_iter->value_node)(t).get(prev).get_vertex());
+ curr_coord_origin=(begin_pos_at_current_time + end_pos_at_current_time)/2;
+ curr_coord_sys[0]=(begin_pos_at_current_time - end_pos_at_current_time).norm();
+ curr_coord_sys[1]=curr_coord_sys[0].perp();
}
-
- /*
- end_cord_origin=(*end_iter->value_node)(end_time).get(prev).get_vertex();
- end_cord_sys[0]=(
- (*begin_iter->value_node)(end_time).get(prev).get_vertex() -
- end_cord_origin
- ).norm();
- end_cord_sys[1]=end_cord_sys[0].perp();
-
- curr_cord_origin=(*end_iter->value_node)(t).get(prev).get_vertex();
- curr_cord_sys[0]=(
- (*begin_iter->value_node)(t).get(prev).get_vertex() -
- curr_cord_origin
- ).norm();
- curr_cord_sys[1]=curr_cord_sys[0].perp();
- */
-
- // Convert start point
- Point a;
- Vector at1,at2;
+
+ /* The code that was here before used just end_iter as the origin, rather than the mid-point */
+
+ // For each of the 3 coordinate systems we've just defined, we convert a point and tangent(s) into that system
+
+ // Convert point where vertex is fully 'off'
+ Point trans_off_point;
+ Vector trans_off_t1,trans_off_t2;
{
- Point tmp(ret.get_vertex()-begin_cord_origin);
- a[0]=tmp*begin_cord_sys[0];
- a[1]=tmp*begin_cord_sys[1];
+ Point tmp(ret.get_vertex()-off_coord_origin);
+ trans_off_point[0]=tmp*off_coord_sys[0];
+ trans_off_point[1]=tmp*off_coord_sys[1];
#define COORD_SYS_RADIAL_TAN_INTERP 1
-
+
#ifdef COORD_SYS_RADIAL_TAN_INTERP
- tmp=ret.get_tangent1()+ret.get_vertex()-begin_cord_origin;
- at1[0]=tmp*begin_cord_sys[0];
- at1[1]=tmp*begin_cord_sys[1];
-
+ tmp=ret.get_tangent1()+ret.get_vertex()-off_coord_origin;
+ trans_off_t1[0]=tmp*off_coord_sys[0];
+ trans_off_t1[1]=tmp*off_coord_sys[1];
+
if(curr.get_split_tangent_flag())
{
- tmp=ret.get_tangent2()+ret.get_vertex()-begin_cord_origin;
- at2[0]=tmp*begin_cord_sys[0];
- at2[1]=tmp*begin_cord_sys[1];
+ tmp=ret.get_tangent2()+ret.get_vertex()-off_coord_origin;
+ trans_off_t2[0]=tmp*off_coord_sys[0];
+ trans_off_t2[1]=tmp*off_coord_sys[1];
}
#endif
}
-
- // Convert finish point
- Point b;
- Vector bt1,bt2;
+
+ // Convert point where vertex is fully 'on'
+ Point trans_on_point;
+ Vector trans_on_t1,trans_on_t2;
{
- Point tmp(curr.get_vertex()-end_cord_origin);
- b[0]=tmp*end_cord_sys[0];
- b[1]=tmp*end_cord_sys[1];
-
+ Point tmp(curr.get_vertex()-on_coord_origin);
+ trans_on_point[0]=tmp*on_coord_sys[0];
+ trans_on_point[1]=tmp*on_coord_sys[1];
+
#ifdef COORD_SYS_RADIAL_TAN_INTERP
- tmp=curr.get_tangent1()+curr.get_vertex()-end_cord_origin;
- bt1[0]=tmp*end_cord_sys[0];
- bt1[1]=tmp*end_cord_sys[1];
+ tmp=curr.get_tangent1()+curr.get_vertex()-on_coord_origin;
+ trans_on_t1[0]=tmp*on_coord_sys[0];
+ trans_on_t1[1]=tmp*on_coord_sys[1];
if(curr.get_split_tangent_flag())
{
- tmp=curr.get_tangent2()+curr.get_vertex()-end_cord_origin;
- bt2[0]=tmp*end_cord_sys[0];
- bt2[1]=tmp*end_cord_sys[1];
+ tmp=curr.get_tangent2()+curr.get_vertex()-on_coord_origin;
+ trans_on_t2[0]=tmp*on_coord_sys[0];
+ trans_on_t2[1]=tmp*on_coord_sys[1];
}
#endif
}
-
+
// Convert current point
- Point c;
- Vector ct1,ct2;
+ Point trans_curr_point;
+ Vector trans_curr_t1,trans_curr_t2;
{
// Transpose (invert)
- swap(curr_cord_sys[0][1],curr_cord_sys[1][0]);
+ swap(curr_coord_sys[0][1],curr_coord_sys[1][0]);
- Point tmp((b-a)*amount+a);
- c[0]=tmp*curr_cord_sys[0];
- c[1]=tmp*curr_cord_sys[1];
- c+=curr_cord_origin;
+ Point tmp((trans_on_point-trans_off_point)*amount+trans_off_point);
+ trans_curr_point[0]=tmp*curr_coord_sys[0];
+ trans_curr_point[1]=tmp*curr_coord_sys[1];
+ trans_curr_point+=curr_coord_origin;
#define INTERP_FUNCTION radial_interpolation
//#define INTERP_FUNCTION linear_interpolation
-
+
#ifdef COORD_SYS_RADIAL_TAN_INTERP
- tmp=INTERP_FUNCTION(at1,bt1,amount);
- ct1[0]=tmp*curr_cord_sys[0];
- ct1[1]=tmp*curr_cord_sys[1];
- ct1+=curr_cord_origin;
- ct1-=c;
-
+ tmp=INTERP_FUNCTION(trans_off_t1,trans_on_t1,amount);
+ trans_curr_t1[0]=tmp*curr_coord_sys[0];
+ trans_curr_t1[1]=tmp*curr_coord_sys[1];
+ trans_curr_t1+=curr_coord_origin;
+ trans_curr_t1-=trans_curr_point;
+
if(curr.get_split_tangent_flag())
{
- tmp=INTERP_FUNCTION(at2,bt2,amount);
- ct2[0]=tmp*curr_cord_sys[0];
- ct2[1]=tmp*curr_cord_sys[1];
- ct2+=curr_cord_origin;
- ct2-=c;
+ tmp=INTERP_FUNCTION(trans_off_t2,trans_on_t2,amount);
+ trans_curr_t2[0]=tmp*curr_coord_sys[0];
+ trans_curr_t2[1]=tmp*curr_coord_sys[1];
+ trans_curr_t2+=curr_coord_origin;
+ trans_curr_t2-=trans_curr_point;
}
#endif
}
- ret.set_vertex(c);
+ ret.set_vertex(trans_curr_point);
#ifndef COORD_SYS_RADIAL_TAN_INTERP
ret.set_tangent1(radial_interpolation(ret.get_tangent1(),curr.get_tangent1(),amount));
ret.set_split_tangent_flag(curr.get_split_tangent_flag());
if(ret.get_split_tangent_flag())
ret.set_tangent2(radial_interpolation(ret.get_tangent2(),curr.get_tangent2(),amount));
-#else
- ret.set_tangent1(ct1);
+#else
+ ret.set_tangent1(trans_curr_t1);
ret.set_split_tangent_flag(curr.get_split_tangent_flag());
if(ret.get_split_tangent_flag())
- ret.set_tangent2(ct2);
+ ret.set_tangent2(trans_curr_t2);
#endif
}
-
+
ret.set_origin(curr.get_origin());
ret.set_width((curr.get_width()-ret.get_width())*amount+ret.get_width());
-
+
// Handle the case where we are the first vertex
if(first_flag)
prev=ret;
}
}
-
+
if(next_scale!=1.0f)
{
ret_list.back().set_split_tangent_flag(true);
);
}
*/
-
+
if(list.empty())
synfig::warning(string("ValueNode_BLine::operator()():")+_("No entries in list"));
else
{
assert(i>=0 && (unsigned)i<list.size());
return etl::strprintf(_("Vertex %03d"),i+1);
-}
+}
ValueNode*
ValueNode_BLine::clone(const GUID& deriv_guid)const
ValueNode_BLine* ret=new ValueNode_BLine();
ret->set_guid(get_guid()^deriv_guid);
-
+
std::vector<ListEntry>::const_iterator iter;
for(iter=list.begin();iter!=list.end();++iter)
}
}
ret->set_loop(get_loop());
-
+
return ret;
}
projects / synfig.git / blobdiff