Bonjour,
Je travaille sous Scilab (équivalent de Matlab). Je post-traite des maillages issu d'un logiciel type mécanique des fluides.
Le maillage est rempli progressivement au cours du temps (écoulement). Ainsi, je souhaite connaître les indices des éléments (tétraèdres) pour lesquels il y a présence de matière.
Les éléments sont décrits de la façon suivante :
...
ind_element ind_node_1 ind_node_2 ind_node_3 ind_node_4
...
et les résultats (dans la région ou il y a de la matière), pression, température,... :
...
ind_node data_node
...
Je souhaite donc connaître tous les ind_element tels que les ind_node_i sont présents dans les résultats.
Pour chaque colonnes il y a pas unicité de ind_node_i et c'est justement à ce niveau là que les choses m'embètent car dans Scilab il y a bien une fonction intersect mais qui ne donne que l'intersection des éléments uniques...
Je cherche à tout prix à éviter des boucles car la durée du traitement serait excessive.
Merci.
PS : on m'a déjà proposé une solution [dummy, index] = find(kron(a, ones(length(b),1)) - kron(b', ones(a)) == 0) mais elle semble nécessiter des ressources en mémoire trop importantes.
Partager