Beefy Boxes and Bandwidth Generously Provided by pair Networks
There's more than one way to do things
 
PerlMonks  

Re^5: Parsing Guassian '03 Log Files

by GrandFather (Saint)
on Jan 31, 2008 at 20:25 UTC ( [id://665433]=note: print w/replies, xml ) Need Help??


in reply to Re^4: Parsing Guassian '03 Log Files
in thread Parsing Guassian '03 Log Files

You can clean that up a little:

#!/usr/bin/perl use strict; use warnings; # Parse Gaussian '03 output files for the lengths of the bonds in # Diaminopolymethine Dyes strictly between Carbon and Nitrogen or Carb +on and # Carbon chomp (my $infile = <>); chomp (my $outfile = <>); open FILE, '<', $infile or die "Unable to open $infile: $!"; my @inlog = grep {/\d/ && ! /[a-zA-Z]/ && $_ >= 1} map {split (/\s/, $_)} grep {!/h|c{3}?|nc{2}?|(estimate)|c{2}?n{1}?/} grep {/^\s?!{1}?\s*(c|n)/} <FILE>; close FILE; open LOG, '>', $outfile or die "Unable to create $outfile: $!"; print LOG join "\n", @inlog; close LOG;

Your version had some rather odd constructs. The pop loops to clear out arrays were perhaps the strangest. Much better to:

@array = ();

The while loop to slurp the file is better as:

@array = <FILE>;

The construct my $infile, my $outfile; is odd. Either my ($var1, $var1); or use two separate statements:

my $var1; my $var2;

However you should declare your variables as close to their first use as possible so it is not often that you need to declare a bunch of variables in one place like that anyway. See too the file name variable declarations in my version of your code above.

You should always use the three parameter version of open and you should always check the result.


Perl is environmentally friendly - it saves trees

Replies are listed 'Best First'.
Re^6: Parsing Guassian '03 Log Files
by Andrew_Levenson (Hermit) on Jan 31, 2008 at 21:31 UTC
    Ooh, @array = ();? Nice. Yeah, I completely forgot how to empty an array, so I... made my own way. I guess we know now why I'm a chemistry major now instead of Computer Science, haha.

    If its okay, i'm going to test your script and replace mine in the bin? Or not, I'll figure it out when I get back to the lab.
    Again, thanks.
    C(qw/74 97 104 112/);sub C{while(@_){$c**=$C;print (map{chr($C!=$c?shift:pop)}$_),$C+=@_%2!=1?1:0}}

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://665433]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others pondering the Monastery: (4)
As of 2024-04-19 04:21 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found