[SOLVED] Use of python to organize and assemble .comm files

rene_epfl

Hello,

I have several analyses to make that have some parts of the .comm file in common. So I would like to write things that are common to each .comm file in seperate files and then call these separate files from the .comm file.

The idea is that if I make a change in lets say the definition of a list, of a material, a calculation stage etc. I don't have to make the change in every .comm file manually, but only in one python script that can be called when used.

In the folder containing the python script I have defined an empty file called __init__.py as well as the files I want to call (e.g myproject_materials.py, calc_step0.py, etc.)

I have tried something like this:

# defining the path of the project
import sys
sys.path.append('/home/code_aster_work/python')

# or should I use: "sys.path.insert(0,'/home/code_aster_work/python')"?

# here I define some materials
import myproject_materials

(...)

# here I define the calculation of step 0
import calc_step_0

But it doesn't work..

the error message that I get is the following:

ERREUR A L'INTERPRETATION DANS ACCAS - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR phase d'initialisation
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! erreur non prevue et non traitee prevenir la maintenance fort.1            !
   ! Traceback (most recent call last):                                         !
   !    File "/tmp/carpaij-LMSPC55-interactif.6364-LMSPC55/proc.0/Python/Noyau/ !
   ! N_JDC.py", line 209, in exec_compile                                       !
   !     exec self.proc_compile in self.g_context                               !
   !    File "fort.1", line 9, in <module>                                      !
   !     import calc_stage_0.py                                                 !
   !  ImportError: No module named calc_stage_0.py                              !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR phase d'initialisation

does anyone know what to do or do you have a functioning example?

Thank you,

René

reiteres

Hello Rene,

