Salut, On avait trouvé quelques erreurs dans la règle TACPoids, voici un code qui fonctionne. il y a seulement le aimedMetier (lignes 146:151 et lignes 182:183) qui chez moi fait boguer la simulation Si je ne m'abuse le fait de le supprimer ne change pas les resultats, je l ai donc mis en commentaire si qqn peut regarder si ces lignes marchent avec sa pecherie ... merci Sigrid rappel : les metiers interdits sont ceux qui ont l espece visee en capture principale l effort de ces metiers est reporte sur les metiers autorises de la strategie au prorata de l effort deja existant sur ces metiers ou si il n y a pas d autre metier possible l effort est mis dans non activité les autres metiers continuent a la capturer mais rejettent avec une survie = prop survie -- Sigrid LEHUTA ~ ><> ~ Ecologie et Modèles pour l'Halieutique IFREMer, rue de l'ile d'Yeu BP 21105 44311 Nantes Cedex 03 package rules; import java.util.HashSet; import java.util.List; import java.util.Set; import static org.codelutin.i18n.I18n._; import org.apache.commons.logging.Log; import org.apache.commons.logging.LogFactory; import org.codelutin.math.matrix.MatrixFactory; import org.codelutin.math.matrix.MatrixIterator; import org.codelutin.math.matrix.MatrixND; import scripts.ResultName; import scripts.RuleUtil; import scripts.SiMatrix; import fr.ifremer.isisfish.entities.EffortDescription; import fr.ifremer.isisfish.entities.Metier; import fr.ifremer.isisfish.entities.Population; import fr.ifremer.isisfish.entities.Species; import fr.ifremer.isisfish.entities.Strategy; import fr.ifremer.isisfish.entities.StrategyMonthInfo; import fr.ifremer.isisfish.entities.TargetSpecies; import fr.ifremer.isisfish.rule.AbstractRule; import fr.ifremer.isisfish.simulator.MetierMonitor; import fr.ifremer.isisfish.simulator.PopulationMonitor; import fr.ifremer.isisfish.simulator.SimulationContext; import fr.ifremer.isisfish.types.Date; import fr.ifremer.isisfish.types.Month; /** * TAC.java * * Created: 7 septembre 2006 * * @author anonymous <anonymous@labs.libre-entreprise.org> * @version $Revision: 1.3 $ * * Last update: $Date: 290607 $ * by : $Author: sigrid $ */ /** * TAC peut-etre utilisé pour les différents TAC, en proportion des effectifs * et/ou avec survie ou non. * * <li>Pour utiliser le tac proportionnel, il faut mettre dans le parametre propTac * une valeur > 0, le TAC sera alors recalculé a chaque mois de janvier. * <li>Pour utiliser la survie il faut mettre dans le paramètre propSurvie * une valeur > 0, automatiquement les suvie seront ajouté aux effectifs * */ public class TACpoids extends AbstractRule { /** to use log facility, just put in your code: log.info("..."); */ static private Log log = LogFactory.getLog(TACpoids.class); public Species param_species = null; public Date param_beginDate = new Date(0); public Date param_endDate = new Date(90); public double param_propSurvie = 0; public double param_propTac = 0; /** TAC in tonnes */ public double param_tacInTons = 1000000; boolean affectation = false; protected String [] necessaryResult = { // put here all necessary result for this rule // example: // ResultName.MATRIX_BIOMASS, // ResultName.MATRIX_NET_VALUE_OF_LANDINGS_PER_STRATEGY_MET, }; /** * @return the necessaryResult */ public String[] getNecessaryResult() { return this.necessaryResult; } /** * Permet d'afficher a l'utilisateur une aide sur la regle. * @return L'aide ou la description de la regle */ public String getDescription() { return _("TAC weight in tons.\nIf you want survival discard use propSurvie other than 0.\nIf you wish TAC computed as a proportion of the biomass use propTac other than 0."); } /** * Appelé au démarrage de la simulation, cette méthode permet d'initialiser * des valeurs * @param simulation La simulation pour lequel on utilise cette regle */ public void init(SimulationContext context) throws Exception { } /** * La condition qui doit etre vrai pour faire les actions * @param simulation La simulation pour lequel on utilise cette regle * @return vrai si on souhaite que les actions soit faites */ public boolean condition(SimulationContext context, Date date, Metier metier) throws Exception { // read species in current session param_species = (Species)context.getDB().findByTopiaId(param_species.getTopiaId()); // on fait le calcul du tac si nécessaire if (param_propTac > 0 && date.getMonth().equals(Month.JANUARY)){ PopulationMonitor popMon = context.getPopulationMonitor(); param_tacInTons = popMon.getBiomass(param_species) * param_propTac; } boolean result = false; if(date.before(param_beginDate)) { result = false; } else if(date.after(param_endDate)) { result = false; } else { TargetSpecies ts = metier.getMetierSeasonInfo(date.getMonth()).getSpeciesTargetSpecies(param_species); if (ts != null) { double catchTons = RuleUtil.getTotalCatchTons(context, param_species); log.info("[TAC] catchTons = " + catchTons + " >= param_tacInTons:" + param_tacInTons ); if (catchTons >= param_tacInTons) { result = true; } } } return result; } /** * Si la condition est vrai alors cette action est executée avant le pas * de temps de la simulation. * @param simulation La simulation pour lequel on utilise cette regle */ public void preAction(SimulationContext context, Date date, Metier metier) throws Exception { affectation = false; log.info("[TAC] preAction for: " + metier); TargetSpecies ts = metier.getMetierSeasonInfo(date.getMonth()).getSpeciesTargetSpecies(param_species); if(ts != null && ts.getPrimaryCatch()) { // recupere tous les metiers qui ont l'espece en capture principale =>metiers vises // aimedMetiers ne fonctionne pas je ne sais pas pourquoi ! mais au final forbiddenMetier aura le meme effet /*List<Metier> aimedMetiers = context.getDB().find("select distinct(metierSeasonInfo.metier) " + "from fr.ifremer.isisfish.entities.TargetSpecies" + "where species=? and primaryCatch=true", param_species); log.info("aimed Metier: " + aimedMetiers);*/ context.getMetierMonitor().addforbiddenMetier(metier); //récupere toutes les stratégies pratiquant le métier et pour lesquelles la proportion !=0 SiMatrix siMatrix = SiMatrix.getSiMatrix(context); Set<Strategy> strs = new HashSet<Strategy>(); for (Strategy str : siMatrix.getStrategies(date)) { double prop = str.getStrategyMonthInfo(date.getMonth()).getProportionMetier(metier); if (prop != 0) { strs.add(str); } } for (Strategy str : strs){ StrategyMonthInfo smi = str.getStrategyMonthInfo(date.getMonth()); // 1er cas de figure: l'effort est reporte sur un metier de la // meme strategie, n'ayant pas l'espece comme capture principale // et pechant avec le meme engin Set<Metier> possibleMetierCase1 = new HashSet<Metier>(); // second cas de figure: on cherche un metier de substitution // sans condition sur les engins, mais qui soit pratique Set<Metier> possibleMetierCase2 = new HashSet<Metier>(); // 3 eme cas de figure: on cherche des metiers non vises, // sans consideration sur les engins, et pour lesquels la // proportion peut etre nulle Set<Metier> possibleMetierCase3 = new HashSet<Metier>(); for (EffortDescription effort : str.getSetOfVessels().getPossibleMetiers()) { Metier newMetier = effort.getPossibleMetiers(); if ( /*!aimedMetiers.contains(newMetier) &&*/ !metier.getName().equalsIgnoreCase("nonActiviy") && !metier.getName().equalsIgnoreCase("nonActivie") && !metier.getName().equalsIgnoreCase("non Activité") && !context.getMetierMonitor().getForbiddenMetier().contains(newMetier) ) { possibleMetierCase3.add(newMetier); if (smi.getProportionMetier(newMetier) != 0) { possibleMetierCase2.add(newMetier); if (metier.getGear().equals(newMetier.getGear())) { possibleMetierCase1.add(newMetier); } } } } Set<Metier> possibleMetier = null; if (possibleMetierCase1.size() != 0){ log.info("[TAC] Use case 1"); possibleMetier = possibleMetierCase1; } else if (possibleMetierCase2.size() != 0) { log.info("[TAC] Use case 2"); possibleMetier = possibleMetierCase2; } else if (possibleMetierCase3.size() != 0){ log.info("[TAC] Use case 3"); possibleMetier = possibleMetierCase3; } if (possibleMetier != null) { // on repartit maintenant l'effort entre les differents metiers // possibles dans la meme strategie si un metier possible existe // bien la repartion est proportionnelle a l'effort deja alloue // dans la strategie double somme=0; for (Metier met : possibleMetier) { somme += smi.getProportionMetier(met); } for (Metier met : possibleMetier) { double newProportion = smi.getProportionMetier(met) + (smi.getProportionMetier(metier) * smi.getProportionMetier(met) / somme); smi.setProportionMetier(met, newProportion); } smi.setProportionMetier(metier, 0); //le metier vise a alors une proportion nulle log.info("[TAC] il y a des metiers possibles"); } else{ log.info("[TAC] Use no activity"); // sinon on met tout dans le metier nonActivite MetierMonitor metierMon = context.getMetierMonitor(); MatrixND mat = metierMon.getOrCreateNoActivity(date, ResultName.MATRIX_NO_ACTIVITY, siMatrix.getStrategies(date), siMatrix.getMetiers(date)); mat.setValue(str, metier, smi.getProportionMetier(metier)); smi.getProportionMetier().setValue(metier, 0); } } } } /** * Si la condition est vrai alors cette action est executée apres le pas * de temps de la simulation. * @param simulation La simulation pour lequel on utilise cette regle */ public void postAction(SimulationContext context, Date date, Metier metier) throws Exception { log.info("[TAC] postAction for: " + metier); TargetSpecies ts = metier.getMetierSeasonInfo(date.getMonth()).getSpeciesTargetSpecies(param_species); if(ts != null){ if (!affectation){ // ATTENTION // les captures pour cette metapop ne sont plus du qu'au metier // pour qui l'espece est secondaire: elles sont affectees aux // rejets //pb : ne se fait pas par metier // il faut une matrice pour chaques pas de temps qui stocke les // rejets par metier, par metapop et par classes d'age (comme // pour les captures) ////// PopulationMonitor popMon = context.getPopulationMonitor(); log.info("popMon biomass" + popMon.getBiomass(param_species)); for (Population pop : param_species.getPopulation()) { if (!pop.getName().equals("Population_new")){ log.info("pop : " + pop.getName()); // si on a deja une matrice rejet on le vide (elle vient // forcement de la regle taille minimale or si le tac est // atteint, tout va dorenavent dans les rejets et on mais // TOUTES les captures dans les rejets MatrixND discard = popMon.getDiscard(date, pop); log.info("discard : " + discard); if (discard != null) { discard.mults(0); } log.info("catch = "+ popMon.getCatch(pop)); discard = popMon.getCatch(pop).copy(); discard.setName(ResultName.MATRIX_DISCARDS_PER_STR_MET); popMon.addDiscard(date, pop, discard); log.info("[TAC] add discard for " + pop + ": " + discard); if (param_propSurvie > 0) { MatrixND eff = popMon.getN(pop); //on réajoute les survivants aux effectifs for (MatrixIterator i=discard.iterator(); i.next();){ Object [] coord = i.getSemanticsCoordinates(); eff.setValue(coord[2],coord[3], eff.getValue(coord[2],coord[3])+i.getValue()*param_propSurvie); } } } // on a affecte une fois cette meta pop au rejet il ne faut pas // le refaire affectation=true; } } } } }