#if 0
#elif 0
 aa_mul(f,x,y);				/* xy= x*y */
 aa_sub(f,f,z);				/* f=f-z=xy-z */
#elif 0
 aa_sqr(x2,x);				/* x2= x*x */
aa_debug(x2);
 aa_sqr(y2,y);				/* y2= y*y */
aa_debug(y2);
 aa_add(f,x2,y2);			/* f=x2+y2 */
aa_debug(f);
 aa_sub(f,f,z);				/* f=f-z=x2+y2-z */
aa_debug(f);
#elif 1
 /* (x^2 + (y-1)^2 + z^2 - 2.25) * (x^2 + (y+1)^2 + z^2 - 0.6) */
 aa_sqr(x2,x);
 aa_sqr(z2,z);
 aa_trans(y2,y,-1);
 aa_sqr(y2,y2);
 aa_add(f,x2,y2);
 aa_add(f,f,z2);
 aa_trans(f,f,-2.25);
 aa_trans(y2,y,1);
 aa_sqr(y2,y2);
 aa_add(a,x2,y2);
 aa_add(a,a,z2);
 aa_trans(a,a,-0.6);
 aa_mul(f,f,a);
#elif 0
 aa_sqr(x2,x);				/* x2= x*x */
 aa_sqr(y2,y);				/* y2= y*y */
 aa_mul(xy,x,y);                        /* xy= x*y */
 aa_sqr(xy2,xy);			/* xy2= x*y*x*y */
 aa_scale(xy2,xy2,0.5);
 aa_add(f,x2,y2);                       /* f=x2+y2 */
 aa_add(f,f,xy);                        /* f=f+xy=x2+y2+xy */
 aa_sub(f,f,xy2);                       /* f=f-xy2=x2+y2+xy-(xy)2/2 */
 aa_sub(f,f,z);				/* f=z-f=z-x2-y2 */
 f[0]-=0.25;
#elif 0
 aa_sqr(x2,x);				/* x2= x*x */
 aa_sqr(y2,y);				/* y2= y*y */
 aa_mul(x3,x,x2);                       /* x3= x*x2 */
 aa_sub(f,y2,x3);                       /* f=y2-x3 */
 aa_add(f,f,x);                         /* f=f+x=y2-x3+x */
 aa_sub(f,f,z);				/* f=f-z */
#elif 1					/* heart */
 aa_sqr(x2,x);
 aa_sqr(y2,y);
 aa_sqr(z2,z);
 aa_mul(z3,z2,z);
 aa_scale(a,x2,2);
 aa_add(a,a,y2);
 aa_add(a,a,z2);
 aa_trans(a,a,-1);
 aa_sqr(f,a);
 aa_mul(f,f,a);
 aa_mul(a,x2,z3);
 aa_scale(a,a,.1);
 aa_sub(f,f,a);
 aa_mul(a,y2,z3);
 aa_sub(f,f,a);
#else
#error no function defined
#endif