doc/user/images/intersect_2d2d.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 import math
  10
  11 salome.salome_init()
  12 theStudy = salome.myStudy
  13
  14 import iparameters
  15 ipar = iparameters.IParameters(salome.myStudy.GetCommonParameters("Interface Applicative", 1), True)
  16
  17 #Set up visual properties:
  18 ipar.setProperty("AP_ACTIVE_VIEW", "VTKViewer_0_0")
  19 ipar.setProperty("AP_WORKSTACK_INFO", "0000000100000000000000020100000001000003a0000000040000000100000000000000080000001a00560054004b005600690065007700650072005f0030005f00300000000002")
  20 ipar.setProperty("AP_ACTIVE_MODULE", "Mesh")
  21 ipar.setProperty("AP_SAVEPOINT_NAME", "GUI state: 2")
  22 #Set up lists:
  23 # fill list AP_VIEWERS_LIST
  24 ipar.append("AP_VIEWERS_LIST", "VTKViewer_1")
  25 # fill list VTKViewer_1
  26 ipar.append("VTKViewer_1", "VTK scene:1 - viewer:1")
  27 ipar.append("VTKViewer_1", """<?xml version="1.0"?>
  28 <ViewState>
  29     <Position X="5.78221" Y="-0.0156473" Z="86.6339"/>
  30     <FocalPoint X="5.78221" Y="-0.0156473" Z="0"/>
  31     <ViewUp X="0" Y="1" Z="0"/>
  32     <ViewScale Parallel="8.81461" X="1" Y="1" Z="1"/>
  33     <DisplayCubeAxis Show="0"/>
  34     <GraduatedAxis Axis="X">
  35         <Title isVisible="1" Text="X" Font="0" Bold="0" Italic="0" Shadow="0">
  36             <Color R="1" G="0" B="0"/>
  37         </Title>
  38         <Labels isVisible="1" Number="3" Offset="2" Font="0" Bold="0" Italic="0" Shadow="0">
  39             <Color R="1" G="0" B="0"/>
  40         </Labels>
  41         <TickMarks isVisible="1" Length="5"/>
  42     </GraduatedAxis>
  43     <GraduatedAxis Axis="Y">
  44         <Title isVisible="1" Text="Y" Font="0" Bold="0" Italic="0" Shadow="0">
  45             <Color R="0" G="1" B="0"/>
  46         </Title>
  47         <Labels isVisible="1" Number="3" Offset="2" Font="0" Bold="0" Italic="0" Shadow="0">
  48             <Color R="0" G="1" B="0"/>
  49         </Labels>
  50         <TickMarks isVisible="1" Length="5"/>
  51     </GraduatedAxis>
  52     <GraduatedAxis Axis="Z">
  53         <Title isVisible="1" Text="Z" Font="0" Bold="0" Italic="0" Shadow="0">
  54             <Color R="0" G="0" B="1"/>
  55         </Title>
  56         <Labels isVisible="1" Number="3" Offset="2" Font="0" Bold="0" Italic="0" Shadow="0">
  57             <Color R="0" G="0" B="1"/>
  58         </Labels>
  59         <TickMarks isVisible="1" Length="5"/>
  60     </GraduatedAxis>
  61     <Trihedron isShown="0" Size="100"/>
  62     <Background Value="bt=1;fn=;tm=0;ts=false;c1=#ffffff;c2=#000000;gt=-1;gr="/>
  63 </ViewState>
  64 """)
  65 # fill list AP_MODULES_LIST
  66 ipar.append("AP_MODULES_LIST", "Mesh")
  67
  68 if sys.platform == "win32":
  69     from MEDCouplingCompat import *
  70 else:
  71     from MEDCoupling import *
  72
  73 coordsArr=DataArrayDouble(range(6))
  74 mesh2d=MEDCouplingCMesh("mesh2")
  75 mesh2d.setCoords(coordsArr,coordsArr[:2])
  76 mesh2d=mesh2d.buildUnstructured()
  77
  78 mesh1d=mesh2d.deepCopy()
  79 mesh1d.rotate( [2.3,0], math.radians( 25 ))
  80 #mesh1d.translate( [0.2,0.4] )
  81 mesh1d.setName("mesh1")
  82
  83 m2d,a2d,a1d=MEDCouplingUMesh.Intersect2DMeshes( mesh2d, mesh1d, 1e-12 )
  84
  85 from MEDLoader import WriteMesh
  86 WriteMesh("mesh1.med",mesh2d,True)
  87 WriteMesh("mesh1.med",mesh1d,False)
  88 WriteMesh("mesh1.med",m2d,False)
  89
  90 ###
  91 ### SMESH component
  92 ###
  93
  94 import  SMESH, SALOMEDS
  95 from salome.smesh import smeshBuilder
  96
  97 smesh = smeshBuilder.New(theStudy)
  98 ([m_2d, mesh_1d, mesh_2d], status) = smesh.CreateMeshesFromMED(r'mesh1.med')
  99
 100 m_2d.TranslateObject( m_2d, [6,0,0], False )
 101
 102
 103 ### Store presentation parameters of displayed objects
 104 import iparameters
 105 ipar = iparameters.IParameters(theStudy.GetModuleParameters("Interface Applicative", "SMESH", 1))
 106
 107 #Set up entries:
 108 # set up entry SMESH_4 (mesh1) parameters
 109 ipar.setParameter("SMESH_4", "VTKViewer_0_Visibility", "On")
 110 ipar.setParameter("SMESH_4", "VTKViewer_0_Representation", "1")
 111 ipar.setParameter("SMESH_4", "VTKViewer_0_IsShrunk", "0")
 112 ipar.setParameter("SMESH_4", "VTKViewer_0_Entities", "e:0:f:1:v:0:0d:0:b:0")
 113 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")
 114 ipar.setParameter("SMESH_4", "VTKViewer_0_Sizes", "line:1:outline:1:elem0d:5:ball:10:1:shrink:0.75:orientation:0.1:0")
 115 ipar.setParameter("SMESH_4", "VTKViewer_0_PointMarker", "std:1:9")
 116 ipar.setParameter("SMESH_4", "VTKViewer_0_Opacity", "1")
 117 ipar.setParameter("SMESH_4", "VTKViewer_0_ClippingPlane", "Off")
 118 # set up entry SMESH_3 (Intersect2D) parameters
 119 ipar.setParameter("SMESH_3", "VTKViewer_0_Visibility", "On")
 120 ipar.setParameter("SMESH_3", "VTKViewer_0_Representation", "1")
 121 ipar.setParameter("SMESH_3", "VTKViewer_0_IsShrunk", "0")
 122 ipar.setParameter("SMESH_3", "VTKViewer_0_Entities", "e:0:f:1:v:0:0d:0:b:0")
 123 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")
 124 ipar.setParameter("SMESH_3", "VTKViewer_0_Sizes", "line:1:outline:1:elem0d:5:ball:10:1:shrink:0.75:orientation:0.1:0")
 125 ipar.setParameter("SMESH_3", "VTKViewer_0_PointMarker", "std:1:9")
 126 ipar.setParameter("SMESH_3", "VTKViewer_0_Opacity", "1")
 127 ipar.setParameter("SMESH_3", "VTKViewer_0_ClippingPlane", "Off")
 128 # set up entry SMESH_5 (mesh2) parameters
 129 ipar.setParameter("SMESH_5", "VTKViewer_0_Visibility", "On")
 130 ipar.setParameter("SMESH_5", "VTKViewer_0_Representation", "1")
 131 ipar.setParameter("SMESH_5", "VTKViewer_0_IsShrunk", "0")
 132 ipar.setParameter("SMESH_5", "VTKViewer_0_Entities", "e:0:f:1:v:0:0d:0:b:0")
 133 ipar.setParameter("SMESH_5", "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")
 134 ipar.setParameter("SMESH_5", "VTKViewer_0_Sizes", "line:1:outline:1:elem0d:5:ball:10:1:shrink:0.75:orientation:0.1:0")
 135 ipar.setParameter("SMESH_5", "VTKViewer_0_PointMarker", "std:1:9")
 136 ipar.setParameter("SMESH_5", "VTKViewer_0_Opacity", "1")
 137 ipar.setParameter("SMESH_5", "VTKViewer_0_ClippingPlane", "Off")
 138
 139
 140 if salome.sg.hasDesktop():
 141   salome.sg.updateObjBrowser(True)
 142   iparameters.getSession().restoreVisualState(1)