347 |
|
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . "_???.flt"; |
348 |
|
print "{\n"; |
349 |
|
system "$cmd" || die "Failure running rcalc"; |
350 |
+ |
for (my $i = ($tensortree==3)*$ns*$ns*$ns/2; $i-- > 0; ) { |
351 |
+ |
print "0\n"; |
352 |
+ |
} |
353 |
|
print "}\n"; |
354 |
|
} |
355 |
|
print |
383 |
|
$cmd .= "$td/" . ($fmodnm,$bmodnm)[$forw] . "_???.flt"; |
384 |
|
print "{\n"; |
385 |
|
system "$cmd" || die "Failure running rcalc"; |
386 |
+ |
for (my $i = ($tensortree==3)*$ns*$ns*$ns/2; $i-- > 0; ) { |
387 |
+ |
print "0\n"; |
388 |
+ |
} |
389 |
|
print "}\n"; |
390 |
|
} |
391 |
|
print |
404 |
|
$tcal = ' |
405 |
|
DEGREE : PI/180; |
406 |
|
sq(x) : x*x; |
407 |
< |
Kpola(r) : select(r+1, -5, 5, 15, 25, 35, 45, 55, 65, 75, 90); |
407 |
> |
Kpola(r) : select(r+1, 0, 5, 15, 25, 35, 45, 55, 65, 75, 90); |
408 |
|
Knaz(r) : select(r, 1, 8, 16, 20, 24, 24, 24, 16, 12); |
409 |
|
Kaccum(r) : if(r-.5, Knaz(r) + Kaccum(r-1), 0); |
410 |
|
Kmax : Kaccum(Knaz(0)); |
497 |
|
<AngleBasis> |
498 |
|
<AngleBasisName>LBNL/Klems Full</AngleBasisName> |
499 |
|
<AngleBasisBlock> |
500 |
< |
<Theta>0</Theta> |
501 |
< |
<nPhis>1</nPhis> |
502 |
< |
<ThetaBounds> |
503 |
< |
<LowerTheta>0</LowerTheta> |
504 |
< |
<UpperTheta>5</UpperTheta> |
505 |
< |
</ThetaBounds> |
506 |
< |
</AngleBasisBlock> |
507 |
< |
<AngleBasisBlock> |
508 |
< |
<Theta>10</Theta> |
509 |
< |
<nPhis>8</nPhis> |
510 |
< |
<ThetaBounds> |
511 |
< |
<LowerTheta>5</LowerTheta> |
512 |
< |
<UpperTheta>15</UpperTheta> |
513 |
< |
</ThetaBounds> |
514 |
< |
</AngleBasisBlock> |
515 |
< |
<AngleBasisBlock> |
516 |
< |
<Theta>20</Theta> |
517 |
< |
<nPhis>16</nPhis> |
518 |
< |
<ThetaBounds> |
519 |
< |
<LowerTheta>15</LowerTheta> |
520 |
< |
<UpperTheta>25</UpperTheta> |
521 |
< |
</ThetaBounds> |
522 |
< |
</AngleBasisBlock> |
523 |
< |
<AngleBasisBlock> |
524 |
< |
<Theta>30</Theta> |
525 |
< |
<nPhis>20</nPhis> |
526 |
< |
<ThetaBounds> |
527 |
< |
<LowerTheta>25</LowerTheta> |
528 |
< |
<UpperTheta>35</UpperTheta> |
529 |
< |
</ThetaBounds> |
530 |
< |
</AngleBasisBlock> |
531 |
< |
<AngleBasisBlock> |
532 |
< |
<Theta>40</Theta> |
533 |
< |
<nPhis>24</nPhis> |
534 |
< |
<ThetaBounds> |
535 |
< |
<LowerTheta>35</LowerTheta> |
536 |
< |
<UpperTheta>45</UpperTheta> |
537 |
< |
</ThetaBounds> |
538 |
< |
</AngleBasisBlock> |
539 |
< |
<AngleBasisBlock> |
540 |
< |
<Theta>50</Theta> |
541 |
< |
<nPhis>24</nPhis> |
542 |
< |
<ThetaBounds> |
543 |
< |
<LowerTheta>45</LowerTheta> |
544 |
< |
<UpperTheta>55</UpperTheta> |
545 |
< |
</ThetaBounds> |
546 |
< |
</AngleBasisBlock> |
547 |
< |
<AngleBasisBlock> |
548 |
< |
<Theta>60</Theta> |
549 |
< |
<nPhis>24</nPhis> |
550 |
< |
<ThetaBounds> |
551 |
< |
<LowerTheta>55</LowerTheta> |
552 |
< |
<UpperTheta>65</UpperTheta> |
553 |
< |
</ThetaBounds> |
554 |
< |
</AngleBasisBlock> |
555 |
< |
<AngleBasisBlock> |
556 |
< |
<Theta>70</Theta> |
557 |
< |
<nPhis>16</nPhis> |
558 |
< |
<ThetaBounds> |
559 |
< |
<LowerTheta>65</LowerTheta> |
560 |
< |
<UpperTheta>75</UpperTheta> |
561 |
< |
</ThetaBounds> |
562 |
< |
</AngleBasisBlock> |
563 |
< |
<AngleBasisBlock> |
564 |
< |
<Theta>82.5</Theta> |
565 |
< |
<nPhis>12</nPhis> |
566 |
< |
<ThetaBounds> |
567 |
< |
<LowerTheta>75</LowerTheta> |
568 |
< |
<UpperTheta>90</UpperTheta> |
569 |
< |
</ThetaBounds> |
500 |
> |
<Theta>0</Theta> |
501 |
> |
<nPhis>1</nPhis> |
502 |
> |
<ThetaBounds> |
503 |
> |
<LowerTheta>0</LowerTheta> |
504 |
> |
<UpperTheta>5</UpperTheta> |
505 |
> |
</ThetaBounds> |
506 |
> |
</AngleBasisBlock> |
507 |
> |
<AngleBasisBlock> |
508 |
> |
<Theta>10</Theta> |
509 |
> |
<nPhis>8</nPhis> |
510 |
> |
<ThetaBounds> |
511 |
> |
<LowerTheta>5</LowerTheta> |
512 |
> |
<UpperTheta>15</UpperTheta> |
513 |
> |
</ThetaBounds> |
514 |
> |
</AngleBasisBlock> |
515 |
> |
<AngleBasisBlock> |
516 |
> |
<Theta>20</Theta> |
517 |
> |
<nPhis>16</nPhis> |
518 |
> |
<ThetaBounds> |
519 |
> |
<LowerTheta>15</LowerTheta> |
520 |
> |
<UpperTheta>25</UpperTheta> |
521 |
> |
</ThetaBounds> |
522 |
> |
</AngleBasisBlock> |
523 |
> |
<AngleBasisBlock> |
524 |
> |
<Theta>30</Theta> |
525 |
> |
<nPhis>20</nPhis> |
526 |
> |
<ThetaBounds> |
527 |
> |
<LowerTheta>25</LowerTheta> |
528 |
> |
<UpperTheta>35</UpperTheta> |
529 |
> |
</ThetaBounds> |
530 |
> |
</AngleBasisBlock> |
531 |
> |
<AngleBasisBlock> |
532 |
> |
<Theta>40</Theta> |
533 |
> |
<nPhis>24</nPhis> |
534 |
> |
<ThetaBounds> |
535 |
> |
<LowerTheta>35</LowerTheta> |
536 |
> |
<UpperTheta>45</UpperTheta> |
537 |
> |
</ThetaBounds> |
538 |
> |
</AngleBasisBlock> |
539 |
> |
<AngleBasisBlock> |
540 |
> |
<Theta>50</Theta> |
541 |
> |
<nPhis>24</nPhis> |
542 |
> |
<ThetaBounds> |
543 |
> |
<LowerTheta>45</LowerTheta> |
544 |
> |
<UpperTheta>55</UpperTheta> |
545 |
> |
</ThetaBounds> |
546 |
> |
</AngleBasisBlock> |
547 |
> |
<AngleBasisBlock> |
548 |
> |
<Theta>60</Theta> |
549 |
> |
<nPhis>24</nPhis> |
550 |
> |
<ThetaBounds> |
551 |
> |
<LowerTheta>55</LowerTheta> |
552 |
> |
<UpperTheta>65</UpperTheta> |
553 |
> |
</ThetaBounds> |
554 |
> |
</AngleBasisBlock> |
555 |
> |
<AngleBasisBlock> |
556 |
> |
<Theta>70</Theta> |
557 |
> |
<nPhis>16</nPhis> |
558 |
> |
<ThetaBounds> |
559 |
> |
<LowerTheta>65</LowerTheta> |
560 |
> |
<UpperTheta>75</UpperTheta> |
561 |
> |
</ThetaBounds> |
562 |
> |
</AngleBasisBlock> |
563 |
> |
<AngleBasisBlock> |
564 |
> |
<Theta>82.5</Theta> |
565 |
> |
<nPhis>12</nPhis> |
566 |
> |
<ThetaBounds> |
567 |
> |
<LowerTheta>75</LowerTheta> |
568 |
> |
<UpperTheta>90</UpperTheta> |
569 |
> |
</ThetaBounds> |
570 |
|
</AngleBasisBlock> |
571 |
|
</AngleBasis> |
572 |
|
</DataDefinition> |