2626// DEALINGS IN THE SOFTWARE.
2727//
2828
29- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
29+
3030#include < teem/air.h>
3131#include < teem/ten.h>
32- # endif
32+
3333
3434#include < Core/Parser/ArrayMathFunctionCatalog.h>
3535#include < Core/Math/MiscMath.h>
@@ -290,9 +290,9 @@ bool eigvec1_t(SCIRun::ArrayMathProgramCode& pc)
290290 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
291291 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
292292 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
293- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
293+
294294 tenEigensolve_d (eval, evec, ten);
295- # endif
295+
296296 data0[0 ] = evec[0 ];
297297 data0[1 ] = evec[1 ];
298298 data0[2 ] = evec[2 ];
@@ -317,9 +317,9 @@ bool eigvec2_t(SCIRun::ArrayMathProgramCode& pc)
317317 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
318318 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
319319 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
320- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
320+
321321 tenEigensolve_d (eval, evec, ten);
322- # endif
322+
323323 data0[0 ] = evec[3 ];
324324 data0[1 ] = evec[4 ];
325325 data0[2 ] = evec[5 ];
@@ -344,9 +344,9 @@ bool eigvec3_t(SCIRun::ArrayMathProgramCode& pc)
344344 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
345345 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
346346 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
347- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
347+
348348 tenEigensolve_d (eval, evec, ten);
349- # endif
349+
350350
351351 data0[0 ] = evec[6 ];
352352 data0[1 ] = evec[7 ];
@@ -372,9 +372,9 @@ bool eigval1_t(SCIRun::ArrayMathProgramCode& pc)
372372 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
373373 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
374374 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
375- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
375+
376376 tenEigensolve_d (eval, evec, ten);
377- # endif
377+
378378
379379 data0[0 ] = eval[0 ];
380380
@@ -398,9 +398,9 @@ bool eigval2_t(SCIRun::ArrayMathProgramCode& pc)
398398 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
399399 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
400400 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
401- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
401+
402402 tenEigensolve_d (eval, evec, ten);
403- # endif
403+
404404
405405 data0[0 ] = eval[1 ];
406406
@@ -424,9 +424,9 @@ bool eigval3_t(SCIRun::ArrayMathProgramCode& pc)
424424 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
425425 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
426426 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
427- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
427+
428428 tenEigensolve_d (eval, evec, ten);
429- # endif
429+
430430
431431 data0[0 ] = eval[2 ];
432432
@@ -467,9 +467,9 @@ bool trace_t(SCIRun::ArrayMathProgramCode& pc)
467467 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
468468 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
469469
470- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
470+
471471 tenEigensolve_d (eval, evec, ten);
472- # endif
472+
473473
474474 data0[0 ] = eval[0 ]+eval[1 ]+eval[2 ];
475475
@@ -493,9 +493,9 @@ bool det_t(SCIRun::ArrayMathProgramCode& pc)
493493 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
494494 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
495495 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
496- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
496+
497497 tenEigensolve_d (eval, evec, ten);
498- # endif
498+
499499
500500 data0[0 ] = eval[0 ]*eval[1 ]*eval[2 ];
501501
@@ -520,9 +520,9 @@ bool B_t(SCIRun::ArrayMathProgramCode& pc)
520520 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
521521 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
522522 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
523- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
523+
524524 tenEigensolve_d (eval, evec, ten);
525- # endif
525+
526526
527527 data0[0 ] = eval[0 ]*eval[1 ]+eval[0 ]*eval[2 ]+eval[2 ]*eval[1 ];
528528
@@ -546,9 +546,9 @@ bool S_t(SCIRun::ArrayMathProgramCode& pc)
546546 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
547547 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
548548 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
549- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
549+
550550 tenEigensolve_d (eval, evec, ten);
551- # endif
551+
552552
553553 data0[0 ] = eval[0 ]*eval[0 ]+eval[1 ]*eval[1 ]+eval[2 ]*eval[2 ];
554554
@@ -572,9 +572,9 @@ bool quality_t(SCIRun::ArrayMathProgramCode& pc)
572572 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
573573 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
574574 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
575- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
575+
576576 tenEigensolve_d (eval, evec, ten);
577- # endif
577+
578578
579579 data0[0 ] = ((eval[0 ]*eval[0 ]+eval[1 ]*eval[1 ]+eval[2 ]*eval[2 ]) -
580580 (eval[0 ]*eval[1 ]+eval[0 ]*eval[2 ]+eval[2 ]*eval[1 ]))/9.0 ;
@@ -599,9 +599,9 @@ bool frobenius_t(SCIRun::ArrayMathProgramCode& pc)
599599 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
600600 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
601601 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
602- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
602+
603603 tenEigensolve_d (eval, evec, ten);
604- # endif
604+
605605
606606 data0[0 ] = ::sqrt (eval[0 ]*eval[0 ]+eval[1 ]*eval[1 ]+eval[2 ]*eval[2 ]);
607607
@@ -625,9 +625,9 @@ bool frobenius2_t(SCIRun::ArrayMathProgramCode& pc)
625625 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
626626 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
627627 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
628- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
628+
629629 tenEigensolve_d (eval, evec, ten);
630- # endif
630+
631631
632632 data0[0 ] = eval[0 ]*eval[0 ]+eval[1 ]*eval[1 ]+eval[2 ]*eval[2 ];
633633
@@ -652,9 +652,9 @@ bool fracanisotropy_t(SCIRun::ArrayMathProgramCode& pc)
652652 ten[0 ] = 1.0 ; ten[1 ] = data1[0 ]; ten[2 ] = data1[1 ];
653653 ten[3 ] = data1[2 ]; ten[4 ] = data1[3 ];
654654 ten[5 ] = data1[4 ]; ten[6 ] = data1[5 ];
655- # ifdef SCIRUN4_CODE_TO_BE_ENABLED_LATER
655+
656656 tenEigensolve_d (eval, evec, ten);
657- # endif
657+
658658
659659 double S = eval[0 ]*eval[0 ]+eval[1 ]*eval[1 ]+eval[2 ]*eval[2 ];
660660 double B = eval[0 ]*eval[1 ]+eval[0 ]*eval[2 ]+eval[2 ]*eval[1 ];
0 commit comments