416 |
|
rev_hessian(hesscol); |
417 |
|
add2hessian(hessian, hessrow[j], hessdia, hesscol, backg); |
418 |
|
} |
419 |
< |
if (gradient != NULL) { |
419 |
> |
if (gradrow != NULL) { |
420 |
|
comp_gradient(graddia, &fftr, hp->rp->ron); |
421 |
|
rev_gradient(gradcol); |
422 |
|
add2gradient(gradient, gradrow[j], graddia, gradcol, backg); |
478 |
|
VSUB(vd, ap->p, hp->rp->rop); |
479 |
|
/* brightness over cosine factor */ |
480 |
|
gfact = colval(ap->v,CIEY) / DOT(hp->rp->ron, vd); |
481 |
< |
/* -sine = -proj_radius/vd_length */ |
482 |
< |
dgsum[0] += DOT(uv[1], vd) * gfact; |
483 |
< |
dgsum[1] -= DOT(uv[0], vd) * gfact; |
481 |
> |
/* sine = proj_radius/vd_length */ |
482 |
> |
dgsum[0] -= DOT(uv[1], vd) * gfact; |
483 |
> |
dgsum[1] += DOT(uv[0], vd) * gfact; |
484 |
|
} |
485 |
|
dg[0] = dgsum[0] / (hp->ns*hp->ns); |
486 |
|
dg[1] = dgsum[1] / (hp->ns*hp->ns); |