29 | *
30 | */
31 | public class PUMASimilarityTable {
32 |
33 | protected static Logger logger = Logger.getLogger("com.pb.models.synpopV3");
34 | protected TableDataSet table;
35 |
36 | public PUMASimilarityTable(TableDataReader designTableReader){
37 | table=designTableReader.getTable("PUMASimilarityTable.csv");
38 | }
39 |
40 | /**
41 | * Given a PUMA get its similar PUMA of a given rank.
42 | * @param puma represents the given PUMA.
43 | * @param rank represents the desired rank of the similar PUMA. If rank=1, similar PUMA is PUMA itself.
44 | * @return
45 | */
46 | public int getSimilarPUMA(int puma, int rank) {
47 | //Initialize similar PUMA as PUMA itself.
48 | //If a similar PUMA of a desired rank is not found, then return PUMA itself as its similar PUMA.
49 | int result = puma;
50 |
51 | int NRows=table.getRowCount();
52 | int[] currentPUMA=table.getColumnAsInt("PUMA_1");
53 | int[] similarPUMA=table.getColumnAsInt("PUMA_2");
54 | int[] ranks=table.getColumnAsInt("SIMILARRAN");
55 |
56 | //find the similar PUMA
57 | for (int i = 0; i < NRows; i++) {
58 | if ((puma == currentPUMA[i]) && (rank == ranks[i])) {
59 | result = similarPUMA[i];
60 | break;
61 | }
62 | }
63 |
64 | return result;
65 | }
66 |
67 | public static void main(String [] args){
68 | TableDataReader designReader=new TableDataReader("design");
69 | PUMASimilarityTable table=new PUMASimilarityTable(designReader);
70 | logger.info("ok, I am done.");
71 | }
72 |
73 | }
74 |
--------------------------------------------------------------------------------
/population_synthesis/2__synpop/model/synpop_source_code/src/java/com/pb/models/synpopV3/PUMSAttrs.java:
--------------------------------------------------------------------------------
1 | /*
2 | * Copyright 2006 PB Consult Inc.
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package com.pb.models.synpopV3;
17 |
18 | import java.util.Vector;
19 | import org.apache.log4j.Logger;
20 |
21 | /**
22 | * Company: PB Consult, Parsons Brinckerhoff
23 | * @author Wu Sun
24 | * @version 1.1, modified on Nov. 15, 2004
25 | *
26 | * Represents attributes of PUMSHH, DerivedHH, PUMSPerson, and DerivedPerson.
27 | * Attributes are read in from property file.
28 | *
29 | */
30 |
31 | public class PUMSAttrs {
32 |
33 | protected static Logger logger;
34 | protected static Vector HHAttrs;
35 | protected static Vector DerivedHHAttrs;
36 | protected static Vector PersonAttrs;
37 | protected static Vector DerivedPersonAttrs;
38 |
39 | static {
40 | logger= Logger.getLogger("com.pb.models.synpopV3");
41 | HHAttrs=PropertyParser.getPropertyElementsByName("pums.hhattrs",",");
42 | DerivedHHAttrs=PropertyParser.getPropertyElementsByName("pums.derivedHHAttrs",",");
43 | PersonAttrs=PropertyParser.getPropertyElementsByName("pums.pattrs",",");
44 | DerivedPersonAttrs=PropertyParser.getPropertyElementsByName("pums.derivedPersonAttrs",",");
45 | }
46 |
47 | public static Vector getHHAttrs(){
48 | return HHAttrs;
49 | }
50 |
51 | public static Vector getPersonAttrs(){
52 | return PersonAttrs;
53 | }
54 |
55 | public static Vector getDerivedHHAttrs(){
56 | return DerivedHHAttrs;
57 | }
58 |
59 | public static Vector getDerivedPersonAttrs(){
60 | return DerivedPersonAttrs;
61 | }
62 |
63 | //for testing purpose only
64 | //successfully tested on March 23, 2005
65 | public static void main(String [] args){
66 | Vector v1=PUMSAttrs.getHHAttrs();
67 | Vector v2=PUMSAttrs.getPersonAttrs();
68 | Vector v3=PUMSAttrs.getDerivedHHAttrs();
69 | System.out.println("ok, I am done");
70 | }
71 | }
--------------------------------------------------------------------------------
/population_synthesis/2__synpop/model/synpop_source_code/src/java/com/pb/models/synpopV3/PUMSBucketBin.java:
--------------------------------------------------------------------------------
1 | /*
2 | * Copyright 2006 PB Consult Inc.
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package com.pb.models.synpopV3;
17 |
18 | import java.util.Vector;
19 |
20 | /**
21 | * Company: PB Consult, Parsons Brinckerhoff
22 | * @author Wu Sun
23 | * @version 1.0, Jan. 8, 2004
24 | */
25 |
26 | public class PUMSBucketBin {
27 |
28 | //HHs in this bin (a vector of DerivedHH objects)
29 | protected Vector hhs;
30 | //number of HHs in this bin
31 | protected int NBinMembers;
32 | //times all HHs in this been has been drawn
33 | protected int NEmptied;
34 |
35 |
36 | public PUMSBucketBin() {
37 | hhs=new Vector();
38 | }
39 |
40 | public void update(){
41 | if(getNDrawn()==NBinMembers){
42 | NEmptied++;
43 | //set drawn status of all HHs in this bin to "false"
44 | setHHsDrawnSta(false);
45 | }
46 | }
47 |
48 | public void reset(){
49 | DerivedHH hh=null;
50 | NEmptied=0;
51 | for(int i=0; iCompany: PB Consult, Parsons Brinckerhoff
24 | * @author Wu Sun
25 | * @version 1.1, modified on Nov. 15, 2004
26 | *
27 | * Create a PUMSPerson object from raw string record of a person in PUMS data.
28 | */
29 |
30 | public class PUMSPerson {
31 |
32 | //raw person record as a string
33 | protected String pRecord_raw;
34 | //PUMS person attributes
35 | protected Vector pattrs;
36 | //person attributes map, attribute name<--->value pairs
37 | protected HashMap pAttrsMap;
38 | //PUMS person record parser
39 | protected PUMSRecordParser recordParser;
40 | //data dictionary
41 | protected DataDictionary dd;
42 |
43 | /**
44 | * constructor
45 | * must be defined so that subclass DerivedPerson compiles
46 | */
47 | public PUMSPerson(){
48 | }
49 | /**
50 | * Constructor
51 | * @param pRecord_raw represents a raw string person record
52 | * @param dd represents a PUMS data dictionary
53 | */
54 | public PUMSPerson(String pRecord_raw, DataDictionary dd) {
55 | this.pRecord_raw=pRecord_raw;
56 | this.dd=dd;
57 | recordParser=new PUMSRecordParser(pRecord_raw, dd);
58 | pattrs=PUMSAttrs.getPersonAttrs();
59 | makeAttrsMap();
60 | }
61 |
62 | /**
63 | * Givan a person attribute name, return its value
64 | * @param attr represents a person attribute name
65 | * @return
66 | */
67 | public int getPAttr(String attr){
68 | int attrVal=((Integer)pAttrsMap.get(attr)).intValue();
69 | return attrVal;
70 | }
71 |
72 | /**
73 | * get all person attributes
74 | * @return
75 | */
76 | public HashMap getPAttrs(){
77 | return pAttrsMap;
78 | }
79 |
80 | private void makeAttrsMap(){
81 |
82 | String attr;
83 | int attrVal;
84 | pAttrsMap=new HashMap();
85 |
86 | for(int i=0; iCompany: PB Consult, Parsons Brinckerhoff
23 | * @author Wu Sun
24 | * @version 1.2, modified on Nov. 17, 2004
25 | *
26 | * Use this class to parse PUMS HH record.
27 | */
28 |
29 | public class PUMSRecordParser {
30 |
31 | protected String pumsDir;
32 | protected String pumsDD;
33 | protected DataDictionary dd;
34 | protected String record;
35 |
36 | /**
37 | * constructor
38 | * @param record represents a raw record string
39 | * @param dd represents a data dictionary
40 | */
41 | public PUMSRecordParser(String record, DataDictionary dd){
42 | this.record=record;
43 | this.dd=dd;
44 | pumsDir=PropertyParser.getPropertyByName("pums.directory");
45 | pumsDD=PropertyParser.getPropertyByName("pums.dictionary");
46 | }
47 |
48 | /**
49 | * Given a string PUMSHH record and a PUMSHH variable name, return value of this variable
50 | * @param PUMSVariable represents PUMSHH variable name
51 | * @return
52 | */
53 | public int getPUMSHHDataValue(String PUMSVariable) {
54 | int startCol=dd.getStartCol(dd.HHAttribs, PUMSVariable);
55 | int endCol=dd.getLastCol(dd.HHAttribs, PUMSVariable);
56 | String temp = record.substring(startCol,endCol);
57 |
58 | String emptyStr=" ";
59 | for(int i=0; iCompany: PB Consult, Parsons Brinckerhoff
24 | * @author Wu Sun
25 | * @version 1.0, Jan. 9, 2004
26 | */
27 |
28 | public class RandomSeedMatrix {
29 | protected static Logger logger = Logger.getLogger("com.pb.models.synpopV3");
30 | protected static int NTAZ;
31 | protected static int NHHCat;
32 | //randome seed matrix, 1st dimension number of HHCat, 2nd dimension number of TAZ
33 | protected static long [][] matrix;
34 | protected static long seed;
35 |
36 | static {
37 | NTAZ=PopulationSynthesizer.numberOfInternalZones;
38 | NHHCat=HHCatTable.NHHCat;
39 | matrix=new long [NHHCat][NTAZ];
40 | seed=(new Long(PropertyParser.getPropertyByName("RandomSeed"))).longValue();
41 | logger.info("Random seed set to: " + seed);
42 | setSeedMatrix();
43 | }
44 |
45 | /**
46 | * Get random seed matrix. First dimension TAZ, second dimension HHCat.
47 | * @return
48 | */
49 | public static long [][] getSeedMatrix(){
50 | return matrix;
51 | }
52 |
53 | /**
54 | * Get a random seed from one cell in the random matrix.
55 | * @param TAZ represents the internal TAZ ID (starts from 0).
56 | * @param HHCat represents the internal HHCat ID (starts from 0).
57 | * @return
58 | */
59 | public static long getSeed(int HHCat, int TAZ){
60 | return matrix[HHCat][TAZ];
61 | }
62 |
63 | /**
64 | * Populate random seed matrix.
65 | */
66 | private static void setSeedMatrix(){
67 | Random generator=new Random(seed);
68 | for(int i=0; i
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
--------------------------------------------------------------------------------
/population_synthesis/3__hh_taz_to_maz/cmap_hh_maz.R:
--------------------------------------------------------------------------------
1 | #Allocate SynPop HHs to MAZs
2 | #Ben Stabler, stabler@pbworld.com, 05/24/13
3 | #Revised to work with TG files by NMP, 5/19/16
4 |
5 | source("../get_dir.R") # Intelligently create PopSynDir variable
6 | inDir1 = file.path(PopSynDir, "2__synpop/outputs")
7 | inDir2 = file.path(PopSynDir, "3__hh_taz_to_maz/inputs")
8 | outDir = file.path(PopSynDir, "3__hh_taz_to_maz/outputs")
9 |
10 | #Read input files
11 | setwd(inDir1)
12 | hh = read.csv("ForecastHHFile.csv")
13 | setwd(inDir2)
14 | #mazs = read.csv("SubzoneData.csv")
15 | hh_in = read.table("HH_IN.TXT", header=FALSE, sep=",",
16 | col.names=c("subzone09","hshld","ahh","whh","chh",
17 | "iqi","ahi","automs","pef"))
18 | hh_in = subset(hh_in, select=c("subzone09","hshld"))
19 |
20 | geog_in = read.table(file="GEOG_IN.TXT", header=FALSE, sep=",",
21 | col.names=c("subzone09","cofips","coname","state","puma1",
22 | "puma5","zone09","chi","cbd","rowcol","sqmi"))
23 | geog_in = subset(geog_in, select=c("subzone09","zone09"))
24 |
25 | mazs = merge(geog_in, hh_in)
26 |
27 | #Get mazs (and their num hhs) by taz
28 | mazsByTaz = tapply(mazs$subzone09, mazs$zone09, function(x) x)
29 | mazHHsByTaz = tapply(mazs$hshld, mazs$zone09, function(x) x)
30 |
31 | #Randomly assign an MAZ to a TAZ using the MAZ num HHs as the probability weight
32 | assignMaz = function(taz) {
33 | mazs = mazsByTaz[taz][[1]]
34 | hhs = mazHHsByTaz[taz][[1]]
35 | ifelse(length(mazs)>1, sample(mazs, 1, T, hhs), mazs)
36 | }
37 |
38 | #Get hh tazs and assign mazs
39 | tazs = as.character(hh$TAZ)
40 | mazs = sapply(tazs, assignMaz)
41 |
42 | #Write result
43 | hh$maz = mazs
44 | setwd(outDir)
45 | write.csv(hh, "ForecastHHFile_maz.csv", quote=F, row.names=F)
46 |
--------------------------------------------------------------------------------
/population_synthesis/3__hh_taz_to_maz/distribute_hh_to_maz.bat:
--------------------------------------------------------------------------------
1 | SET Rscript="C:\Program Files\R\R-3.2.0\bin\x64\Rscript.exe"
2 | CALL %Rscript% cmap_hh_maz.R
3 | PAUSE
--------------------------------------------------------------------------------
/population_synthesis/README.TXT:
--------------------------------------------------------------------------------
1 | HOW TO RUN ABM POPULATION SYNTHESIS
2 | -----------------------------------
3 |
4 | 0. OPTIONAL: If control totals are being run updated with new PUMS data for
5 | the first time, uncomment the large section of code in 1__control_files/
6 | createNewControls.R that creates countsPUMS.csv and distPUMS.csv in
7 | 1__control_files/inputs. These new CSVs can then be reused for any future
8 | population synthesis based on the same PUMS data. (NOTE: the 2007-2011 ACS
9 | was the last one where all of the PUMAs match the 2000 PUMAs that the
10 | current model was estimated with, so it is not recommended to use newer
11 | PUMS data without first updating the PopSyn model itself.)
12 |
13 | 1. a) FOR BASE YEAR (2010) POPSYN:
14 | Update any inputs in 1__control_files (replacing oldControls.csv with the
15 | previous version of forecastControls.csv from 2__synpop/inputs/datatables),
16 | then submit update_controls.bat. The output file will be saved to
17 | 2__synpop/inputs/datatables/forecastControls.csv.
18 |
19 | b) FOR FUTURE YEAR POPSYN:
20 | Apply base year control variables' rates to future year's household counts,
21 | and save the results to 2__synpop/inputs/datatables/forecastControls.csv.
22 |
23 | 2. Submit 2__synpop/model/RunCMAPPopSynForecast.bat. This requires the Java
24 | Runtime Environment to be installed, so it is recommended that this be done
25 | from one of the modeling servers that have it already (e.g. cmap-mod04).
26 |
27 | 3. Replace the GEOG_IN.TXT and HH_IN.TXT files in 3__hh_taz_to_maz/inputs if
28 | they are out of date. These are taken directly from the 4-step Trip
29 | Generation model.
30 |
31 | 4. Submit 3__hh_taz_to_maz/distribute_hh_to_maz.bat. This will assign the
32 | zone-level households produced by ARCPopSyn to subzones.
33 |
34 | 5. Copy the following CSVs to the main ABM inputs folder for new model runs:
35 | - 2__synpop/outputs/ForecastPersonFile.csv
36 | - 3__hh_taz_to_maz/outputs/ForecastHHFile_MAZ.csv
37 |
--------------------------------------------------------------------------------
/population_synthesis/get_dir.R:
--------------------------------------------------------------------------------
1 | # When this script is called with "source()", it will save its parent directory
2 | # (i.e. the root population synthesis directory) to a variable called PopSynDir
3 | # that can subsequently be used to define other paths in a relative manner.
4 | PopSynDir <- normalizePath(dirname(parent.frame(2)$ofile))
--------------------------------------------------------------------------------
/results_analysis/abm/__init__.py:
--------------------------------------------------------------------------------
1 | from abm import *
2 | from comparison import *
3 | import to_csv
4 | import print_summary
5 |
--------------------------------------------------------------------------------
/results_analysis/abm/comparison.py:
--------------------------------------------------------------------------------
1 | #!/usr/bin/env python
2 | '''
3 | comparison.py
4 | Author: npeterson
5 | Revised: 3/16/16
6 | ---------------------------------------------------------------------------
7 | A class for comparing the contents of two ABM object databases.
8 |
9 | '''
10 | import copy
11 | import abm
12 |
13 |
14 | class Comparison(object):
15 | ''' A class for comparing two ABM objects. Initialized with a base and a
16 | test ABM object. All comparisons are relative to the base. '''
17 |
18 | # --- Init ---
19 | def __init__(self, base_abm, test_abm):
20 | self.base = base_abm
21 | self.test = test_abm
22 | return None
23 |
24 |
25 | def __str__(self):
26 | return '[Comparison: BASE {0}; TEST {1}]'.format(self.base, self.test)
27 |
28 |
29 | # --- Instance methods ---
30 | def _open_dbs(self):
31 | ''' Open base & test ABM database connections. '''
32 | self.base._open_db()
33 | self.test._open_db()
34 | return None
35 |
36 |
37 | def _close_dbs(self):
38 | ''' Close base & test ABM database connections. '''
39 | self.base._close_db()
40 | self.test._close_db()
41 | return None
42 |
--------------------------------------------------------------------------------
/results_analysis/results_analysis.py:
--------------------------------------------------------------------------------
1 | #!/usr/bin/env python
2 | '''
3 | results_analysis.py
4 | Author: npeterson
5 | Revised: 3/16/16
6 | ---------------------------------------------------------------------------
7 | A script for printing a series of summary statistics for base and test ABM
8 | objects and a Comparison object created from them.
9 |
10 | '''
11 | from abm import *
12 |
13 | # Set Base ABM parameters
14 | BASE_DIR = r'X:\new_server_tests\test_5pct'
15 | BASE_PCT = 0.05
16 |
17 | # Set Test ABM parameters
18 | TEST_DIR = r'X:\new_server_tests\test_50pct_newsocec_3sp'
19 | TEST_PCT = 0.5
20 |
21 |
22 | # Initialize Base ABM object & print summaries
23 | print '\n{0:*^50}'.format(' P R O C E S S I N G ')
24 | print '\n{0:=^50}\n'.format(' BASE NETWORK ')
25 | base = ABM(BASE_DIR, BASE_PCT)
26 | print base
27 | print_summary.mode_share(base)
28 | print_summary.transit_stats(base)
29 | print_summary.ptrips_by_class(base)
30 | print_summary.vmt_by_speed(base)
31 | print ' '
32 |
33 |
34 | # Initialize Test ABM object & print summaries
35 | print '\n{0:=^50}\n'.format(' TEST NETWORK ')
36 | test = ABM(TEST_DIR, TEST_PCT)
37 | print test
38 | print_summary.mode_share(test)
39 | print_summary.transit_stats(test)
40 | print_summary.ptrips_by_class(test)
41 | print_summary.vmt_by_speed(test)
42 | print ' '
43 |
44 |
45 | # Initialize Test vs. Base Comparison object
46 | print '\n{0:=^50}\n'.format(' COMPARISON ')
47 | comp = Comparison(base, test)
48 | print comp
49 | print ' '
50 |
51 |
52 | # Print Comparison summaries
53 | print '\n{0:*^50}'.format(' R E S U L T S ')
54 | print_summary.mode_share_change(comp)
55 | print_summary.transit_stats_change(comp)
56 | print_summary.ptrips_by_class_change(comp)
57 | print_summary.vmt_by_speed_change(comp)
58 | print_summary.new_all(comp)
59 | print_summary.new_auto(comp)
60 | print_summary.new_dtt(comp)
61 | print_summary.new_wtt(comp)
62 | print_summary.new_other(comp)
63 |
--------------------------------------------------------------------------------
/runCTRAMP-SingleProcess.bat:
--------------------------------------------------------------------------------
1 |
2 | rem ############ PARAMETERS ############
3 | set RUNTIME=Y:/{{{TEMPLATE}}}/cmap_abm
4 | set JAVA_PATH=C:/Program Files/Java/jdk1.7.0_13
5 |
6 | rem ########## Run CTRAMP ##############
7 | "%JAVA_PATH%\bin\java" -Xmx120000m -cp "%RUNTIME%;%RUNTIME%/exec/jxl.jar;%RUNTIME%/exec/jppf-5.0.4-lib/*;%RUNTIME%/exec/*;config" -Dlog4j.configuration=log4j.xml -server -Djppf.config=jppf-guiLocal.properties com.pb.cmap.tourBased.CmapTourBasedModel "cmap" -iteration 1 -sampleRate %1 -sampleSeed 0
8 |
9 | rem ########## Run CTRAMP Remote Debug ##
10 | rem "%JAVA_PATH%\bin\java" -Xdebug -Xrunjdwp:transport=dt_socket,address=1044,server=y,suspend=y -Xmx120000m -cp "%RUNTIME%;%RUNTIME%/exec/jxl.jar;%RUNTIME%/exec/jppf-5.0.4-lib/*;%RUNTIME%/exec/*;config" -Dlog4j.configuration=log4j.xml -server -Djppf.config=jppf-guiLocal.properties com.pb.cmap.tourBased.CmapTourBasedModel "cmap" -iteration 1 -sampleRate %1 -sampleSeed 0
11 |
--------------------------------------------------------------------------------
/runCTRAMP.bat:
--------------------------------------------------------------------------------
1 |
2 | rem ############ PARAMETERS ############
3 | set RUNTIME=Y:/{{{TEMPLATE}}}/cmap_abm
4 | set JAVA_PATH=C:/Program Files/Java/jdk1.7.0_13
5 |
6 | rem ########## Run CTRAMP ##############
7 | "%JAVA_PATH%\bin\java" -Xmx8000m -cp "%RUNTIME%;%RUNTIME%/exec/jxl.jar;%RUNTIME%/exec/jppf-5.0.4-lib/*;%RUNTIME%/exec/*;config" -Dlog4j.configuration=log4j.xml -server -Djppf.config=jppf-guiDistributed.properties com.pb.cmap.tourBased.CmapTourBasedModel "cmap" -iteration 1 -sampleRate %1 -sampleSeed 0
8 |
--------------------------------------------------------------------------------
/runMAZSkimsInitial.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Run MAZ Skimming Initial Setup
3 | # Ben Stabler, stabler@pbworld.com, 10/29/12
4 | # Roshan Kumar, kumarr1@pbworld.com
5 | # "C:\Python26\P64\python.exe" runMAZSkimsInitial
6 | #############################################################################
7 |
8 | #inputs
9 | isBaseYear = True
10 | navteqShpFile = "inputs/Chicago_streets.shp"
11 | #SPThreshold = 15*5280 #max dist in feet
12 | SPThreshold = 5*5280 #max dist in feet
13 | maxMAZID = 70000
14 | minTAPID = 60000
15 | maxCenID = 16819
16 | rsThreshold = 9
17 | landUseTypes = 12
18 | centroidsMAZFile = "inputs/CMAP_MAZ_cents.txt"
19 | rsToSZInput = "inputs/cmapRStoSubZone.txt"
20 | tapFile = "inputs/tap_attributes.csv"
21 | spMAZ2TAP = "outputs/SP_MAZ_to_MAZ_and_TAP.txt"
22 | nodesinMAZFile = "inputs/nodes_in_MAZ_CMAP.txt"
23 |
24 | #additional parameters for future year maz skims
25 | nodeDistFromHwy = "inputs/cmapDistFromHwy.txt" #for future build MAZs
26 | mazLUFile = "inputs/cmapMAZLandUseShare.txt" #for future build MAZs
27 | spInsideRS8 = "outputs/SP_with_TAP_clean_RS8_LU_Share_Dist.txt" #for future build MAZs
28 | spOutsideRS8 = "outputs/SP_with_TAP_clean_outside_RS8_LU_Share_Dist.txt" #for future build MAZs
29 |
30 | ############################################################################
31 |
32 | #intermediate settings
33 | import os, csv, math, shutil
34 | root = os.getcwd()
35 | int_output = "outputs"
36 | connectorEDist = 2 * SPThreshold
37 | minTAPid = 2 * minTAPID
38 | SPThresholdConn = 3 * SPThreshold
39 | SPThresholdAct = SPThreshold
40 | number_of_cores = int(os.environ["NUMBER_OF_PROCESSORS"])
41 | linkFile_int = "outputs/cmap_NavTechLinks_no_freeway.txt" #without connectors
42 | nodeFile_int = "outputs/cmap_NavTechNodes_no_freeway.txt" #without centroids
43 | outputSPFile = "outputs/centToCentDist.txt"
44 | tapConnectorFile = "outputs/tapConnectors.txt"
45 | connectorFile = "outputs/CMAP_connectors.txt"
46 | linkFile = "outputs/linkFile_Final.txt" #with connectors
47 | nodeFile = "outputs/nodeFile_Final.txt" #with centroids
48 |
49 | ############################################################################
50 |
51 | def main():
52 |
53 | #copy itself to the scripts folder for other scripts to reference
54 | shutil.copyfile("runMAZSkimsInitial.py", "scripts/parameters.py")
55 |
56 | print("run initial network build procedures")
57 | from scripts import cmapTransportationNetwork
58 | from scripts import cmapTapConnector
59 | from scripts import cmapCentroidsConnectors
60 | from scripts import cmapInputFileGen
61 | if isBaseYear:
62 | from scripts import cmapShortestPath_NX
63 | else:
64 | from scripts import cmapShortestPath_NX_future
65 |
66 | print("read shapefile and generate csv file of navteq links")
67 | cmapTransportationNetwork.getLinkDirectionality()
68 | print("read shapefile and generate csv file of navteq nodes")
69 | cmapTransportationNetwork.createAllPoints()
70 | print("create MAZ centroid connectors")
71 | cmapCentroidsConnectors.createNetworkWithCentroidConnectors()
72 | print("create TAP connectors")
73 | cmapTapConnector.createNetworkWithCentroidConnectors()
74 | print("merge node and link files to one node and one link file containing all elements")
75 | cmapInputFileGen.createNetworkInput()
76 | if isBaseYear:
77 | print("pickle network for later use")
78 | cmapShortestPath_NX.pickleNetwork()
79 | else:
80 | print("pickle network for later use")
81 | cmapShortestPath_NX_future.pickleNetwork()
82 | print("initial shortest path network build complete")
83 |
84 | if __name__ == '__main__':
85 | main()
86 |
--------------------------------------------------------------------------------
/runModel.bat:
--------------------------------------------------------------------------------
1 | :: runModel.bat
2 | :: Written by BTS, 06/24/13
3 | :: Modified by NMP, 04/14/16 - run accessibilities after skimming
4 | :: 08/12/16 - supply scenario as parameter to scripts
5 | :: ----------------------------------------------------------------------------
6 | :: Run the CMAP ABM from start to finish.
7 | :: Call from command line, e.g.: runModel.bat > blog.txt
8 |
9 | :: Parameters
10 | SET projDir=Y:/{{{TEMPLATE}}}/cmap_abm
11 | SET sampleRate=0.50
12 | SET scen=100
13 |
14 | :: Set Python paths
15 | CALL EmmeConfig.bat
16 | SET emmepy="%EMMEPATH%\Python27\python.exe"
17 | SET py64="C:\Python27\python.exe"
18 |
19 | :: Optional parameters for distributed skimming/assignment.
20 | :: If enabled, call script with username and password,
21 | :: e.g.: runModel.bat npeterson P@ssw0rd > blog.txt
22 | REM @SET user=%1
23 | REM @SET pwd=%2
24 | REM SET un=cmap\%user%
25 | REM SET IP1=\\10.10.1.52
26 | REM SET IP2=\\10.10.1.63
27 |
28 | :: Initialize timestamp log
29 | if exist model_run_timestamp.txt (del model_run_timestamp.txt /Q)
30 | @ECHO ====================================================== >> model_run_timestamp.txt
31 | @ECHO BEGIN CMAP REGIONAL MODEL RUN - SAMPLE RATE %sampleRate% >> model_run_timestamp.txt
32 | @ECHO Model Run Start Time: %date% %time% >> model_run_timestamp.txt
33 | @ECHO ====================================================== >> model_run_timestamp.txt
34 |
35 | :: Create MAZ skims
36 | @ECHO Create MAZ skims: %date% %time% >> model_run_timestamp.txt
37 | %py64% runMAZSkimsInitial.py
38 | %py64% scripts/SPwrapper.py
39 | %py64% scripts/cmapPostProcess.py
40 |
41 | :: Create networks and skim
42 | @ECHO Build networks: %date% %time% >> model_run_timestamp.txt
43 | %emmepy% runBuildNetworks.py 1,1,1,1,1,1,1,1 %scen%
44 | @ECHO Run initial skims: %date% %time% >> model_run_timestamp.txt
45 | %emmepy% runInitialSkims.py 1,1,1,1,1,1,1,1 %scen%
46 |
47 | :: Run accessibilities script
48 | @ECHO Calculate accessibilities: %date% %time% >> model_run_timestamp.txt
49 | cd accessibilities
50 | CALL runAccessibilities.bat %scen%
51 | copy outputs\accessibility_maz.csv ..\inputs
52 | cd ..
53 |
54 | :: Create TAP lines file and run CT-RAMP
55 | @ECHO Create TAP lines: %date% %time% >> model_run_timestamp.txt
56 | %emmepy% runTapLines.py %scen%
57 | @ECHO Run CT-RAMP: %date% %time% >> model_run_timestamp.txt
58 | CALL runCTRAMP-SingleProcess.bat %sampleRate%
59 |
60 | :: Final skimming and assignments
61 | @ECHO Run final assignments: %date% %time% >> model_run_timestamp.txt
62 | %emmepy% runFinalAssignments.py 1,1,1,1,1,1,1,1 %scen%
63 |
64 | @ECHO ====================================================== >> model_run_timestamp.txt
65 | @ECHO END CMAP REGIONAL MODEL RUN >> model_run_timestamp.txt
66 | @ECHO Model Run End Time: %date% %time% >> model_run_timestamp.txt
67 | @ECHO ====================================================== >> model_run_timestamp.txt
--------------------------------------------------------------------------------
/runTVPB.bat:
--------------------------------------------------------------------------------
1 |
2 | rem #Run CMAP Transit Virtual Path Builder
3 | rem #Ben Stabler, stabler@pbworld.com, 02/12/13
4 |
5 | rem ########## Parameters ############
6 | set RUNTIME=%CD%
7 | set HOST_IP=localhost
8 |
9 | rem ########## Run TVPB ##############
10 | java -Xmx40000m -cp "%RUNTIME%;%RUNTIME%/exec/jxl.jar;%RUNTIME%/exec/*;exec" -Dlog4j.configuration=config/log4j.xml -server com.pb.cmap.tvpb.TransitVirtualPathBuilder "cmap"
11 |
12 | rem ########## Run TVPB Remote Debug ##
13 | rem java -Xdebug -Xrunjdwp:transport=dt_socket,address=1044,server=y,suspend=y -Xmx40000m -cp "%RUNTIME%;%RUNTIME%/exec/jxl.jar;%RUNTIME%/exec/*;exec" -Dlog4j.configuration=config/log4j.xml -server com.pb.cmap.tvpb.TransitVirtualPathBuilder "cmap"
14 |
15 |
--------------------------------------------------------------------------------
/runTapLines.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Create lines by tap file for CT-RAMP
3 | # Ben Stabler, stabler@pbworld.com, 02/18/13
4 | # "C:\Program Files (x86)\INRO\Emme\Emme 4\Emme-4.0.3\Python26\python.exe" runTapLines.py
5 | #############################################################################
6 |
7 | #load libraries
8 | import inro.modeller as m
9 | import inro.emme.desktop.app as d
10 | import inro.emme.prompt as p
11 | import os, csv, math, pickle, datetime, sys
12 |
13 | #EMME project file
14 | empFile = "CMAP-ABM/CMAP-ABM.emp"
15 |
16 | #scenarios
17 | transitImport = int(sys.argv[1]) #100
18 | tods = [1,2,3,4,5,6,7,8]
19 | highwayScenarios = [i for i in tods]
20 | transitScenarios = [transitImport + i for i in tods]
21 |
22 | #settings
23 | tapFile = "inputs/tap_attributes.csv"
24 | tapLinks = "inputs/tap2node.csv"
25 | tapLinesOutFileName = "outputs/tapLines.csv"
26 |
27 | ############################################################################
28 |
29 | #start EMME desktop and attach a modeller session
30 | desktop = d.start_dedicated(True, "cmap", empFile)
31 | m = m.Modeller(desktop)
32 |
33 | ############################################################################
34 |
35 | #create tap virtual zones
36 | taps = []
37 | with open(tapFile, 'rb') as csvfile:
38 | tapreader = csv.reader(csvfile, skipinitialspace=True)
39 | for row in tapreader:
40 | taps.append(row)
41 | taps_col_names = taps.pop(0)
42 |
43 | #create access links
44 | access_links = []
45 | with open(tapLinks, 'rb') as csvfile:
46 | tapreader = csv.reader(csvfile, skipinitialspace=True)
47 | for row in tapreader:
48 | access_links.append(row)
49 | access_links_col_names = access_links.pop(0)
50 |
51 | ############################################################################
52 |
53 | #generates lines by tap file to trim tap set for use in CT-RAMP
54 |
55 | #get all transit nodes from EMME
56 | linesByNode = dict()
57 | for tranScen in transitScenarios:
58 |
59 | #get transit segments for scenario
60 | network = m.emmebank.scenario(tranScen).get_network()
61 | tsegs = network.transit_segments()
62 |
63 | for tseg in tsegs:
64 | line = tseg.line.id
65 | from_node = tseg.i_node.id
66 |
67 | if from_node not in linesByNode:
68 | linesByNode[from_node] = set()
69 | linesByNode[from_node].add(line)
70 |
71 | to_node = tseg.j_node
72 | if to_node != None:
73 | to_node = tseg.j_node.id
74 | if to_node not in linesByNode:
75 | linesByNode[to_node] = set()
76 | linesByNode[to_node].add(line)
77 |
78 | #get nodes connected to each tap
79 | nodesByTap = dict()
80 | for linkRow in access_links:
81 | tap = linkRow[0]
82 | node = linkRow[1]
83 |
84 | if tap not in nodesByTap:
85 | nodesByTap[tap] = set()
86 | nodesByTap[tap].add(node)
87 |
88 | #get lines for each tap
89 | linesByTap = dict()
90 | for tapRow in taps:
91 | tap = tapRow[0]
92 | nodes = nodesByTap[tap]
93 |
94 | if tap not in linesByTap:
95 | linesByTap[tap] = set()
96 | for node in nodes:
97 | if node in linesByNode:
98 | lines = linesByNode[node]
99 | for line in lines:
100 | linesByTap[tap].add(line)
101 |
102 | #write out tapLines file for CT-RAMP
103 | f = file(tapLinesOutFileName,"wt")
104 | f.write("TAP,LINES\n")
105 | for tap in linesByTap.keys():
106 | lines = " ".join(list(linesByTap[tap]))
107 | if lines != "":
108 | f.write("%s,%s\n" % (tap,lines))
109 | f.close()
110 |
111 | #log results
112 | print("TapLines Complete %s" % (datetime.datetime.now()))
113 |
--------------------------------------------------------------------------------
/scripts/HwayMatIn.mac:
--------------------------------------------------------------------------------
1 | ~o=295
2 | reports=nul
3 | ~/Directory where the source emmebank is in
4 | ~t1=%1%
5 | ~/Tranfer matrices from the external emmebank
6 | 1.31
7 | ~+|1|%t1%
8 | ~+|2|3|2|mf01,mf99|y|y
9 | 4
10 | 9
--------------------------------------------------------------------------------
/scripts/SPwrapper.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Shortest Path code for CMAP non-motorized modes
3 | # Roshan Kumar, kumarr1@pbworld.com, 12/20/2012
4 | #############################################################################
5 |
6 | #Should always be placed in the same folder as cmapPostProcess.py
7 | #Wrapper code that distributes the shortest path computations to the different cores.
8 | #The total # of parallel processes equal to the number of processors
9 |
10 | import parameters, time
11 | from multiprocessing import Process
12 |
13 | if parameters.isBaseYear:
14 | import cmapShortestPath_NX
15 | else:
16 | import cmapShortestPath_NX_future
17 |
18 | ################################################################
19 |
20 | #parallel distribution of SP computation to different processors
21 | def runSPinParallel():
22 | #alpha = 1
23 | if __name__ == '__main__':
24 |
25 | for num in range(parameters.number_of_cores):
26 |
27 | if parameters.isBaseYear:
28 | p = Process(target=cmapShortestPath_NX.cenToCenSP, args=(num,))
29 | p.start()
30 | else:
31 | p = Process(target=cmapShortestPath_NX_future.cenToCenSP, args=(num,))
32 | p.start()
33 |
34 | ################################################################
35 |
36 | print "SP code start time: ", time.time()
37 | runSPinParallel()
38 |
--------------------------------------------------------------------------------
/scripts/TranMatIn.mac:
--------------------------------------------------------------------------------
1 | ~o=295
2 | reports=nul
3 | ~/Input transit trip tables for time period
4 | / Directory
5 | ~t1=../../inputs/
6 | ~x=%1%
7 | ~+|batchin=%t1%Conv%x%_%x%467.txt|3.11|2
8 | ~+|batchin=%t1%Prem%x%_%x%468.txt|3.11|2
9 | ~/TOD Period %x% done!
10 | batchin=
11 | reports=
12 | q
--------------------------------------------------------------------------------
/scripts/__init__.py:
--------------------------------------------------------------------------------
https://raw.githubusercontent.com/CMAP-REPOS/cmap_abm/e0b6bd57e34e5bd3d95ecb159feb1f4745baf0c6/scripts/__init__.py
--------------------------------------------------------------------------------
/scripts/cmapInputFileGen.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Shortest Path code for CMAP non-motorized modes
3 | # Roshan Kumar, kumarr1@pbworld.com, 12/20/2012
4 | #############################################################################
5 |
6 | #Should always be placed in the same folder as SPwrapper.py
7 |
8 | #Combines the CSVs generated into one node CSV and one link CSV. Both are stored in ...\SP\final_output
9 |
10 | import parameters
11 | import csv
12 |
13 | #############################################################################
14 |
15 | def createNetworkInput():
16 | linkCSV = open(parameters.linkFile,"w")
17 | linkCSV.write('"FNODE_","TNODE_","LENGTH"\n')
18 | nodeCSV = open(parameters.nodeFile,"w")
19 | nodeCSV.write('"ID","X","Y"\n')
20 |
21 | for record in csv.DictReader(open(parameters.centroidsMAZFile, "r")):
22 | maz = int(record["MAZ_ID_9"])
23 | x = float(record["X"])
24 | y = float(record["Y"])
25 | nodeCSV.write("%d,%f,%f\n" % (maz, x, y))
26 |
27 | for record in csv.DictReader(open(parameters.tapFile, "r")):
28 | tap = int(record["tap_id"])
29 | x = float(record["x"])
30 | y = float(record["y"])
31 | nodeCSV.write("%d,%f,%f\n" % (tap, x, y))
32 |
33 | for record in csv.DictReader(open(parameters.nodeFile_int, "r")):
34 | nid = int(record["ID"])
35 | x = float(record["X"])
36 | y = float(record["Y"])
37 | nodeCSV.write("%d,%f,%f\n" % (nid, x, y))
38 |
39 | for record in csv.DictReader(open(parameters.connectorFile, "r")):
40 | fnode = int(record["FNODE_"])
41 | tnode = int(record["TNODE_"])
42 | length= float(record["LENGTH"])
43 | linkCSV.write("%d,%d,%f\n" % (fnode, tnode, length))
44 |
45 | for record in csv.DictReader(open(parameters.tapConnectorFile, "r")):
46 | fnode = int(record["FNODE_"])
47 | tnode = int(record["TNODE_"])
48 | length= float(record["LENGTH"])
49 | linkCSV.write("%d,%d,%f\n" % (fnode, tnode, length))
50 |
51 | for record in csv.DictReader(open(parameters.linkFile_int, "r")):
52 | fnode = int(record["FNODE_"])
53 | tnode = int(record["TNODE_"])
54 | length= float(record["LENGTH"])
55 | linkCSV.write("%d,%d,%f\n" % (fnode, tnode, length))
56 |
57 | linkCSV.close()
58 | nodeCSV.close()
59 |
60 | #############################################################################
61 |
--------------------------------------------------------------------------------
/scripts/cmapMAZLUShare.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Shortest Path code for CMAP non-motorized modes
3 | # Roshan Kumar, kumarr1@pbworld.com, 12/20/2012
4 | #############################################################################
5 |
6 | #Should always be placed in the same folder as SPwrapper.py
7 |
8 | import math
9 | #import random
10 | import sys
11 | #import pickle
12 | import csv
13 | import os
14 | #import pdb
15 | from collections import defaultdict
16 | from heapq import heapify, heappop, heappush
17 | from array import *
18 | #import numpy as np
19 | import shapefile
20 | import parameters
21 |
22 | #############################################################################
23 |
24 | def getLandUseShare():
25 |
26 | result = {}
27 | #count1 = {}
28 | for i in range(16820):
29 | result[i] = []
30 | #count1[i] = 0
31 | fileName = parameters.mazLUFile
32 | line = 0
33 | for fields in csv.DictReader(open(fileName, "r")):
34 | line += 1
35 | mazId = int(fields["MAZ_ID"])
36 | lu1 = float(fields["LU_1"])
37 | lu2 = float(fields["LU_2"])
38 | lu3 = float(fields["LU_3"])
39 | lu4 = float(fields["LU_4"])
40 | lu5 = float(fields["LU_5"])
41 | lu6 = float(fields["LU_6"])
42 | lu7 = float(fields["LU_7"])
43 | lu8 = float(fields["LU_8"])
44 | lu9 = float(fields["LU_9"])
45 | lu10 = float(fields["LU_10"])
46 | lu11 = float(fields["LU_11"])
47 | lu12 = float(fields["LU_12"])
48 |
49 |
50 | if mazId == 0:
51 | continue
52 |
53 | result[mazId].append(lu1)
54 | result[mazId].append(lu2)
55 | result[mazId].append(lu3)
56 | result[mazId].append(lu4)
57 | result[mazId].append(lu5)
58 | result[mazId].append(lu6)
59 | result[mazId].append(lu7)
60 | result[mazId].append(lu8)
61 | result[mazId].append(lu9)
62 | result[mazId].append(lu10)
63 | result[mazId].append(lu11)
64 | result[mazId].append(lu12)
65 |
66 | #count1[mazId] = count1[mazId]+1
67 |
68 | result = dict([(k,v) for k,v in result.items() if len(v)>0])
69 |
70 | #print count1[2504]
71 | #print result[2504][1]
72 |
73 | return result
74 |
75 | #############################################################################
76 |
77 | def getMAZaroundNodes():
78 |
79 | result = {}
80 | #count1 = {}
81 | #for i in range(95141329):
82 | #result[i] = []
83 | #count1[i] = 0
84 | fileName = parameters.nodesinMAZFile
85 | line = 0
86 | for fields in csv.DictReader(open(fileName, "r")):
87 |
88 | mazId = int(fields["MAZ_ID_9"])
89 | nodes = int(fields["Nodes"])
90 | result[nodes] = []
91 |
92 | for fields in csv.DictReader(open(fileName, "r")):
93 |
94 | mazId = int(fields["MAZ_ID_9"])
95 | nodes = int(fields["Nodes"])
96 |
97 |
98 |
99 | if mazId == 0:
100 | continue
101 |
102 | result[nodes].append(mazId)
103 | #count1[mazId] = count1[mazId]+1
104 |
105 | result = dict([(k,v) for k,v in result.items() if len(v)>0])
106 | #print result[904259084][1]
107 |
108 | return result
109 |
110 | getLandUseShare()
111 | getMAZaroundNodes()
112 |
--------------------------------------------------------------------------------
/scripts/cmapParallelPP.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Shortest Path code for CMAP non-motorized modes
3 | # Roshan Kumar, kumarr1@pbworld.com, 12/20/2012
4 | #############################################################################
5 |
6 | #Should always be placed in the same folder as SPwrapper.py
7 |
8 | #The multiple CSVs generated by the SP script are conflated into one CSV
9 |
10 | import parameters
11 | import os
12 |
13 | #############################################################################
14 |
15 | def parallelPostProcess():
16 | try:
17 | f = open(parameters.outputSPFile, "w")
18 | #f.write('"Origin","Destin","SP","Euclid","Ratio","nearHwy","nearHwyDist","LU_1","LU_2","LU_3","LU_4","LU_5","LU_6","LU_7","LU_8","LU_9","LU_10","LU_11","LU_12"\n')
19 | f.write('"Origin","Destin","SP","Euclid","Ratio"\n')
20 |
21 | processNo = 1
22 | while processNo<=parameters.number_of_cores:
23 | fileNo = "%s_" % (processNo)
24 | fileN = fileNo+"_test_parallel_SP.txt"
25 | #print fileN
26 | inputFile = os.path.join(parameters.root, parameters.int_output, fileN)
27 | #print inputFile
28 | with open(inputFile, "r") as file1:
29 | old = file1.read()
30 | f.write(old)
31 | processNo = processNo+1
32 | f.close()
33 | except IOError:
34 | pass
35 |
--------------------------------------------------------------------------------
/scripts/cmapPostProcess.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Shortest Path code for CMAP non-motorized modes
3 | # Roshan Kumar, kumarr1@pbworld.com, 12/20/2012
4 | #############################################################################
5 |
6 | #Should always be placed in the same folder as SPwrapper.py
7 |
8 | #Wrapper code that post processes the SP output
9 |
10 | import parameters
11 | import cmapParallelPP
12 | if parameters.isBaseYear:
13 | import cmapRStoSubZone
14 | else:
15 | import cmapRStoSubZone_future
16 |
17 | cmapParallelPP.parallelPostProcess()
18 | cmapRStoSubZone.rsToSubZone()
19 |
--------------------------------------------------------------------------------
/scripts/cmapRStoSubZone.py:
--------------------------------------------------------------------------------
1 | #############################################################################
2 | # Shortest Path code for CMAP non-motorized modes
3 | # Roshan Kumar, kumarr1@pbworld.com, 12/20/2012
4 | #############################################################################
5 |
6 | #Should always be placed in the same folder as SPwrapper.py
7 |
8 | #Post-processes the output to exclude all TAP-TAP paths
9 |
10 | import math
11 | #import random
12 | import sys
13 | #import pickle
14 | import csv
15 | import os
16 | #import pdb
17 | from collections import defaultdict
18 | from heapq import heapify, heappop, heappush
19 | from array import *
20 | #import numpy as np
21 | import shapefile
22 | import parameters
23 |
24 | #############################################################################
25 |
26 | def rsToSubZone():
27 | output = open(parameters.spMAZ2TAP, "w")
28 | output.write('"Origin","Destin","SP","Euclid","Ratio"\n')
29 |
30 | szForTap = {}
31 | for fields in csv.DictReader(open(parameters.tapFile, "r")):
32 | tap_id = int(fields["tap_id"])
33 | sz09 = int(fields["maz09"])
34 | szForTap[tap_id] = sz09
35 | count_ = 1
36 | for fields in csv.DictReader(open(parameters.outputSPFile, "r")):
37 | origin = int(fields["Origin"])
38 | destin = int(fields["Destin"])
39 | sp = float(fields["SP"])
40 | euclid = float(fields["Euclid"])
41 | ratio = float(fields["Ratio"])
42 | if destin > parameters.minTAPID:
43 | try:
44 | if origin == szForTap[destin]:
45 | sp = euclid
46 | count_ = count_+1
47 | except KeyError:
48 | pass
49 |
50 |
51 |
52 | if (origin+destin parameters.minTAPID:
67 | try:
68 | if origin == szForTap[destin]:
69 | sp = euclid
70 | count_ = count_+1
71 | except KeyError:
72 | pass
73 |
74 | if (origin