Perl变量作用域和自定义函数小结

想MS用perl编程，先熟悉perl的语法，没必要深究，看http://perldoc.perl.org/perlintro.html就够了，主要熟悉变量、operator，file and I/O。 然后就是到MS看功能函数或模块，这才是重点。

1. 变量定义问题

Perl默认所有变量为包变量（Package variables），包变量为全局变量，这意味着程序的任何其他部分，甚至在其他文件里定义的子程序，都能影响和修改变量的值。在一定程度上讲，这样是“不安全”的。

 

my变量：自Perl 5以后，增加了新的非全局变量，也称为词法变量、私有变量、局部变量，或称为my变量。注意，不是local变量，local变量有另外的含义。例如my $a;就定义了一个局部变量$a，它的作用域是当前块（Block），通俗地讲就是大括号里面，如果没有大括号就是从定义的地方开始到程序结束。而local变量的意思，是

 

local变量：local和my明显的不同，my创建局部变量，而local作用于包变量。具体地讲，local $x实际做的事是：存储包变量$x的当前值在一个安全的地方，然后用一个新值替换它，假如没有指定新值，就使用undef代 替。当控制离开当前块时，它也会恢复$x的旧值。它影响的是包变量，这个包变量获取了本地值。但包变量总是全局的，local申明的包变量亦无例外。为了 显示其区别，请看这个：

        $lo = 'global'; 
    $m  = 'global'; 
    A();

    sub A { 
      local $lo = 'AAA'; 
      my    $m  = 'AAA'; 
      B(); 
    }

    sub B { 
      print "B ", ($lo eq 'AAA' ? 'can' : 'cannot') , 
            " see the value of lo set by A.\n";

      print "B ", ($m  eq 'AAA' ? 'can' : 'cannot') , 
            " see the value of m  set by A.\n"; 
        }

结果会打印：

    B can see the value of lo set by A. 
    B cannot see the value of m  set by A.

 

发生了什么？在A函数里的local申明，给包变量$lo赋予了一个新的临时值AAA。旧值global会被存储起来，直到A返回，但在这点之前，A调用了B。B访问$lo的内容没有问题，因为$lo是包变量，包变量总是全局可见的，所以它能见到A设置的AAA值。

 

2. 自定义函数

 

用户函数又称子程序(Subroutine)，在Perl中用下面的结构来定义用户函数：

　　　　sub 子程序名{

　　　　　　语句块;
　　　　}

这里的子程序名与变量的取名规则类似。

以显示欢迎词的程序为例：

　　　　sub say_hello{

　　　　　　print "你好，欢迎光临网上学园";
　　　　}

　　用户函数的定义可以位于程序的任何位置，比如说放在文件的未尾。如果两个子程序使用了相同的程序名， 后面的子程序将覆盖前面子程序。

　　用户函数中的变量默认为全局变量，与其他程序共享。

　　用户函数的调用：通过在子程序前加“&”调用，可在任一表达式内调用。 子程序中可以再调用另外的子程序。

　　调用用户函数产生的结果称为返回值(return value)。返回值是每次调用函数中最后一个表达式的计算值。 以加法函数为例：

　　　　sub add_a_b{

　　　　　　$a+$b;
　　　　}

　　函数最后一条表达式为$a+$b，故返回值为$a+$b。以下是调用情况：

　　　　$a=5;
　　　　$b=6;
　　　　$c=&add_a_b; #$c的值为11
　　　　$d=5*&add_a_b; #$d的值为5*11即55

　　上述的函数功能与传统直接写在程序中没什么两样，如果加上参数传递就可以实现全新的功能了。 在Perl中，如果函数调用后面跟着一个用括号括起来的列表，则在函数调用期间该列表将被自动分配给以@_命名的特殊变量。 函数可以访问该变量，从而确定参数的个数及赋值。

　　仍以加法函数为例：

　　　　sub add_a_b{
　　　　　　$_〔0〕+$_〔1〕;
　　　　}
　　　　$c=&add_a_b(5,6); #$c的值为11
　　　　$d=5*&add_a_b(2,3); #d的值为5*5即25

　　如何改变参数的个数呢？我们可以用循环的方式来实现：

　　　　sub add_all{
　　　　　　$sum=0; #将sum初始化
　　　　　　foreach $_(@_) { #遍历参数列表
　　　　　　　$sum+=$_; #累加每个元素
　　　　　　}
　　　　　　$sum; #返回sum即总和的值
　　　　}
　　　　$a=&add_all(3,4,5); #$a的值为3+4+5即12
　　　　$d=2*&add_all(1,2,3,4,5); #d的值为2*15即30

　　既然函数中的变量全为全程变量，那么上述程序中若调用程序中含有$sum变量时将替换，这不是我们所要的。 那么如何解决这一问题呢？

　　答案就是：使用局部变量, 使用local()操作符就可实现此功能。在上面的程序中，只需在第一行$sum=0;前加入：

　　　　local($sum);
　　当函数执行时，$sum的全程变量的值被保留起来，同时建立一个局部变量$sum。

 

用perl进行截断能收敛性测试（Convergence test）的例子，以供大家参考。

 

===================================================================

#!perl

use strict;
use MaterialsScript qw(:all);

# 定义计算文档和结果输出文档

my $myDoc=$Documents{"WO3.xsd"};#要计算的模型文件
my $myStudyTable=Documents->new("energy-encut.std");#新建StudyTable存放计算结果，第一列是截断能，第二列是体系总能
my $mySheet=$myStudyTable->ActiveSheet;
$mySheet->ColumnHeading(0)="Enegy Cutoff(eV)";
$mySheet->ColumnHeading(1)="Final Enegy(eV)";

# castep single calculation
my $castep=Modules->CASTEP;
my $startEnergy=200;#测试起点:200 eV
my $endEnergy=500;#测试终点: 500 eV
my $intervalEnergy=20;#测试点间隔
my $sumIteration=($endEnergy-$startEnergy)/$intervalEnergy;

#循环计算每个截断能测试点的单点能
for(my $counter=0;$counter<=$sumIteration;++$counter){
        my $energyCutoff=$startEnergy+$intervalEnergy*$counter;

        $castep->ChangeSettings(
                Settings(Quality=>"Fine",
                UseCustomEnergyCutoff=>"Yes",
                EnergyCutoff=>$energyCutoff));#这里设置截断能为测试点的截断能
        $castep->Energy->Run($myDoc);
        
        $mySheet->Cell($counter,0)=$energyCutoff;
        
        #read final energy from castep output files
        #下面这一段提取castep文件中的总能，即Final energy后面的结果
        foreach my $line (@{$Documents{"WO3.castep"}->Lines}) {
    if ($line=~/^Final energy/){
            my $finalEnergy = substr($line,31,15);
            print $energyCutoff,"\t",$finalEnergy,"\n";
            $mySheet->Cell($counter,1)=$finalEnergy;#将总能结果放入StudyTable
    }}
}

本文引用地址： http://bbs.sciencenet.cn/home.php?mod=space&uid=548899&do=blog&id=471288