function y = d_abs_cos( x )
  y  = -sin(x) ;
  ix = cos(x) < 0 ;
  y( ix ) = sin( x(ix) ) ;
end