normally this should work. I have for example a Python module for Unit conversion (see here https://github.com/maldun/UnitCalculator)

and it is imported normally like this:

import sys


import os

# Assuming module is in home directory
sys.path.append(os.path.expanduser("~"))

from UnitCalculator import *

But be aware that Code_Aster is not directly Python. Especially Aster Concepts are not real Python objects but more something homebrewn from EDF to communicate with the Fortran Code underneath. So for example if you add some Python function into your code, you can't return a concept.

But there are Macros. You can write your own macros in Python for Code_Aster which can then called as normal Aster commands, where you also can return aster Concepts correctly. But those macros have their own rules to follow.

See this Training Slides for making your own macro: http://code-aster.org/V2/UPLOAD/DOC/Formations/04a_Macro-command.pdf

I have recently started to write macros myself. If you get the gist of it it is quite easy and rewarding.

hope this answers your question.

rene_epfl

Hello Reiteres,

Thank you for your very useful answer and example code!

I will start working on it and when it will (hopefully) work like i want it to, i will post what i have done.

have a nice day,

René

rene_epfl

Hello,

It's partially working :-). I now have a .comm file in which i can define parameters for my calculation and then open an python-script containing information that is also valid for other calculations, for which there is another .comm file.

However variables defined in the .comm file are not recognized during the execution of the python-scripts.

I did the following:

# this is the .comm file, with par_lot='non' to enable execution line by line
DEBUT(	PAR_LOT='NON',)

# add the folder with my .py-scripts
import sys
sys.path.append('/home/my_project/python_scripts')

# here i define a variable named pc
pc=100

(...)

import my_python_script

(...)
FIN();

And the file my_python_script.py looks something like this:

# importing some functions of code aster
from Cata.cata import FORMULE, DEFI_CONSTANTE
(...)
# here i use (in whatever way) the variable pc defined in the .comm file
pres=DEFI_CONSTANTE(VALE=pc)
# or
if pc==3:
    do this
else
    do something else
(...)

The error I get is the following:

ERREUR A L'INTERPRETATION DANS ACCAS - INTERRUPTION >> JDC.py : DEBUT RAPPORT CR phase d'initialisation !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! erreur de syntaxe, Erreur de nom: name 'pc' is not defined ligne 132 ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

fin CR phase d'initialisation

So it seems that the variable 'pc' is not known although it was previously defined. Do i specifically have to import the already defined variables into the python-script? I tried to find a way how to do that by looking through some .py files in the /bibpyt folder of the aster installation but I didn't find a solution. Does someone know how to get access to variables in a .py file? is there some thing like:

import 'list of all variables that were defined'?

Thank you

René

reiteres

Hi!

Be aware that when importing a Python file as module it is loaded as a class and not that the script simply gets executed!

Additionally this is not proper programming to simply use variables somewhere defined (especially in a language where everything could be simply overwritten)

The proper way to do this is to define functions/classes in your Python script file, then import them and provide them your variable as input.

Then the variable is awaited as input.

Example:

#my_python.py


def func(pc):

do_something_with_pc

and in the .comm file

from my_python import func

func(pc)

be aware that Python functions may take python objects and aster concepts, but they don't return aster concepts only python objects are handled correctly!

Yet there is a (unclean) way to achieve what you want to do: the exec statement/function from Python ( see e.g. http://late.am/post/2012/04/30/the-exec-statement-and-a-python-mystery.html )

exec[/b] simple executes a string of code. This means you could import your my_python_script.py as a textfile for example with [b]open[/b] and then convert the whole text to a string and provide it to the [b]exec statement.

Since exec[/b] works like you would magically write the code segment at the position where [b]exec is placed, also variables defined before would be known.

But be careful whan going such a way!

jlf

Hi,

Did you try the INCLUDE command ?

JLF

reiteres

Nice!

Didn't know that!

rene_epfl

Hi reiteres and jlf,

Thank you for your help and inputs!

I've started to write a macro command but I realized that right now I don't have enough time to learn how to do that. I eventually wrote a short python script, that does more or less what i want, to work around the problem. Maybe I'll have some more time during summer to work on this project, if so I'll try to write macro commands for my problems and post them.

In case you're interested in what i finally did or if you want to do something similar, here is the python script I wrote:

# This python scripts assembles different parts of a command file that


# are written as python scripts.

import os
import fileinput

# define the filename and number
file_name='initial_cond_test_'
file_number=46

# define which files you want to put together (in this order):
if file_number==51:
	input_files=[	'a_beginning.py',
			'b_mesh_input.py',
			'c_model_definition.py',
			'd_material_properties_version_3.py',
			'e_material_definition_case_1.py',
			'f_material_fields.py',
			'g_heat_flow_version_2.py',
			'h_time_list.py',
			'i_initial_conditions_case_1.py',
			'j_boundary_conditions_case_1.py',
			'k_calc_step_0.py',
			'k_calc_step_1.py',
			'k_calc_step_2.py',
			'k_calc_step_3.py',
			'k_calc_step_4_version_1.py',
			'y_fin.py',]

elif file_number==52:
	input_files=[	'y_poursuite.py',
			'k_calc_step_5.py',
			'k_calc_step_6.py',
			'y_fin.py',]

elif file_number==43:
	input_files=[	'y_poursuite.py',
			'n_post_processing_fields.py',
			'y_fin.py',]

elif file_number==46:
	input_files=[	'y_poursuite.py',
			'm_post_processing_tables.py',
			'y_fin.py',]

# define the path containing the python scripts with the model-information
input_path='/home/my_name/my_project/3d/python_scripts/'

# define the path where the new command file is to be saved at
output_path='/home/my_name/my_project/3d/comm_files/'

# assembling the output file name
output_file_name=file_name+str(output_file_number)+'.comm'

# adding the path to the filenames
n=len(input_files)
for i in range(0,n):
    input_files[i]=input_path+input_files[i]

# creating a new .comm file where the files listed above are inserted
with open(os.path.join(output_path,output_file_name),'w') as output_file:
	input_lines=fileinput.input(input_files)

output_file.writelines(input_lines)

have a nice day!

René

DanielOliveira

If you want to import all statements of a python script (this include variables) you should use a command like this:

from myscript import *

This will import variables.