Dakota File - Direct Call

[rezvan]

Hello,
First, Thank you so much for the extremely helpful Pelicun module.
I have a few questions and I will be so thankful if you could help with them.

1- As mentioned in workshop#3, in case the Pelicun module is being used directly, it is possible to delete unnecessary parts of the dakota.json file. The file 6_story_BRBF_pelicun.json file prepared for the workshop was an example of the shortened file. I tried to run Pelicun using this file and the response.csv file prepared for the workshop, but it seems there is some missing information in dakota file. I attached the log file and the message on the command prompt, I was wondering if you could share a shortened but correct json file.

2- For performing analysis on a large number of buildings, how is it possible to generate the dakota.json file automatically in python. For example, reading the required information from a csv file and generating the dakota file for each building automatically. As I understood, the information could be changed manually from the file or using PBE software and get dakota file which neither is practical for a large number of cases.

Thank you so much and again I appreciate your precious effort on this helpful and informative module.


P.S. Here is the message after running the DL_calculation.py in command prompt:


C:\Users\pooya\Desktop\Temp\pelicun>python C:\PBE\applications\performDL\pelicun\DL_calculation.py --filenameDL 6_story_BRBF_pelicun.json  --filenameEDP response.csv
2022-02-09T05:16:02Z First line of DL_calculation
2022-02-09T05:16:03Z Initializing pelicun calculation...
Traceback (most recent call last):
  File "C:\PBE\applications\performDL\pelicun\DL_calculation.py", line 316, in <module>
    main(sys.argv[1:])
  File "C:\PBE\applications\performDL\pelicun\DL_calculation.py", line 300, in main
    run_pelicun(
  File "C:\PBE\applications\performDL\pelicun\DL_calculation.py", line 257, in run_pelicun
    A.read_inputs(DL_input_path, EDP_files[s_i], verbose=False) # make DL inputs into array of all BIM files
  File "C:\PBE\applications\performDL\pelicun\pelicunPBE\control.py", line 767, in read_inputs
    super(FEMA_P58_Assessment, self).read_inputs(path_DL_input,
  File "C:\PBE\applications\performDL\pelicun\pelicunPBE\control.py", line 165, in read_inputs
    self._AIM_in = read_SimCenter_DL_input(
  File "C:\PBE\applications\performDL\pelicun\pelicunPBE\file_io.py", line 341, in read_SimCenter_DL_input
    raise ValueError(
ValueError: PlanArea has to be specified in the DL input file to estimate injuries decision variable(s).

[adamzs]

Hi Rezvan,

Thank you for reaching out to us, it is good to hear that you find Pelicun helpful for your work.

Based on the error message, I assume the problem with the shortened file in your first question is that the names of certain attributes under the GeneralInformation part of the file do not follow the standard naming convention we introduced across our tools last summer. I assume the workshop you refer to was held before that date. Since those standard names were introduced, pelicun 2.6 and later versions look for NumberOfStories and PlanArea under GeneralInformation and does not accept other versions of these attributes anymore. Please take a look at the shortened file and edit it if needed so that it follows these conventions. That should help you get past the error.

Several approaches are available to analyze a large number of buildings:
- If this is a regional analysis, i.e., the buildings are in a geographical context with a location assigned to each, then I suggest using our R2D Tool or the rWHALE backend to run the analysis.
- If this is more of a parametric study on a large set of archetypes, then you can do one of the following:
  = As you mentioned, you could prepare a dakota.json file for each building. I know grad students who do this through MatLab (without using PBE at all) by printing out a text file and then running Pelicun as an application directly from MatLab. Nevertheless, I agree with you that this approach is far from efficient.
  = A better way to handle this would be to import pelicun in a Python script and use it as a library rather than as an application. PBE uses the DL_calculation.py script to run pelicun as an application. If you take a look at that script (it's under tools in the pelicun package), you'll see how pelicun is imported and how the various methods in the library are called. Calling these directly will make your code more efficient, but you'd still need to prepare input files and read output files if you are using pelicun 2.6. Also note that the dakota.json input file is just a dictionary under Python. So, you can prepare a dictionary and save it to a json file using the json package in Python.
  = One of the major changes in pelicun 3 is a redesign of how researchers can interact with the library exactly to support the use case that you have. You do not need to prepare an input file and you can get the outputs directly as Python objects to stay within Python for the entire analysis and make large calculations much more efficient. I will present these features next Friday (Feb 18) during our Live Expert Tips session. If you are interested, I encourage you to register and participate in the event. Here is a link: https://designsafe-ci.zoom.us/meeting/register/tJYpdOGuqTgrHt3FR0yM7dxmCYf6kiEx5Btm

Let me know if you have further questions.

Adam

PS. I hope you don't mind if we delete the copy of your question that was posted on the PBE board.

[rezvan]

Dear Adam,

I appreciate your detailed answer and proposed solutions.
I will try all the mentioned points and I am sure I will be able to handle it with your guidance.
By the way, I am extremely grateful for the invitation, without any doubt I will participate.

Pooya

[xub]

Hello,
First, Thank you so much for the extremely helpful Pelicun module.
I have a few questions and I will be so thankful if you could help with them.

File "D:\python\lib\site-packages\numexpr\necompiler.py", line 281, in stringToExpression
    raise ValueError(f'Expression {s} has forbidden control characters.')
ValueError: Expression ((index == b'B1035.001') | (index == b'B2022.001') | (index == b'C1011.001b') | (index == b'C1011.001c') | (index == b'C3011.001b') | (index == b'C3011.002c') | (index == b'C3027.002') | (index == b'C3032.003a') | (index == b'C3032.003b') | (index == b'C3034.002') | (index == b'D1014.021') | (index == b'D3041.011c') | (index == b'D3067.012b') | (index == b'D4011.053a')) has forbidden control characters.
Closing remaining open files:D:\python\lib\site-packages\pelicun/resources/FEMA_P58_2nd_ed.hdf...done