[ripudaman]

Hi,

I am relatively inexperienced in pyFoam. I am using bits and pieces of information that I have gathered from the various pyFoam talks and presentations over the years.

I am having a problem in being able to run a case, generate a specific log file and read the log file for a regExp to modify the output on the screen all in one go. I am able to do all this but I am having to run the solver twice. I want to find a way to do this without having to run the solver twice. Also another tips on this will be very useful.

Here is the snippet of code that I am using

run=AnalyzedRunner(Analyzer(),silent=True,argv=["convergeFracWidthFoam","-case",work.name],logname='CFWF')
run.start()

CONVERGED=LogAnalyzer("CONVERGED")
analyzer=LogAnalyzerApplication(CONVERGED)
analyzer.run(path.join(target,"CFWF.logfile"))

Here is the logAnalyzer.py file that I use to parse the regular expressions in my created log file

import re,sys
from PyFoam.LogAnalysis.FoamLogAnalyzer import FoamLogAnalyzer
from PyFoam.LogAnalysis.LogAnalyzerApplication import LogAnalyzerApplication
from PyFoam.LogAnalysis.LogLineAnalyzer import LogLineAnalyzer
from PyFoam.LogAnalysis.BoundingLogAnalyzer import BoundingLogAnalyzer
class LogAnalyzer(FoamLogAnalyzer):
    def __init__(self,name):
    super(LogAnalyzer,self).__init__()
    self.addAnalyzer(name,LineAnalyzer())
    def obtainWidth(self,name):
    return self.getAnalyzer(name).obtainWidth()
    def obtainHaltRatio(self,name):
    return self.getAnalyzer(name).obtainHaltRatio()

class LineAnalyzer(LogLineAnalyzer):
    def __init__(self):
        LogLineAnalyzer.__init__(self)
 
        self.told=""
        self.exp1=re.compile("^(.+):  Solving for (.+), Initial residual = (.+), Final residual = (.+), No Iterations (.+)$")
     self.exp2=re.compile("Previous Fracture Width = (.+)  Current Fracture Width = (.+)$")
    self.exp3=re.compile("Current Halt Ratio = (.+)  Desired Halt Ratio = (.+)$")
    self.currentwidth=0
    self.haltratio=1

    def doAnalysis(self,line):
        m=self.exp1.match(line)
        if m!=None:
            name=m.group(2)
            resid=m.group(3)
        iters=m.group(5)
            time=self.getTime()
            if time!=self.told:
                self.told=time
                print "\n t = %3g: " % ( float(time) ),
            print "%3s: %3e %5s: %1g" % (name,float(resid),"NIter",float(iters)),
            #sys.stdout.flush()

     m=self.exp2.match(line)
    if m!=None:
        previous=m.group(1)
        current=m.group(2)
        print "\t %5s: %6e" % ("FracWidth",float(current)),
        self.currentwidth=float(current)
        sys.stdout.flush()

     m=self.exp3.match(line)
    if m!=None:
        current=m.group(1)
        target=m.group(2)
        self.haltratio=float(current)

    def obtainWidth(self):
    return self.currentwidth

    def obtainHaltRatio(self):
    return self.haltratio

class Analyzer(BoundingLogAnalyzer):
    def __init__(self):
        BoundingLogAnalyzer.__init__(self)
        self.addAnalyzer("Analyzer",LineAnalyzer())

I apologize for any redundancy here since I do not fully understand the way pyFoam libraries are set and am not very familiar with the class system in python.

Any help here is appreciated.

Thank you.
Regards,
Ripu