doc/user/images/medrenumber.py

   1 # -*- coding: utf-8 -*-
   2
   3 ###
   4 ### This script is intended to be launched in a new SALOME study
   5 ###
   6
   7 import sys
   8 import salome
   9
  10 salome.salome_init()
  11 theStudy = salome.myStudy
  12
  13 import salome_notebook
  14 notebook = salome_notebook.NoteBook(theStudy)
  15 sys.path.insert( 0, r'/misc/dn27/users_Linux/eap/salome/tmp')
  16
  17 import iparameters
  18 ipar = iparameters.IParameters(salome.myStudy.GetCommonParameters("Interface Applicative", 1), True)
  19
  20 #Set up visual properties:
  21 ipar.setProperty("AP_ACTIVE_VIEW", "VTKViewer_0_0")
  22 ipar.setProperty("AP_WORKSTACK_INFO", "0000000100000000000000020100000001000003a0000000040000000200000002000000080000001a004f00430043005600690065007700650072005f0030005f00300000000102000000080000001a00560054004b005600690065007700650072005f0030005f00300000000202")
  23 ipar.setProperty("AP_ACTIVE_MODULE", "Mesh")
  24 ipar.setProperty("AP_SAVEPOINT_NAME", "GUI state: 1")
  25 #Set up lists:
  26 # fill list AP_VIEWERS_LIST
  27 ipar.append("AP_VIEWERS_LIST", "VTKViewer_2")
  28 # fill list VTKViewer_2
  29 ipar.append("VTKViewer_2", "VTK scene:1 - viewer:1")
  30 ipar.append("VTKViewer_2", """<?xml version="1.0"?>
  31 <ViewState>
  32     <Position X="50.1781" Y="3.66147" Z="1055.58"/>
  33     <FocalPoint X="50.1781" Y="-2.45923" Z="0.0142599"/>
  34     <ViewUp X="0" Y="0.999983" Z="-0.0057984"/>
  35     <ViewScale Parallel="189.306" X="1" Y="1" Z="1"/>
  36     <DisplayCubeAxis Show="0"/>
  37     <GraduatedAxis Axis="X">
  38         <Title isVisible="1" Text="X" Font="0" Bold="0" Italic="0" Shadow="0">
  39             <Color R="1" G="0" B="0"/>
  40         </Title>
  41         <Labels isVisible="1" Number="3" Offset="2" Font="0" Bold="0" Italic="0" Shadow="0">
  42             <Color R="1" G="0" B="0"/>
  43         </Labels>
  44         <TickMarks isVisible="1" Length="5"/>
  45     </GraduatedAxis>
  46     <GraduatedAxis Axis="Y">
  47         <Title isVisible="1" Text="Y" Font="0" Bold="0" Italic="0" Shadow="0">
  48             <Color R="0" G="1" B="0"/>
  49         </Title>
  50         <Labels isVisible="1" Number="3" Offset="2" Font="0" Bold="0" Italic="0" Shadow="0">
  51             <Color R="0" G="1" B="0"/>
  52         </Labels>
  53         <TickMarks isVisible="1" Length="5"/>
  54     </GraduatedAxis>
  55     <GraduatedAxis Axis="Z">
  56         <Title isVisible="1" Text="Z" Font="0" Bold="0" Italic="0" Shadow="0">
  57             <Color R="0" G="0" B="1"/>
  58         </Title>
  59         <Labels isVisible="1" Number="3" Offset="2" Font="0" Bold="0" Italic="0" Shadow="0">
  60             <Color R="0" G="0" B="1"/>
  61         </Labels>
  62         <TickMarks isVisible="1" Length="5"/>
  63     </GraduatedAxis>
  64     <Trihedron isShown="0" Size="100"/>
  65     <Background Value="bt=1;fn=;tm=0;ts=false;c1=#ffffff;c2=#000000;gt=-1;gr="/>
  66 </ViewState>
  67 """)
  68 # fill list AP_MODULES_LIST
  69 ipar.append("AP_MODULES_LIST", "Mesh")
  70
  71
  72 ###
  73 ### GEOM component
  74 ###
  75
  76 import GEOM
  77 from salome.geom import geomBuilder
  78 import math
  79 import SALOMEDS
  80
  81
  82 geompy = geomBuilder.New(theStudy)
  83
  84 Face_1 = geompy.MakeFaceHW(100, 100, 1)
  85 geompy.addToStudy( Face_1, 'Face_1' )
  86
  87 ###
  88 ### SMESH component
  89 ###
  90
  91 import  SMESH, SALOMEDS
  92 from salome.smesh import smeshBuilder
  93
  94 smesh = smeshBuilder.New(theStudy)
  95 Mesh_1 = smesh.Mesh(Face_1)
  96 Regular_1D = Mesh_1.Segment()
  97 Number_of_Segments_1 = Regular_1D.NumberOfSegments(3)
  98 MEFISTO_2D = Mesh_1.Triangle(algo=smeshBuilder.MEFISTO)
  99 isDone = Mesh_1.Compute()
 100 isDone = Mesh_1.RemoveElements( range( 1,13 ))
 101 Mesh_1.ExportMED( r'mesh1.med', overwrite=1 )
 102
 103
 104 #from MEDCoupling import *
 105 from MEDLoader import ReadMeshFromFile, WriteMesh
 106 m = ReadMeshFromFile("mesh1.med")
 107
 108 from MEDRenumber import RenumberingFactory
 109 ren=RenumberingFactory("BOOST")
 110 a,b=m.computeNeighborsOfCells()
 111 n2o,_=ren.renumber(a,b)
 112 mrenum=m[n2o]
 113 WriteMesh("mesh2.med",mrenum,True)
 114
 115 ([Mesh_renum], status) = smesh.CreateMeshesFromMED(r'mesh2.med')
 116 Mesh_renum.TranslateObject( Mesh_renum, [ 120, 0, 0 ], 0 )
 117
 118
 119 ### Store presentation parameters of displayed objects
 120 import iparameters
 121 ipar = iparameters.IParameters(theStudy.GetModuleParameters("Interface Applicative", "SMESH", 1))
 122
 123 #Set up entries:
 124 # set up entry SMESH_3 (Mesh_1) parameters
 125 ipar.setParameter("SMESH_3", "VTKViewer_0_Visibility", "On")
 126 ipar.setParameter("SMESH_3", "VTKViewer_0_Representation", "2")
 127 ipar.setParameter("SMESH_3", "VTKViewer_0_IsShrunk", "0")
 128 ipar.setParameter("SMESH_3", "VTKViewer_0_Entities", "e:0:f:1:v:0:0d:0:b:0")
 129 ipar.setParameter("SMESH_3", "VTKViewer_0_Colors", "surface:0:0.666667:1:backsurface:100:volume:1:0:0.666667:-100:edge:0:0.666667:1:node:1:0:0:outline:0:0.27451:0:elem0d:0:1:0:ball:0:0.333333:1:orientation:1:1:1")
 130 ipar.setParameter("SMESH_3", "VTKViewer_0_Sizes", "line:1:outline:1:elem0d:5:ball:10:1:shrink:0.75:orientation:0.1:0")
 131 ipar.setParameter("SMESH_3", "VTKViewer_0_PointMarker", "std:1:9")
 132 ipar.setParameter("SMESH_3", "VTKViewer_0_Opacity", "1")
 133 ipar.setParameter("SMESH_3", "VTKViewer_0_ClippingPlane", "Off")
 134 # set up entry SMESH_4 (Mesh_Renum) parameters
 135 ipar.setParameter("SMESH_4", "VTKViewer_0_Visibility", "On")
 136 ipar.setParameter("SMESH_4", "VTKViewer_0_Representation", "2")
 137 ipar.setParameter("SMESH_4", "VTKViewer_0_IsShrunk", "0")
 138 ipar.setParameter("SMESH_4", "VTKViewer_0_Entities", "e:1:f:1:v:1:0d:1:b:1")
 139 ipar.setParameter("SMESH_4", "VTKViewer_0_Colors", "surface:0:0.666667:1:backsurface:100:volume:1:0:0.666667:-100:edge:0:0.666667:1:node:1:0:0:outline:0:0.27451:0:elem0d:0:1:0:ball:0:0.333333:1:orientation:1:1:1")
 140 ipar.setParameter("SMESH_4", "VTKViewer_0_Sizes", "line:1:outline:1:elem0d:5:ball:10:1:shrink:0.75:orientation:0.1:0")
 141 ipar.setParameter("SMESH_4", "VTKViewer_0_PointMarker", "std:1:9")
 142 ipar.setParameter("SMESH_4", "VTKViewer_0_Opacity", "1")
 143 ipar.setParameter("SMESH_4", "VTKViewer_0_ClippingPlane", "Off")
 144
 145
 146 if salome.sg.hasDesktop():
 147   salome.sg.updateObjBrowser(True)
 148   iparameters.getSession().restoreVisualState(1)