source/blender/render/intern/source/shadeinput.c

   1 /*
   2 * $Id$
   3  *
   4  * ***** BEGIN GPL LICENSE BLOCK *****
   5  *
   6  * This program is free software; you can redistribute it and/or
   7  * modify it under the terms of the GNU General Public License
   8  * as published by the Free Software Foundation; either version 2
   9  * of the License, or (at your option) any later version.
  10  *
  11  * This program is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14  * GNU General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU General Public License
  17  * along with this program; if not, write to the Free Software Foundation,
  18  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  19  *
  20  * The Original Code is Copyright (C) 2006 Blender Foundation
  21  * All rights reserved.
  22  *
  23  * Contributors: Hos, Robert Wenzlaff.
  24  *
  25  * ***** END GPL LICENSE BLOCK *****
  26  */
  27
  28 /** \file blender/render/intern/source/shadeinput.c
  29  *  \ingroup render
  30  */
  31
  32
  33 #include <stdio.h>
  34 #include <math.h>
  35 #include <string.h>
  36
  37
  38 #include "BLI_math.h"
  39 #include "BLI_blenlib.h"
  40 #include "BLI_utildefines.h"
  41
  42 #include "DNA_curve_types.h"
  43 #include "DNA_group_types.h"
  44 #include "DNA_lamp_types.h"
  45 #include "DNA_meshdata_types.h"
  46 #include "DNA_material_types.h"
  47
  48 #include "BKE_colortools.h"
  49
  50 #include "BKE_node.h"
  51
  52 /* local include */
  53 #include "raycounter.h"
  54 #include "renderpipeline.h"
  55 #include "render_types.h"
  56 #include "renderdatabase.h"
  57 #include "rendercore.h"
  58 #include "shadbuf.h"
  59 #include "shading.h"
  60 #include "strand.h"
  61 #include "texture.h"
  62 #include "volumetric.h"
  63 #include "zbuf.h"
  64
  65 /* ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ */
  66 /* defined in pipeline.c, is hardcopy of active dynamic allocated Render */
  67 /* only to be used here in this file, it's for speed */
  68 extern struct Render R;
  69 /* ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ */
  70
  71
  72 #define VECADDISFAC(v1,v3,fac) {*(v1)+= *(v3)*(fac); *(v1+1)+= *(v3+1)*(fac); *(v1+2)+= *(v3+2)*(fac);}
  73
  74
  75
  76 /* Shade Sample order:
  77
  78 - shade_samples_fill_with_ps()
  79         - for each sample
  80                 - shade_input_set_triangle()  <- if prev sample-face is same, use shade_input_copy_triangle()
  81                 - if vlr
  82                         - shade_input_set_viewco()    <- not for ray or bake
  83                         - shade_input_set_uv()        <- not for ray or bake
  84                         - shade_input_set_normals()
  85 - shade_samples()
  86         - if AO
  87                 - shade_samples_do_AO()
  88         - if shading happens
  89                 - for each sample
  90                         - shade_input_set_shade_texco()
  91                         - shade_samples_do_shade()
  92 - OSA: distribute sample result with filter masking
  93
  94         */
  95
  96 /* initialise material variables in shadeinput,
  97  * doing inverse gamma correction where applicable */
  98 void shade_input_init_material(ShadeInput *shi)
  99 {
 100         /* note, keep this synced with render_types.h */
 101         memcpy(&shi->r, &shi->mat->r, 23*sizeof(float));
 102         shi->har= shi->mat->har;
 103 }
 104
 105 /* also used as callback for nodes */
 106 /* delivers a fully filled in ShadeResult, for all passes */
 107 void shade_material_loop(ShadeInput *shi, ShadeResult *shr)
 108 {
 109
 110         shade_lamp_loop(shi, shr);      /* clears shr */
 111
 112         if(shi->translucency!=0.0f) {
 113                 ShadeResult shr_t;
 114                 float fac= shi->translucency;
 115
 116                 shade_input_init_material(shi);
 117
 118                 VECCOPY(shi->vn, shi->vno);
 119                 VECMUL(shi->vn, -1.0f);
 120                 VECMUL(shi->facenor, -1.0f);
 121                 shi->depth++;   /* hack to get real shadow now */
 122                 shade_lamp_loop(shi, &shr_t);
 123                 shi->depth--;
 124
 125                 /* a couple of passes */
 126                 VECADDISFAC(shr->combined, shr_t.combined, fac);
 127                 if(shi->passflag & SCE_PASS_SPEC)
 128                         VECADDISFAC(shr->spec, shr_t.spec, fac);
 129                 if(shi->passflag & SCE_PASS_DIFFUSE)
 130                         VECADDISFAC(shr->diff, shr_t.diff, fac);
 131                 if(shi->passflag & SCE_PASS_SHADOW)
 132                         VECADDISFAC(shr->shad, shr_t.shad, fac);
 133
 134                 VECMUL(shi->vn, -1.0f);
 135                 VECMUL(shi->facenor, -1.0f);
 136         }
 137
 138         /* depth >= 1 when ray-shading */
 139         if(shi->depth==0 || shi->volume_depth > 0) {
 140                 if(R.r.mode & R_RAYTRACE) {
 141                         if(shi->ray_mirror!=0.0f || ((shi->mode & MA_TRANSP) && (shi->mode & MA_RAYTRANSP) && shr->alpha!=1.0f)) {
 142                                 /* ray trace works on combined, but gives pass info */
 143                                 ray_trace(shi, shr);
 144                         }
 145                 }
 146                 /* disable adding of sky for raytransp */
 147                 if((shi->mode & MA_TRANSP) && (shi->mode & MA_RAYTRANSP))
 148                         if((shi->layflag & SCE_LAY_SKY) && (R.r.alphamode==R_ADDSKY))
 149                                 shr->alpha= 1.0f;
 150         }
 151
 152         if(R.r.mode & R_RAYTRACE) {
 153                 if (R.render_volumes_inside.first)
 154                         shade_volume_inside(shi, shr);
 155         }
 156 }
 157
 158
 159 /* do a shade, finish up some passes, apply mist */
 160 void shade_input_do_shade(ShadeInput *shi, ShadeResult *shr)
 161 {
 162         float alpha;
 163
 164         /* ------  main shading loop -------- */
 165 #ifdef RE_RAYCOUNTER
 166         memset(&shi->raycounter, 0, sizeof(shi->raycounter));
 167 #endif
 168
 169         if(shi->mat->nodetree && shi->mat->use_nodes) {
 170                 ntreeShaderExecTree(shi->mat->nodetree, shi, shr);
 171         }
 172         else {
 173                 /* copy all relevant material vars, note, keep this synced with render_types.h */
 174                 shade_input_init_material(shi);
 175
 176                 if (shi->mat->material_type == MA_TYPE_VOLUME) {
 177                         if(R.r.mode & R_RAYTRACE) {
 178                                 shade_volume_outside(shi, shr);
 179                         }
 180                 } else { /* MA_TYPE_SURFACE, MA_TYPE_WIRE */
 181                         shade_material_loop(shi, shr);
 182                 }
 183         }
 184
 185         /* copy additional passes */
 186         if(shi->passflag & (SCE_PASS_VECTOR|SCE_PASS_NORMAL)) {
 187                 QUATCOPY(shr->winspeed, shi->winspeed);
 188                 VECCOPY(shr->nor, shi->vn);
 189         }
 190
 191         /* MIST */
 192         if((shi->passflag & SCE_PASS_MIST) || ((R.wrld.mode & WO_MIST) && (shi->mat->mode & MA_NOMIST)==0))  {
 193                 if(R.r.mode & R_ORTHO)
 194                         shr->mist= mistfactor(-shi->co[2], shi->co);
 195                 else
 196                         shr->mist= mistfactor(len_v3(shi->co), shi->co);
 197         }
 198         else shr->mist= 0.0f;
 199
 200         if((R.wrld.mode & WO_MIST) && (shi->mat->mode & MA_NOMIST)==0 ) {
 201                 alpha= shr->mist;
 202         }
 203         else alpha= 1.0f;
 204
 205         /* add mist and premul color */
 206         if(shr->alpha!=1.0f || alpha!=1.0f) {
 207                 float fac= alpha*(shr->alpha);
 208                 shr->combined[3]= fac;
 209
 210                 if (shi->mat->material_type!= MA_TYPE_VOLUME)
 211                         mul_v3_fl(shr->combined, fac);
 212         }
 213         else
 214                 shr->combined[3]= 1.0f;
 215
 216         /* add z */
 217         shr->z= -shi->co[2];
 218
 219         /* RAYHITS */
 220 /*
 221         if(1 || shi->passflag & SCE_PASS_RAYHITS)
 222         {
 223                 shr->rayhits[0] = (float)shi->raycounter.faces.test;
 224                 shr->rayhits[1] = (float)shi->raycounter.bb.hit;
 225                 shr->rayhits[2] = 0.0;
 226                 shr->rayhits[3] = 1.0;
 227         }
 228  */
 229         RE_RC_MERGE(&re_rc_counter[shi->thread], &shi->raycounter);
 230 }
 231
 232 /* **************************************************************************** */
 233 /*                    ShadeInput                                                */
 234 /* **************************************************************************** */
 235
 236
 237 void vlr_set_uv_indices(VlakRen *vlr, int *i1, int *i2, int *i3)
 238 {
 239         /* to prevent storing new tfaces or vcols, we check a split runtime */
 240         /*              4---3           4---3 */
 241         /*              |\ 1|   or  |1 /| */
 242         /*              |0\ |           |/ 0| */
 243         /*              1---2           1---2   0 = orig face, 1 = new face */
 244
 245         /* Update vert nums to point to correct verts of original face */
 246         if(vlr->flag & R_DIVIDE_24) {
 247                 if(vlr->flag & R_FACE_SPLIT) {
 248                         (*i1)++; (*i2)++; (*i3)++;
 249                 }
 250                 else {
 251                         (*i3)++;
 252                 }
 253         }
 254         else if(vlr->flag & R_FACE_SPLIT) {
 255                 (*i2)++; (*i3)++;
 256         }
 257 }
 258
 259 /* copy data from face to ShadeInput, general case */
 260 /* indices 0 1 2 3 only */
 261 void shade_input_set_triangle_i(ShadeInput *shi, ObjectInstanceRen *obi, VlakRen *vlr, short i1, short i2, short i3)
 262 {
 263         VertRen **vpp= &vlr->v1;
 264
 265         shi->vlr= vlr;
 266         shi->obi= obi;
 267         shi->obr= obi->obr;
 268
 269         shi->v1= vpp[i1];
 270         shi->v2= vpp[i2];
 271         shi->v3= vpp[i3];
 272
 273         shi->i1= i1;
 274         shi->i2= i2;
 275         shi->i3= i3;
 276
 277         /* note, shi->mat is set in node shaders */
 278         shi->mat= shi->mat_override?shi->mat_override:vlr->mat;
 279
 280         shi->osatex= (shi->mat->texco & TEXCO_OSA);
 281         shi->mode= shi->mat->mode_l;            /* or-ed result for all nodes */
 282
 283         /* facenormal copy, can get flipped */
 284         shi->flippednor= 0;
 285         RE_vlakren_get_normal(&R, obi, vlr, shi->facenor);
 286
 287         /* calculate vertexnormals */
 288         if(vlr->flag & R_SMOOTH) {
 289                 VECCOPY(shi->n1, shi->v1->n);
 290                 VECCOPY(shi->n2, shi->v2->n);
 291                 VECCOPY(shi->n3, shi->v3->n);
 292
 293                 if(obi->flag & R_TRANSFORMED) {
 294                         mul_m3_v3(obi->nmat, shi->n1); normalize_v3(shi->n1);
 295                         mul_m3_v3(obi->nmat, shi->n2); normalize_v3(shi->n2);
 296                         mul_m3_v3(obi->nmat, shi->n3); normalize_v3(shi->n3);
 297                 }
 298         }
 299 }
 300
 301 /* note, facenr declared volatile due to over-eager -O2 optimizations
 302  * on cygwin (particularly -frerun-cse-after-loop)
 303  */
 304
 305 /* copy data from face to ShadeInput, scanline case */
 306 void shade_input_set_triangle(ShadeInput *shi, volatile int obi, volatile int facenr, int UNUSED(normal_flip))
 307 {
 308         if(facenr>0) {
 309                 shi->obi= &R.objectinstance[obi];
 310                 shi->obr= shi->obi->obr;
 311                 shi->facenr= (facenr-1) & RE_QUAD_MASK;
 312                 if( shi->facenr < shi->obr->totvlak ) {
 313                         VlakRen *vlr= RE_findOrAddVlak(shi->obr, shi->facenr);
 314
 315                         if(facenr & RE_QUAD_OFFS)
 316                                 shade_input_set_triangle_i(shi, shi->obi, vlr, 0, 2, 3);
 317                         else
 318                                 shade_input_set_triangle_i(shi, shi->obi, vlr, 0, 1, 2);
 319                 }
 320                 else
 321                         shi->vlr= NULL; /* general signal we got sky */
 322         }
 323         else
 324                 shi->vlr= NULL; /* general signal we got sky */
 325 }
 326
 327 /* full osa case: copy static info */
 328 void shade_input_copy_triangle(ShadeInput *shi, ShadeInput *from)
 329 {
 330         /* not so nice, but works... warning is in RE_shader_ext.h */
 331         memcpy(shi, from, sizeof(struct ShadeInputCopy));
 332 }
 333
 334 /* copy data from strand to shadeinput */
 335 void shade_input_set_strand(ShadeInput *shi, StrandRen *strand, StrandPoint *spoint)
 336 {
 337         /* note, shi->mat is set in node shaders */
 338         shi->mat= shi->mat_override? shi->mat_override: strand->buffer->ma;
 339
 340         shi->osatex= (shi->mat->texco & TEXCO_OSA);
 341         shi->mode= shi->mat->mode_l;            /* or-ed result for all nodes */
 342
 343         /* shade_input_set_viewco equivalent */
 344         VECCOPY(shi->co, spoint->co);
 345         VECCOPY(shi->view, shi->co);
 346         normalize_v3(shi->view);
 347
 348         shi->xs= (int)spoint->x;
 349         shi->ys= (int)spoint->y;
 350
 351         if(shi->osatex || (R.r.mode & R_SHADOW)) {
 352                 VECCOPY(shi->dxco, spoint->dtco);
 353                 VECCOPY(shi->dyco, spoint->dsco);
 354         }
 355
 356         /* dxview, dyview, not supported */
 357
 358         /* facenormal, simply viewco flipped */
 359         VECCOPY(shi->facenor, spoint->nor);
 360
 361         /* shade_input_set_normals equivalent */
 362         if(shi->mat->mode & MA_TANGENT_STR) {
 363                 VECCOPY(shi->vn, spoint->tan)
 364         }
 365         else {
 366                 float cross[3];
 367
 368                 cross_v3_v3v3(cross, spoint->co, spoint->tan);
 369                 cross_v3_v3v3(shi->vn, cross, spoint->tan);
 370                 normalize_v3(shi->vn);
 371
 372                 if(INPR(shi->vn, shi->view) < 0.0f)
 373                         negate_v3(shi->vn);
 374         }
 375
 376         VECCOPY(shi->vno, shi->vn);
 377 }
 378
 379 void shade_input_set_strand_texco(ShadeInput *shi, StrandRen *strand, StrandVert *svert, StrandPoint *spoint)
 380 {
 381         StrandBuffer *strandbuf= strand->buffer;
 382         ObjectRen *obr= strandbuf->obr;
 383         StrandVert *sv;
 384         int mode= shi->mode;            /* or-ed result for all nodes */
 385         short texco= shi->mat->texco;
 386
 387         if((shi->mat->texco & TEXCO_REFL)) {
 388                 /* shi->dxview, shi->dyview, not supported */
 389         }
 390
 391         if(shi->osatex && (texco & (TEXCO_NORM|TEXCO_REFL))) {
 392                 /* not supported */
 393         }
 394
 395         if(mode & (MA_TANGENT_V|MA_NORMAP_TANG)) {
 396                 VECCOPY(shi->tang, spoint->tan);
 397                 VECCOPY(shi->nmaptang, spoint->tan);
 398         }
 399
 400         if(mode & MA_STR_SURFDIFF) {
 401                 float *surfnor= RE_strandren_get_surfnor(obr, strand, 0);
 402
 403                 if(surfnor)
 404                         VECCOPY(shi->surfnor, surfnor)
 405                 else
 406                         VECCOPY(shi->surfnor, shi->vn)
 407
 408                 if(shi->mat->strand_surfnor > 0.0f) {
 409                         shi->surfdist= 0.0f;
 410                         for(sv=strand->vert; sv!=svert; sv++)
 411                                 shi->surfdist+=len_v3v3(sv->co, (sv+1)->co);
 412                         shi->surfdist += spoint->t*len_v3v3(sv->co, (sv+1)->co);
 413                 }
 414         }
 415
 416         if(R.r.mode & R_SPEED) {
 417                 float *speed;
 418
 419                 speed= RE_strandren_get_winspeed(shi->obi, strand, 0);
 420                 if(speed)
 421                         QUATCOPY(shi->winspeed, speed)
 422                 else
 423                         shi->winspeed[0]= shi->winspeed[1]= shi->winspeed[2]= shi->winspeed[3]= 0.0f;
 424         }
 425
 426         /* shade_input_set_shade_texco equivalent */
 427         if(texco & NEED_UV) {
 428                 if(texco & TEXCO_ORCO) {
 429                         VECCOPY(shi->lo, strand->orco);
 430                         /* no shi->osatex, orco derivatives are zero */
 431                 }
 432
 433                 if(texco & TEXCO_GLOB) {
 434                         VECCOPY(shi->gl, shi->co);
 435                         mul_m4_v3(R.viewinv, shi->gl);
 436
 437                         if(shi->osatex) {
 438                                 VECCOPY(shi->dxgl, shi->dxco);
 439                                 mul_mat3_m4_v3(R.viewinv, shi->dxgl);
 440                                 VECCOPY(shi->dygl, shi->dyco);
 441                                 mul_mat3_m4_v3(R.viewinv, shi->dygl);
 442                         }
 443                 }
 444
 445                 if(texco & TEXCO_STRAND) {
 446                         shi->strandco= spoint->strandco;
 447
 448                         if(shi->osatex) {
 449                                 shi->dxstrand= spoint->dtstrandco;
 450                                 shi->dystrand= 0.0f;
 451                         }
 452                 }
 453
 454                 if((texco & TEXCO_UV) || (mode & (MA_VERTEXCOL|MA_VERTEXCOLP|MA_FACETEXTURE)))  {
 455                         MCol *mcol;
 456                         float *uv;
 457                         char *name;
 458                         int i;
 459
 460                         shi->totuv= 0;
 461                         shi->totcol= 0;
 462                         shi->actuv= obr->actmtface;
 463                         shi->actcol= obr->actmcol;
 464
 465                         if(mode & (MA_VERTEXCOL|MA_VERTEXCOLP)) {
 466                                 for (i=0; (mcol=RE_strandren_get_mcol(obr, strand, i, &name, 0)); i++) {
 467                                         ShadeInputCol *scol= &shi->col[i];
 468                                         char *cp= (char*)mcol;
 469
 470                                         shi->totcol++;
 471                                         scol->name= name;
 472
 473                                         scol->col[0]= cp[3]/255.0f;
 474                                         scol->col[1]= cp[2]/255.0f;
 475                                         scol->col[2]= cp[1]/255.0f;
 476                                         scol->col[3]= cp[0]/255.0f;
 477                                 }
 478
 479                                 if(shi->totcol) {
 480                                         shi->vcol[0]= shi->col[shi->actcol].col[0];
 481                                         shi->vcol[1]= shi->col[shi->actcol].col[1];
 482                                         shi->vcol[2]= shi->col[shi->actcol].col[2];
 483                                         shi->vcol[3]= shi->col[shi->actcol].col[3];
 484                                 }
 485                                 else {
 486                                         shi->vcol[0]= 0.0f;
 487                                         shi->vcol[1]= 0.0f;
 488                                         shi->vcol[2]= 0.0f;
 489                                         shi->vcol[3]= 0.0f;
 490                                 }
 491                         }
 492
 493                         for (i=0; (uv=RE_strandren_get_uv(obr, strand, i, &name, 0)); i++) {
 494                                 ShadeInputUV *suv= &shi->uv[i];
 495
 496                                 shi->totuv++;
 497                                 suv->name= name;
 498
 499                                 if(strandbuf->overrideuv == i) {
 500                                         suv->uv[0]= -1.0f;
 501                                         suv->uv[1]= spoint->strandco;
 502                                         suv->uv[2]= 0.0f;
 503                                 }
 504                                 else {
 505                                         suv->uv[0]= -1.0f + 2.0f*uv[0];
 506                                         suv->uv[1]= -1.0f + 2.0f*uv[1];
 507                                         suv->uv[2]= 0.0f;       /* texture.c assumes there are 3 coords */
 508                                 }
 509
 510                                 if(shi->osatex) {
 511                                         suv->dxuv[0]= 0.0f;
 512                                         suv->dxuv[1]= 0.0f;
 513                                         suv->dyuv[0]= 0.0f;
 514                                         suv->dyuv[1]= 0.0f;
 515                                 }
 516
 517                                 if((mode & MA_FACETEXTURE) && i==obr->actmtface) {
 518                                         if((mode & (MA_VERTEXCOL|MA_VERTEXCOLP))==0) {
 519                                                 shi->vcol[0]= 1.0f;
 520                                                 shi->vcol[1]= 1.0f;
 521                                                 shi->vcol[2]= 1.0f;
 522                                                 shi->vcol[3]= 1.0f;
 523                                         }
 524                                 }
 525                         }
 526
 527                         if(shi->totuv == 0) {
 528                                 ShadeInputUV *suv= &shi->uv[0];
 529
 530                                 suv->uv[0]= 0.0f;
 531                                 suv->uv[1]= spoint->strandco;
 532                                 suv->uv[2]= 0.0f;       /* texture.c assumes there are 3 coords */
 533
 534                                 if(mode & MA_FACETEXTURE) {
 535                                         /* no tface? set at 1.0f */
 536                                         shi->vcol[0]= 1.0f;
 537                                         shi->vcol[1]= 1.0f;
 538                                         shi->vcol[2]= 1.0f;
 539                                         shi->vcol[3]= 1.0f;
 540                                 }
 541                         }
 542
 543                 }
 544
 545                 if(texco & TEXCO_NORM) {
 546                         shi->orn[0]= -shi->vn[0];
 547                         shi->orn[1]= -shi->vn[1];
 548                         shi->orn[2]= -shi->vn[2];
 549                 }
 550
 551                 if(texco & TEXCO_STRESS) {
 552                         /* not supported */
 553                 }
 554
 555                 if(texco & TEXCO_TANGENT) {
 556                         if((mode & MA_TANGENT_V)==0) {
 557                                 /* just prevent surprises */
 558                                 shi->tang[0]= shi->tang[1]= shi->tang[2]= 0.0f;
 559                                 shi->nmaptang[0]= shi->nmaptang[1]= shi->nmaptang[2]= 0.0f;
 560                         }
 561                 }
 562         }
 563
 564         /* this only avalailable for scanline renders */
 565         if(shi->depth==0) {
 566                 if(texco & TEXCO_WINDOW) {
 567                         shi->winco[0]= -1.0f + 2.0f*spoint->x/(float)R.winx;
 568                         shi->winco[1]= -1.0f + 2.0f*spoint->y/(float)R.winy;
 569                         shi->winco[2]= 0.0f;
 570
 571                         /* not supported */
 572                         if(shi->osatex) {
 573                                 shi->dxwin[0]= 0.0f;
 574                                 shi->dywin[1]= 0.0f;
 575                                 shi->dxwin[0]= 0.0f;
 576                                 shi->dywin[1]= 0.0f;
 577                         }
 578                 }
 579
 580                 if(texco & TEXCO_STICKY) {
 581                         /* not supported */
 582                 }
 583         }
 584
 585         if (shi->do_manage) {
 586                 if(mode & (MA_VERTEXCOL|MA_VERTEXCOLP|MA_FACETEXTURE)) {
 587                         srgb_to_linearrgb_v3_v3(shi->vcol, shi->vcol);
 588                 }
 589         }
 590
 591 }
 592
 593 /* from scanline pixel coordinates to 3d coordinates, requires set_triangle */
 594 void shade_input_calc_viewco(ShadeInput *shi, float x, float y, float z, float *view, float *dxyview, float *co, float *dxco, float *dyco)
 595 {
 596         /* returns not normalized, so is in viewplane coords */
 597         calc_view_vector(view, x, y);
 598
 599         if(shi->mat->material_type == MA_TYPE_WIRE) {
 600                 /* wire cannot use normal for calculating shi->co, so
 601                  * we reconstruct the coordinate less accurate */
 602                 if(R.r.mode & R_ORTHO)
 603                         calc_renderco_ortho(co, x, y, z);
 604                 else
 605                         calc_renderco_zbuf(co, view, z);
 606         }
 607         else {
 608                 /* for non-wire, intersect with the triangle to get the exact coord */
 609                 float fac, dface, v1[3];
 610
 611                 VECCOPY(v1, shi->v1->co);
 612                 if(shi->obi->flag & R_TRANSFORMED)
 613                         mul_m4_v3(shi->obi->mat, v1);
 614
 615                 dface= v1[0]*shi->facenor[0]+v1[1]*shi->facenor[1]+v1[2]*shi->facenor[2];
 616
 617                 /* ortho viewplane cannot intersect using view vector originating in (0,0,0) */
 618                 if(R.r.mode & R_ORTHO) {
 619                         /* x and y 3d coordinate can be derived from pixel coord and winmat */
 620                         float fx= 2.0f/(R.winx*R.winmat[0][0]);
 621                         float fy= 2.0f/(R.winy*R.winmat[1][1]);
 622
 623                         co[0]= (x - 0.5f*R.winx)*fx - R.winmat[3][0]/R.winmat[0][0];
 624                         co[1]= (y - 0.5f*R.winy)*fy - R.winmat[3][1]/R.winmat[1][1];
 625
 626                         /* using a*x + b*y + c*z = d equation, (a b c) is normal */
 627                         if(shi->facenor[2]!=0.0f)
 628                                 co[2]= (dface - shi->facenor[0]*co[0] - shi->facenor[1]*co[1])/shi->facenor[2];
 629                         else
 630                                 co[2]= 0.0f;
 631
 632                         if(dxco && dyco) {
 633                                 dxco[0]= fx;
 634                                 dxco[1]= 0.0f;
 635                                 if(shi->facenor[2]!=0.0f)
 636                                         dxco[2]= -(shi->facenor[0]*fx)/shi->facenor[2];
 637                                 else
 638                                         dxco[2]= 0.0f;
 639
 640                                 dyco[0]= 0.0f;
 641                                 dyco[1]= fy;
 642                                 if(shi->facenor[2]!=0.0f)
 643                                         dyco[2]= -(shi->facenor[1]*fy)/shi->facenor[2];
 644                                 else
 645                                         dyco[2]= 0.0f;
 646
 647                                 if(dxyview) {
 648                                         if(co[2]!=0.0f) fac= 1.0f/co[2]; else fac= 0.0f;
 649                                         dxyview[0]= -R.viewdx*fac;
 650                                         dxyview[1]= -R.viewdy*fac;
 651                                 }
 652                         }
 653                 }
 654                 else {
 655                         float div;
 656
 657                         div= shi->facenor[0]*view[0] + shi->facenor[1]*view[1] + shi->facenor[2]*view[2];
 658                         if (div!=0.0f) fac= dface/div;
 659                         else fac= 0.0f;
 660
 661                         co[0]= fac*view[0];
 662                         co[1]= fac*view[1];
 663                         co[2]= fac*view[2];
 664
 665                         /* pixel dx/dy for render coord */
 666                         if(dxco && dyco) {
 667                                 float u= dface/(div - R.viewdx*shi->facenor[0]);
 668                                 float v= dface/(div - R.viewdy*shi->facenor[1]);
 669
 670                                 dxco[0]= co[0]- (view[0]-R.viewdx)*u;
 671                                 dxco[1]= co[1]- (view[1])*u;
 672                                 dxco[2]= co[2]- (view[2])*u;
 673
 674                                 dyco[0]= co[0]- (view[0])*v;
 675                                 dyco[1]= co[1]- (view[1]-R.viewdy)*v;
 676                                 dyco[2]= co[2]- (view[2])*v;
 677
 678                                 if(dxyview) {
 679                                         if(fac!=0.0f) fac= 1.0f/fac;
 680                                         dxyview[0]= -R.viewdx*fac;
 681                                         dxyview[1]= -R.viewdy*fac;
 682                                 }
 683                         }
 684                 }
 685         }
 686
 687         /* set camera coords - for scanline, it's always 0.0,0.0,0.0 (render is in camera space)
 688          * however for raytrace it can be different - the position of the last intersection */
 689         shi->camera_co[0] = shi->camera_co[1] = shi->camera_co[2] = 0.0f;
 690
 691         /* cannot normalize earlier, code above needs it at viewplane level */
 692         normalize_v3(view);
 693 }
 694
 695 /* from scanline pixel coordinates to 3d coordinates, requires set_triangle */
 696 void shade_input_set_viewco(ShadeInput *shi, float x, float y, float xs, float ys, float z)
 697 {
 698         float *dxyview= NULL, *dxco= NULL, *dyco= NULL;
 699
 700         /* currently in use for dithering (soft shadow), node preview, irregular shad */
 701         shi->xs= (int)xs;
 702         shi->ys= (int)ys;
 703
 704         /* original scanline coordinate without jitter */
 705         shi->scanco[0]= x;
 706         shi->scanco[1]= y;
 707         shi->scanco[2]= z;
 708
 709         /* check if we need derivatives */
 710         if(shi->osatex || (R.r.mode & R_SHADOW)) {
 711                 dxco= shi->dxco;
 712                 dyco= shi->dyco;
 713
 714                 if((shi->mat->texco & TEXCO_REFL))
 715                         dxyview= &shi->dxview;
 716         }
 717
 718         shade_input_calc_viewco(shi, xs, ys, z, shi->view, dxyview, shi->co, dxco, dyco);
 719 }
 720
 721 /* calculate U and V, for scanline (silly render face u and v are in range -1 to 0) */
 722 void shade_input_set_uv(ShadeInput *shi)
 723 {
 724         VlakRen *vlr= shi->vlr;
 725
 726         if((vlr->flag & R_SMOOTH) || (shi->mat->texco & NEED_UV) || (shi->passflag & SCE_PASS_UV)) {
 727                 float v1[3], v2[3], v3[3];
 728
 729                 VECCOPY(v1, shi->v1->co);
 730                 VECCOPY(v2, shi->v2->co);
 731                 VECCOPY(v3, shi->v3->co);
 732
 733                 if(shi->obi->flag & R_TRANSFORMED) {
 734                         mul_m4_v3(shi->obi->mat, v1);
 735                         mul_m4_v3(shi->obi->mat, v2);
 736                         mul_m4_v3(shi->obi->mat, v3);
 737                 }
 738
 739                 /* exception case for wire render of edge */
 740                 if(vlr->v2==vlr->v3) {
 741                         float lend, lenc;
 742
 743                         lend= len_v3v3(v2, v1);
 744                         lenc= len_v3v3(shi->co, v1);
 745
 746                         if(lend==0.0f) {
 747                                 shi->u=shi->v= 0.0f;
 748                         }
 749                         else {
 750                                 shi->u= - (1.0f - lenc/lend);
 751                                 shi->v= 0.0f;
 752                         }
 753
 754                         if(shi->osatex) {
 755                                 shi->dx_u=  0.0f;
 756                                 shi->dx_v=  0.0f;
 757                                 shi->dy_u=  0.0f;
 758                                 shi->dy_v=  0.0f;
 759                         }
 760                 }
 761                 else {
 762                         /* most of this could become re-used for faces */
 763                         float detsh, t00, t10, t01, t11, xn, yn, zn;
 764                         int axis1, axis2;
 765
 766                         /* find most stable axis to project */
 767                         xn= fabs(shi->facenor[0]);
 768                         yn= fabs(shi->facenor[1]);
 769                         zn= fabs(shi->facenor[2]);
 770
 771                         if(zn>=xn && zn>=yn) { axis1= 0; axis2= 1; }
 772                         else if(yn>=xn && yn>=zn) { axis1= 0; axis2= 2; }
 773                         else { axis1= 1; axis2= 2; }
 774
 775                         /* compute u,v and derivatives */
 776                         t00= v3[axis1]-v1[axis1]; t01= v3[axis2]-v1[axis2];
 777                         t10= v3[axis1]-v2[axis1]; t11= v3[axis2]-v2[axis2];
 778
 779                         detsh= (t00*t11-t10*t01);
 780                         detsh= (detsh != 0.0f)? 1.0f/detsh: 0.0f;
 781                         t00*= detsh; t01*=detsh;
 782                         t10*=detsh; t11*=detsh;
 783
 784                         shi->u= (shi->co[axis1]-v3[axis1])*t11-(shi->co[axis2]-v3[axis2])*t10;
 785                         shi->v= (shi->co[axis2]-v3[axis2])*t00-(shi->co[axis1]-v3[axis1])*t01;
 786                         if(shi->osatex) {
 787                                 shi->dx_u=  shi->dxco[axis1]*t11- shi->dxco[axis2]*t10;
 788                                 shi->dx_v=  shi->dxco[axis2]*t00- shi->dxco[axis1]*t01;
 789                                 shi->dy_u=  shi->dyco[axis1]*t11- shi->dyco[axis2]*t10;
 790                                 shi->dy_v=  shi->dyco[axis2]*t00- shi->dyco[axis1]*t01;
 791                         }
 792
 793                         /* u and v are in range -1 to 0, we allow a little bit extra but not too much, screws up speedvectors */
 794                         CLAMP(shi->u, -2.0f, 1.0f);
 795                         CLAMP(shi->v, -2.0f, 1.0f);
 796                 }
 797         }
 798 }
 799
 800 void shade_input_set_normals(ShadeInput *shi)
 801 {
 802         float u= shi->u, v= shi->v;
 803         float l= 1.0f+u+v;
 804
 805         shi->flippednor = 0;
 806
 807         /* test flip normals to viewing direction */
 808         if(!(shi->vlr->flag & R_TANGENT)) {
 809                 if(dot_v3v3(shi->facenor, shi->view) < 0.0f) {
 810                         negate_v3(shi->facenor);
 811                         shi->flippednor= 1;
 812                 }
 813         }
 814
 815         /* calculate vertexnormals */
 816         if(shi->vlr->flag & R_SMOOTH) {
 817                 float *n1= shi->n1, *n2= shi->n2, *n3= shi->n3;
 818
 819                 if(shi->flippednor) {
 820                         negate_v3(n1);
 821                         negate_v3(n2);
 822                         negate_v3(n3);
 823                 }
 824
 825                 shi->vn[0]= l*n3[0]-u*n1[0]-v*n2[0];
 826                 shi->vn[1]= l*n3[1]-u*n1[1]-v*n2[1];
 827                 shi->vn[2]= l*n3[2]-u*n1[2]-v*n2[2];
 828
 829                 // use unnormalized normal (closer to games)
 830                 VECCOPY(shi->nmapnorm, shi->vn);
 831
 832                 normalize_v3(shi->vn);
 833         }
 834         else
 835         {
 836                 VECCOPY(shi->vn, shi->facenor);
 837                 VECCOPY(shi->nmapnorm, shi->vn);
 838         }
 839
 840         /* used in nodes */
 841         VECCOPY(shi->vno, shi->vn);
 842
 843         /* flip normals to viewing direction */
 844         if(!(shi->vlr->flag & R_TANGENT))
 845                 if(dot_v3v3(shi->facenor, shi->view) < 0.0f)
 846                         shade_input_flip_normals(shi);
 847 }
 848
 849 /* XXX shi->flippednor messes up otherwise */
 850 void shade_input_set_vertex_normals(ShadeInput *shi)
 851 {
 852         float u= shi->u, v= shi->v;
 853         float l= 1.0f+u+v;
 854
 855         /* calculate vertexnormals */
 856         if(shi->vlr->flag & R_SMOOTH) {
 857                 float *n1= shi->n1, *n2= shi->n2, *n3= shi->n3;
 858
 859                 shi->vn[0]= l*n3[0]-u*n1[0]-v*n2[0];
 860                 shi->vn[1]= l*n3[1]-u*n1[1]-v*n2[1];
 861                 shi->vn[2]= l*n3[2]-u*n1[2]-v*n2[2];
 862
 863                 // use unnormalized normal (closer to games)
 864                 VECCOPY(shi->nmapnorm, shi->vn);
 865
 866                 normalize_v3(shi->vn);
 867         }
 868         else
 869         {
 870                 VECCOPY(shi->vn, shi->facenor);
 871                 VECCOPY(shi->nmapnorm, shi->vn);
 872         }
 873
 874         /* used in nodes */
 875         VECCOPY(shi->vno, shi->vn);
 876 }
 877
 878
 879 /* use by raytrace, sss, bake to flip into the right direction */
 880 void shade_input_flip_normals(ShadeInput *shi)
 881 {
 882         shi->facenor[0]= -shi->facenor[0];
 883         shi->facenor[1]= -shi->facenor[1];
 884         shi->facenor[2]= -shi->facenor[2];
 885
 886         shi->vn[0]= -shi->vn[0];
 887         shi->vn[1]= -shi->vn[1];
 888         shi->vn[2]= -shi->vn[2];
 889
 890         shi->vno[0]= -shi->vno[0];
 891         shi->vno[1]= -shi->vno[1];
 892         shi->vno[2]= -shi->vno[2];
 893
 894         shi->nmapnorm[0] = -shi->nmapnorm[0];
 895         shi->nmapnorm[1] = -shi->nmapnorm[1];
 896         shi->nmapnorm[2] = -shi->nmapnorm[2];
 897
 898         shi->flippednor= !shi->flippednor;
 899 }
 900
 901 void shade_input_set_shade_texco(ShadeInput *shi)
 902 {
 903         ObjectInstanceRen *obi= shi->obi;
 904         ObjectRen *obr= shi->obr;
 905         VertRen *v1= shi->v1, *v2= shi->v2, *v3= shi->v3;
 906         float u= shi->u, v= shi->v;
 907         float l= 1.0f+u+v, dl;
 908         int mode= shi->mode;            /* or-ed result for all nodes */
 909         short texco= shi->mat->texco;
 910
 911         /* calculate dxno */
 912         if(shi->vlr->flag & R_SMOOTH) {
 913
 914                 if(shi->osatex && (texco & (TEXCO_NORM|TEXCO_REFL)) ) {
 915                         float *n1= shi->n1, *n2= shi->n2, *n3= shi->n3;
 916
 917                         dl= shi->dx_u+shi->dx_v;
 918                         shi->dxno[0]= dl*n3[0]-shi->dx_u*n1[0]-shi->dx_v*n2[0];
 919                         shi->dxno[1]= dl*n3[1]-shi->dx_u*n1[1]-shi->dx_v*n2[1];
 920                         shi->dxno[2]= dl*n3[2]-shi->dx_u*n1[2]-shi->dx_v*n2[2];
 921                         dl= shi->dy_u+shi->dy_v;
 922                         shi->dyno[0]= dl*n3[0]-shi->dy_u*n1[0]-shi->dy_v*n2[0];
 923                         shi->dyno[1]= dl*n3[1]-shi->dy_u*n1[1]-shi->dy_v*n2[1];
 924                         shi->dyno[2]= dl*n3[2]-shi->dy_u*n1[2]-shi->dy_v*n2[2];
 925
 926                 }
 927         }
 928
 929         /* calc tangents */
 930         if (mode & (MA_TANGENT_V|MA_NORMAP_TANG) || R.flag & R_NEED_TANGENT) {
 931                 float *tangent, *s1, *s2, *s3;
 932                 float tl, tu, tv;
 933
 934                 if(shi->vlr->flag & R_SMOOTH) {
 935                         tl= l;
 936                         tu= u;
 937                         tv= v;
 938                 }
 939                 else {
 940                         /* qdn: flat faces have tangents too,
 941                            could pick either one, using average here */
 942                         tl= 1.0f/3.0f;
 943                         tu= -1.0f/3.0f;
 944                         tv= -1.0f/3.0f;
 945                 }
 946
 947                 shi->tang[0]= shi->tang[1]= shi->tang[2]= 0.0f;
 948                 shi->nmaptang[0]= shi->nmaptang[1]= shi->nmaptang[2]= 0.0f;
 949
 950                 if(mode & MA_TANGENT_V) {
 951                         s1 = RE_vertren_get_tangent(obr, v1, 0);
 952                         s2 = RE_vertren_get_tangent(obr, v2, 0);
 953                         s3 = RE_vertren_get_tangent(obr, v3, 0);
 954
 955                         if(s1 && s2 && s3) {
 956                                 shi->tang[0]= (tl*s3[0] - tu*s1[0] - tv*s2[0]);
 957                                 shi->tang[1]= (tl*s3[1] - tu*s1[1] - tv*s2[1]);
 958                                 shi->tang[2]= (tl*s3[2] - tu*s1[2] - tv*s2[2]);
 959
 960                                 if(obi->flag & R_TRANSFORMED)
 961                                         mul_m3_v3(obi->nmat, shi->tang);
 962
 963                                 normalize_v3(shi->tang);
 964                                 VECCOPY(shi->nmaptang, shi->tang);
 965                         }
 966                 }
 967
 968                 if(mode & MA_NORMAP_TANG || R.flag & R_NEED_TANGENT) {
 969                         tangent= RE_vlakren_get_nmap_tangent(obr, shi->vlr, 0);
 970
 971                         if(tangent) {
 972                                 int j1= shi->i1, j2= shi->i2, j3= shi->i3;
 973                                 float c0[3], c1[3], c2[3];
 974
 975                                 vlr_set_uv_indices(shi->vlr, &j1, &j2, &j3);
 976
 977                                 VECCOPY(c0, &tangent[j1*4]);
 978                                 VECCOPY(c1, &tangent[j2*4]);
 979                                 VECCOPY(c2, &tangent[j3*4]);
 980
 981                                 // keeping tangents normalized at vertex level
 982                                 // corresponds better to how it's done in game engines
 983                                 if(obi->flag & R_TRANSFORMED)
 984                                 {
 985                                         mul_mat3_m4_v3(obi->mat, c0); normalize_v3(c0);
 986                                         mul_mat3_m4_v3(obi->mat, c1); normalize_v3(c1);
 987                                         mul_mat3_m4_v3(obi->mat, c2); normalize_v3(c2);
 988                                 }
 989
 990                                 // we don't normalize the interpolated TBN tangent
 991                                 // corresponds better to how it's done in game engines
 992                                 shi->nmaptang[0]= (tl*c2[0] - tu*c0[0] - tv*c1[0]);
 993                                 shi->nmaptang[1]= (tl*c2[1] - tu*c0[1] - tv*c1[1]);
 994                                 shi->nmaptang[2]= (tl*c2[2] - tu*c0[2] - tv*c1[2]);
 995
 996                                 // the sign is the same for all 3 vertices of any
 997                                 // non degenerate triangle.
 998                                 shi->nmaptang[3]= tangent[j1*4+3];
 999                         }
1000                 }
1001         }
1002
1003         if(mode & MA_STR_SURFDIFF) {
1004                 float *surfnor= RE_vlakren_get_surfnor(obr, shi->vlr, 0);
1005
1006                 if(surfnor) {
1007                         VECCOPY(shi->surfnor, surfnor)
1008                         if(obi->flag & R_TRANSFORMED)
1009                                 mul_m3_v3(obi->nmat, shi->surfnor);
1010                 }
1011                 else
1012                         VECCOPY(shi->surfnor, shi->vn)
1013
1014                 shi->surfdist= 0.0f;
1015         }
1016
1017         if(R.r.mode & R_SPEED) {
1018                 float *s1, *s2, *s3;
1019
1020                 s1= RE_vertren_get_winspeed(obi, v1, 0);
1021                 s2= RE_vertren_get_winspeed(obi, v2, 0);
1022                 s3= RE_vertren_get_winspeed(obi, v3, 0);
1023                 if(s1 && s2 && s3) {
1024                         shi->winspeed[0]= (l*s3[0] - u*s1[0] - v*s2[0]);
1025                         shi->winspeed[1]= (l*s3[1] - u*s1[1] - v*s2[1]);
1026                         shi->winspeed[2]= (l*s3[2] - u*s1[2] - v*s2[2]);
1027                         shi->winspeed[3]= (l*s3[3] - u*s1[3] - v*s2[3]);
1028                 }
1029                 else {
1030                         shi->winspeed[0]= shi->winspeed[1]= shi->winspeed[2]= shi->winspeed[3]= 0.0f;
1031                 }
1032         }
1033
1034         /* pass option forces UV calc */
1035         if(shi->passflag & SCE_PASS_UV)
1036                 texco |= (NEED_UV|TEXCO_UV);
1037
1038         /* texture coordinates. shi->dxuv shi->dyuv have been set */
1039         if(texco & NEED_UV) {
1040
1041                 if(texco & TEXCO_ORCO) {
1042                         if(v1->orco) {
1043                                 float *o1, *o2, *o3;
1044
1045                                 o1= v1->orco;
1046                                 o2= v2->orco;
1047                                 o3= v3->orco;
1048
1049                                 shi->lo[0]= l*o3[0]-u*o1[0]-v*o2[0];
1050                                 shi->lo[1]= l*o3[1]-u*o1[1]-v*o2[1];
1051                                 shi->lo[2]= l*o3[2]-u*o1[2]-v*o2[2];
1052
1053                                 if(shi->osatex) {
1054                                         dl= shi->dx_u+shi->dx_v;
1055                                         shi->dxlo[0]= dl*o3[0]-shi->dx_u*o1[0]-shi->dx_v*o2[0];
1056                                         shi->dxlo[1]= dl*o3[1]-shi->dx_u*o1[1]-shi->dx_v*o2[1];
1057                                         shi->dxlo[2]= dl*o3[2]-shi->dx_u*o1[2]-shi->dx_v*o2[2];
1058                                         dl= shi->dy_u+shi->dy_v;
1059                                         shi->dylo[0]= dl*o3[0]-shi->dy_u*o1[0]-shi->dy_v*o2[0];
1060                                         shi->dylo[1]= dl*o3[1]-shi->dy_u*o1[1]-shi->dy_v*o2[1];
1061                                         shi->dylo[2]= dl*o3[2]-shi->dy_u*o1[2]-shi->dy_v*o2[2];
1062                                 }
1063                         }
1064
1065                         VECCOPY(shi->duplilo, obi->dupliorco);
1066                 }
1067
1068                 if(texco & TEXCO_GLOB) {
1069                         VECCOPY(shi->gl, shi->co);
1070                         mul_m4_v3(R.viewinv, shi->gl);
1071                         if(shi->osatex) {
1072                                 VECCOPY(shi->dxgl, shi->dxco);
1073                                 mul_mat3_m4_v3(R.viewinv, shi->dxgl);
1074                                 VECCOPY(shi->dygl, shi->dyco);
1075                                 mul_mat3_m4_v3(R.viewinv, shi->dygl);
1076                         }
1077                 }
1078
1079                 if(texco & TEXCO_STRAND) {
1080                         shi->strandco= (l*v3->accum - u*v1->accum - v*v2->accum);
1081                         if(shi->osatex) {
1082                                 dl= shi->dx_u+shi->dx_v;
1083                                 shi->dxstrand= dl*v3->accum-shi->dx_u*v1->accum-shi->dx_v*v2->accum;
1084                                 dl= shi->dy_u+shi->dy_v;
1085                                 shi->dystrand= dl*v3->accum-shi->dy_u*v1->accum-shi->dy_v*v2->accum;
1086                         }
1087                 }
1088
1089                 if((texco & TEXCO_UV) || (mode & (MA_VERTEXCOL|MA_VERTEXCOLP|MA_FACETEXTURE)))  {
1090                         VlakRen *vlr= shi->vlr;
1091                         MTFace *tface;
1092                         MCol *mcol;
1093                         char *name;
1094                         int i, j1=shi->i1, j2=shi->i2, j3=shi->i3;
1095
1096                         /* uv and vcols are not copied on split, so set them according vlr divide flag */
1097                         vlr_set_uv_indices(vlr, &j1, &j2, &j3);
1098
1099                         shi->totuv= 0;
1100                         shi->totcol= 0;
1101                         shi->actuv= obr->actmtface;
1102                         shi->actcol= obr->actmcol;
1103
1104                         if(mode & (MA_VERTEXCOL|MA_VERTEXCOLP)) {
1105                                 for (i=0; (mcol=RE_vlakren_get_mcol(obr, vlr, i, &name, 0)); i++) {
1106                                         ShadeInputCol *scol= &shi->col[i];
1107                                         char *cp1, *cp2, *cp3;
1108
1109                                         shi->totcol++;
1110                                         scol->name= name;
1111
1112                                         cp1= (char *)(mcol+j1);
1113                                         cp2= (char *)(mcol+j2);
1114                                         cp3= (char *)(mcol+j3);
1115
1116                                         /* alpha value */
1117                                         scol->col[3]= (l*((float)cp3[0]) - u*((float)cp1[0]) - v*((float)cp2[0]))/255.0f;
1118
1119                                         /* try to prevent invalid color sampling of zero alpha points */
1120                                         if (!cp1[0]) cp1 = cp2; if (!cp1[0]) cp1 = cp3;
1121                                         if (!cp2[0]) cp2 = cp1; if (!cp2[0]) cp2 = cp3;
1122                                         if (!cp3[0]) cp3 = cp1; if (!cp3[0]) cp3 = cp2;
1123
1124                                         /* sample color value */
1125                                         scol->col[0]= (l*((float)cp3[3]) - u*((float)cp1[3]) - v*((float)cp2[3]))/255.0f;
1126                                         scol->col[1]= (l*((float)cp3[2]) - u*((float)cp1[2]) - v*((float)cp2[2]))/255.0f;
1127                                         scol->col[2]= (l*((float)cp3[1]) - u*((float)cp1[1]) - v*((float)cp2[1]))/255.0f;
1128                                 }
1129
1130                                 if(shi->totcol) {
1131                                         shi->vcol[0]= shi->col[shi->actcol].col[0];
1132                                         shi->vcol[1]= shi->col[shi->actcol].col[1];
1133                                         shi->vcol[2]= shi->col[shi->actcol].col[2];
1134                                         shi->vcol[3]= shi->col[shi->actcol].col[3];
1135                                 }
1136                                 else {
1137                                         shi->vcol[0]= 0.0f;
1138                                         shi->vcol[1]= 0.0f;
1139                                         shi->vcol[2]= 0.0f;
1140                                         shi->vcol[3]= 1.0f;
1141                                 }
1142                         }
1143
1144                         for (i=0; (tface=RE_vlakren_get_tface(obr, vlr, i, &name, 0)); i++) {
1145                                 ShadeInputUV *suv= &shi->uv[i];
1146                                 float *uv1, *uv2, *uv3;
1147
1148                                 shi->totuv++;
1149                                 suv->name= name;
1150
1151                                 uv1= tface->uv[j1];
1152                                 uv2= tface->uv[j2];
1153                                 uv3= tface->uv[j3];
1154
1155                                 suv->uv[0]= -1.0f + 2.0f*(l*uv3[0]-u*uv1[0]-v*uv2[0]);
1156                                 suv->uv[1]= -1.0f + 2.0f*(l*uv3[1]-u*uv1[1]-v*uv2[1]);
1157                                 suv->uv[2]= 0.0f;       /* texture.c assumes there are 3 coords */
1158
1159                                 if(shi->osatex) {
1160                                         float duv[2];
1161
1162                                         dl= shi->dx_u+shi->dx_v;
1163                                         duv[0]= shi->dx_u;
1164                                         duv[1]= shi->dx_v;
1165
1166                                         suv->dxuv[0]= 2.0f*(dl*uv3[0]-duv[0]*uv1[0]-duv[1]*uv2[0]);
1167                                         suv->dxuv[1]= 2.0f*(dl*uv3[1]-duv[0]*uv1[1]-duv[1]*uv2[1]);
1168
1169                                         dl= shi->dy_u+shi->dy_v;
1170                                         duv[0]= shi->dy_u;
1171                                         duv[1]= shi->dy_v;
1172
1173                                         suv->dyuv[0]= 2.0f*(dl*uv3[0]-duv[0]*uv1[0]-duv[1]*uv2[0]);
1174                                         suv->dyuv[1]= 2.0f*(dl*uv3[1]-duv[0]*uv1[1]-duv[1]*uv2[1]);
1175                                 }
1176
1177                                 if((mode & MA_FACETEXTURE) && i==obr->actmtface) {
1178                                         if((mode & (MA_VERTEXCOL|MA_VERTEXCOLP))==0) {
1179                                                 shi->vcol[0]= 1.0f;
1180                                                 shi->vcol[1]= 1.0f;
1181                                                 shi->vcol[2]= 1.0f;
1182                                                 shi->vcol[3]= 1.0f;
1183                                         }
1184                                         if(tface && tface->tpage)
1185                                                 render_realtime_texture(shi, tface->tpage);
1186                                 }
1187
1188
1189                         }
1190
1191                         shi->dupliuv[0]= -1.0f + 2.0f*obi->dupliuv[0];
1192                         shi->dupliuv[1]= -1.0f + 2.0f*obi->dupliuv[1];
1193                         shi->dupliuv[2]= 0.0f;
1194
1195                         if(shi->totuv == 0) {
1196                                 ShadeInputUV *suv= &shi->uv[0];
1197
1198                                 suv->uv[0]= 2.0f*(u+.5f);
1199                                 suv->uv[1]= 2.0f*(v+.5f);
1200                                 suv->uv[2]= 0.0f;       /* texture.c assumes there are 3 coords */
1201
1202                                 if(mode & MA_FACETEXTURE) {
1203                                         /* no tface? set at 1.0f */
1204                                         shi->vcol[0]= 1.0f;
1205                                         shi->vcol[1]= 1.0f;
1206                                         shi->vcol[2]= 1.0f;
1207                                         shi->vcol[3]= 1.0f;
1208                                 }
1209                         }
1210                 }
1211
1212                 if(texco & TEXCO_NORM) {
1213                         shi->orn[0]= -shi->vn[0];
1214                         shi->orn[1]= -shi->vn[1];
1215                         shi->orn[2]= -shi->vn[2];
1216                 }
1217
1218                 if(texco & TEXCO_STRESS) {
1219                         float *s1, *s2, *s3;
1220
1221                         s1= RE_vertren_get_stress(obr, v1, 0);
1222                         s2= RE_vertren_get_stress(obr, v2, 0);
1223                         s3= RE_vertren_get_stress(obr, v3, 0);
1224                         if(s1 && s2 && s3) {
1225                                 shi->stress= l*s3[0] - u*s1[0] - v*s2[0];
1226                                 if(shi->stress<1.0f) shi->stress-= 1.0f;
1227                                 else shi->stress= (shi->stress-1.0f)/shi->stress;
1228                         }
1229                         else shi->stress= 0.0f;
1230                 }
1231
1232                 if(texco & TEXCO_TANGENT) {
1233                         if((mode & MA_TANGENT_V)==0) {
1234                                 /* just prevent surprises */
1235                                 shi->tang[0]= shi->tang[1]= shi->tang[2]= 0.0f;
1236                                 shi->nmaptang[0]= shi->nmaptang[1]= shi->nmaptang[2]= 0.0f;
1237                         }
1238                 }
1239         }
1240
1241         /* this only avalailable for scanline renders */
1242         if(shi->depth==0) {
1243                 float x= shi->xs;
1244                 float y= shi->ys;
1245
1246                 if(texco & TEXCO_WINDOW) {
1247                         shi->winco[0]= -1.0f + 2.0f*x/(float)R.winx;
1248                         shi->winco[1]= -1.0f + 2.0f*y/(float)R.winy;
1249                         shi->winco[2]= 0.0f;
1250                         if(shi->osatex) {
1251                                 shi->dxwin[0]= 2.0f/(float)R.winx;
1252                                 shi->dywin[1]= 2.0f/(float)R.winy;
1253                                 shi->dxwin[1]= shi->dxwin[2]= 0.0f;
1254                                 shi->dywin[0]= shi->dywin[2]= 0.0f;
1255                         }
1256                 }
1257
1258                 if(texco & TEXCO_STICKY) {
1259                         float *s1, *s2, *s3;
1260
1261                         s1= RE_vertren_get_sticky(obr, v1, 0);
1262                         s2= RE_vertren_get_sticky(obr, v2, 0);
1263                         s3= RE_vertren_get_sticky(obr, v3, 0);
1264
1265                         if(s1 && s2 && s3) {
1266                                 float obwinmat[4][4], winmat[4][4], ho1[4], ho2[4], ho3[4];
1267                                 float Zmulx, Zmuly;
1268                                 float hox, hoy, l, dl, u, v;
1269                                 float s00, s01, s10, s11, detsh;
1270
1271                                 /* old globals, localized now */
1272                                 Zmulx=  ((float)R.winx)/2.0f; Zmuly=  ((float)R.winy)/2.0f;
1273
1274                                 zbuf_make_winmat(&R, winmat);
1275                                 if(shi->obi->flag & R_TRANSFORMED)
1276                                         mul_m4_m4m4(obwinmat, obi->mat, winmat);
1277                                 else
1278                                         copy_m4_m4(obwinmat, winmat);
1279
1280                                 zbuf_render_project(obwinmat, v1->co, ho1);
1281                                 zbuf_render_project(obwinmat, v2->co, ho2);
1282                                 zbuf_render_project(obwinmat, v3->co, ho3);
1283
1284                                 s00= ho3[0]/ho3[3] - ho1[0]/ho1[3];
1285                                 s01= ho3[1]/ho3[3] - ho1[1]/ho1[3];
1286                                 s10= ho3[0]/ho3[3] - ho2[0]/ho2[3];
1287                                 s11= ho3[1]/ho3[3] - ho2[1]/ho2[3];
1288
1289                                 detsh= s00*s11-s10*s01;
1290                                 detsh= (detsh != 0.0f)? 1.0f/detsh: 0.0f;
1291                                 s00*= detsh; s01*=detsh;
1292                                 s10*=detsh; s11*=detsh;
1293
1294                                 /* recalc u and v again */
1295                                 hox= x/Zmulx -1.0f;
1296                                 hoy= y/Zmuly -1.0f;
1297                                 u= (hox - ho3[0]/ho3[3])*s11 - (hoy - ho3[1]/ho3[3])*s10;
1298                                 v= (hoy - ho3[1]/ho3[3])*s00 - (hox - ho3[0]/ho3[3])*s01;
1299                                 l= 1.0f+u+v;
1300
1301                                 shi->sticky[0]= l*s3[0]-u*s1[0]-v*s2[0];
1302                                 shi->sticky[1]= l*s3[1]-u*s1[1]-v*s2[1];
1303                                 shi->sticky[2]= 0.0f;
1304
1305                                 if(shi->osatex) {
1306                                         float dxuv[2], dyuv[2];
1307                                         dxuv[0]=  s11/Zmulx;
1308                                         dxuv[1]=  - s01/Zmulx;
1309                                         dyuv[0]=  - s10/Zmuly;
1310                                         dyuv[1]=  s00/Zmuly;
1311
1312                                         dl= dxuv[0] + dxuv[1];
1313                                         shi->dxsticky[0]= dl*s3[0] - dxuv[0]*s1[0] - dxuv[1]*s2[0];
1314                                         shi->dxsticky[1]= dl*s3[1] - dxuv[0]*s1[1] - dxuv[1]*s2[1];
1315                                         dl= dyuv[0] + dyuv[1];
1316                                         shi->dysticky[0]= dl*s3[0] - dyuv[0]*s1[0] - dyuv[1]*s2[0];
1317                                         shi->dysticky[1]= dl*s3[1] - dyuv[0]*s1[1] - dyuv[1]*s2[1];
1318                                 }
1319                         }
1320                 }
1321         } /* else {
1322          Note! For raytracing winco is not set, important because thus means all shader input's need to have their variables set to zero else in-initialized values are used
1323         */
1324         if (shi->do_manage) {
1325                 if(mode & (MA_VERTEXCOL|MA_VERTEXCOLP|MA_FACETEXTURE)) {
1326                         srgb_to_linearrgb_v3_v3(shi->vcol, shi->vcol);
1327                 }
1328         }
1329
1330 }
1331
1332 /* ****************** ShadeSample ************************************** */
1333
1334 /* initialize per part, not per pixel! */
1335 void shade_input_initialize(ShadeInput *shi, RenderPart *pa, RenderLayer *rl, int sample)
1336 {
1337
1338         memset(shi, 0, sizeof(ShadeInput));
1339
1340         shi->sample= sample;
1341         shi->thread= pa->thread;
1342         shi->do_preview= (R.r.scemode & R_MATNODE_PREVIEW) != 0;
1343         shi->do_manage= (R.r.color_mgt_flag & R_COLOR_MANAGEMENT);
1344         shi->lay= rl->lay;
1345         shi->layflag= rl->layflag;
1346         shi->passflag= rl->passflag;
1347         shi->combinedflag= ~rl->pass_xor;
1348         shi->mat_override= rl->mat_override;
1349         shi->light_override= rl->light_override;
1350 //      shi->rl= rl;
1351         /* note shi.depth==0  means first hit, not raytracing */
1352
1353 }
1354
1355 /* initialize per part, not per pixel! */
1356 void shade_sample_initialize(ShadeSample *ssamp, RenderPart *pa, RenderLayer *rl)
1357 {
1358         int a, tot;
1359
1360         tot= R.osa==0?1:R.osa;
1361
1362         for(a=0; a<tot; a++) {
1363                 shade_input_initialize(&ssamp->shi[a], pa, rl, a);
1364                 memset(&ssamp->shr[a], 0, sizeof(ShadeResult));
1365         }
1366
1367         get_sample_layers(pa, rl, ssamp->rlpp);
1368 }
1369
1370 /* Do AO or (future) GI */
1371 void shade_samples_do_AO(ShadeSample *ssamp)
1372 {
1373         ShadeInput *shi;
1374         int sample;
1375
1376         if(!(R.r.mode & R_SHADOW))
1377                 return;
1378         if(!(R.r.mode & R_RAYTRACE) && !(R.wrld.ao_gather_method == WO_AOGATHER_APPROX))
1379                 return;
1380
1381         if(R.wrld.mode & (WO_AMB_OCC|WO_ENV_LIGHT|WO_INDIRECT_LIGHT)) {
1382                 shi= &ssamp->shi[0];
1383
1384                 if(((shi->passflag & SCE_PASS_COMBINED) && (shi->combinedflag & (SCE_PASS_AO|SCE_PASS_ENVIRONMENT|SCE_PASS_INDIRECT)))
1385                         || (shi->passflag & (SCE_PASS_AO|SCE_PASS_ENVIRONMENT|SCE_PASS_INDIRECT)))
1386                         for(sample=0, shi= ssamp->shi; sample<ssamp->tot; shi++, sample++)
1387                                 if(!(shi->mode & MA_SHLESS))
1388                                         ambient_occlusion(shi);         /* stores in shi->ao[] */
1389         }
1390 }
1391
1392
1393 void shade_samples_fill_with_ps(ShadeSample *ssamp, PixStr *ps, int x, int y)
1394 {
1395         ShadeInput *shi;
1396         float xs, ys;
1397
1398         ssamp->tot= 0;
1399
1400         for(shi= ssamp->shi; ps; ps= ps->next) {
1401                 shade_input_set_triangle(shi, ps->obi, ps->facenr, 1);
1402
1403                 if(shi->vlr) {  /* NULL happens for env material or for 'all z' */
1404                         unsigned short curmask= ps->mask;
1405
1406                         /* full osa is only set for OSA renders */
1407                         if(shi->vlr->flag & R_FULL_OSA) {
1408                                 short shi_cp= 0, samp;
1409
1410                                 for(samp=0; samp<R.osa; samp++) {
1411                                         if(curmask & (1<<samp)) {
1412                                                 /* zbuffer has this inverse corrected, ensures xs,ys are inside pixel */
1413                                                 xs= (float)x + R.jit[samp][0] + 0.5f;
1414                                                 ys= (float)y + R.jit[samp][1] + 0.5f;
1415
1416                                                 if(shi_cp)
1417                                                         shade_input_copy_triangle(shi, shi-1);
1418
1419                                                 shi->mask= (1<<samp);
1420 //                                              shi->rl= ssamp->rlpp[samp];
1421                                                 shi->samplenr= R.shadowsamplenr[shi->thread]++; /* this counter is not being reset per pixel */
1422                                                 shade_input_set_viewco(shi, x, y, xs, ys, (float)ps->z);
1423                                                 shade_input_set_uv(shi);
1424                                                 if(shi_cp==0)
1425                                                         shade_input_set_normals(shi);
1426                                                 else  /* XXX shi->flippednor messes up otherwise */
1427                                                         shade_input_set_vertex_normals(shi);
1428
1429                                                 shi_cp= 1;
1430                                                 shi++;
1431                                         }
1432                                 }
1433                         }
1434                         else {
1435                                 if(R.osa) {
1436                                         short b= R.samples->centmask[curmask];
1437                                         xs= (float)x + R.samples->centLut[b & 15] + 0.5f;
1438                                         ys= (float)y + R.samples->centLut[b>>4] + 0.5f;
1439                                 }
1440                                 else {
1441                                         xs= (float)x + 0.5f;
1442                                         ys= (float)y + 0.5f;
1443                                 }
1444
1445                                 shi->mask= curmask;
1446                                 shi->samplenr= R.shadowsamplenr[shi->thread]++;
1447                                 shade_input_set_viewco(shi, x, y, xs, ys, (float)ps->z);
1448                                 shade_input_set_uv(shi);
1449                                 shade_input_set_normals(shi);
1450                                 shi++;
1451                         }
1452
1453                         /* total sample amount, shi->sample is static set in initialize */
1454                         if(shi!=ssamp->shi)
1455                                 ssamp->tot= (shi-1)->sample + 1;
1456                 }
1457         }
1458 }
1459
1460 /* shades samples, returns true if anything happened */
1461 int shade_samples(ShadeSample *ssamp, PixStr *ps, int x, int y)
1462 {
1463         shade_samples_fill_with_ps(ssamp, ps, x, y);
1464
1465         if(ssamp->tot) {
1466                 ShadeInput *shi= ssamp->shi;
1467                 ShadeResult *shr= ssamp->shr;
1468                 int samp;
1469
1470                 /* if shadow or AO? */
1471                 shade_samples_do_AO(ssamp);
1472
1473                 /* if shade (all shadepinputs have same passflag) */
1474                 if(ssamp->shi[0].passflag & ~(SCE_PASS_Z|SCE_PASS_INDEXOB|SCE_PASS_INDEXMA)) {
1475
1476                         for(samp=0; samp<ssamp->tot; samp++, shi++, shr++) {
1477                                 shade_input_set_shade_texco(shi);
1478                                 shade_input_do_shade(shi, shr);
1479                         }
1480                 }
1481                 else if(shi->passflag & SCE_PASS_Z) {
1482                         for(samp=0; samp<ssamp->tot; samp++, shi++, shr++)
1483                                 shr->z= -shi->co[2];
1484                 }
1485
1486                 return 1;
1487         }
1488         return 0;
1489 }
1490