How to use a code to get the corresponding value of a loop in another file

I have a particular challenge in a project I am trying to round up. I have two files with loops on both files which I need to interact correspondingly that is to say both files contain the same loops. As such what I am looking is a case where the primary (main) file start with a value of 1 goes to the other file picks the value of one as well run the process have the value comes back thru the loop pick the next value of 2 goes all over again to the secondary file pick 2 comes back till the end of the loop.

The second challenges I am having is the identification of files generated by the second file as that is the information needed by the first file

% This is the first file
function D_fat = scr_response(x)
  ;;;;;;;;
  ;;;;;;;;
  environment.WaveDirection = 0.0;
  for environment_wavedirection = 0:45:360
    if environment_wavedirection < 360
      environment.WaveDirection = environment_wavedirection;
    end 
    ;;;;;;;;;;;;;;
    ;;;;;;;;;;;;;;
    for scores = [caseload01(x),caseload02(x),caseload03(x)]
      ;;;;;
      ;;;;;
      D_fat = xlsread('filename.xlsx',1,'B8')
      D_total = sum(D_fat)
    end
  end
end

% This is the second file

function seastateload1 = case01(x)
  ;;;;;;;
  ;;;;;;;
  ;;;;;;;
  for environment_wavedirection = 0:45:360
    if environment_wavedirection < 360
      environment.WaveDirection = environment_wavedirection;
    end 
    model.CalculateStatics; 
    model.RunSimulation;                    %Run the Dynamic simulation
    model.SaveSimulation(sprintf('loadcase1Hs%d.sim',environment.WaveDirection));       
    %I need the name of the file generated as that is what is needed in file 1
    seastateload1 = 'generated file name'   
  end 
end

Answer 1

[dpb]

Your outer loop range is 0:45:360 but you only set a value for the struct environment for those values less than 360. What is it supposed to do at 360? As coded, it will rerun with the 315 value left over from the previous iteration. Looks like you should probably just quit? Or adjust the upper limit of loop to 360-45.

For the filename, just save the name you created as a variable instead of burying it inside the argument of the .SaveSimulation method--

function seastateload1 = case01(x)
    ...
    ...
    seastateload1=sprintf('loadcase1Hs%d.sim',environment.WaveDirection);
    model.SaveSimulation(sprintf(seastateload1 ,environment.WaveDirection));       
  end 
end

[joshua Abam]

Thanks for your response what i intend to achieve by the range is to capture 0,45,90,135,180,225,270 and 315 while at 360 it will quit. I had to add the if statment to the for loop. to ensure it does'nt capture the 360. although the main reason for keeping the 360 is that if i consider using difference of 20 or 30 that will enable the code capture up to 340 or 330 as the case may be.

[dpb]

The reason for not burying magic numbers in the code..use variables instead.

The loop as written, however, does NOT quit but runs for environment_wavedirection = 360 the last loop but all the if clause does is to leave environment.WaveDirection at the previously set value. So it runs the last case twice.

To fix that, use something like

function D_fat = scr_response(x,dTheta)
  ;;;;;;;;
  ;;;;;;;;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;
    ...

You could make dTheta optional argument and have a default of 45 or whatever, but this shows the idea. Don't mix up potentially changeable parameters into the code itself...then have to edit the code to make changes instead of changing just inputs.

[joshua Abam]

Hi DPD

I appreaciate your prompt responce. I have one more challenge which is to creat the link that would enable the main file to communicate and extract the information from the secondary file at a particular dTheta. let say I run the scr_response and it start by running thru the the first for loop where it starts with environment.WaveDirection = zero then go ahead to the next for loop where it calls and run the caseload01(x) at environment.WaveDirection = zero then it goes all over again but this time at environment.WaveDirection = 45 and so on till it comes to the end.

 % This is the first file
function D_fat = scr_response(x)
  ;;;;;;;;
  ;;;;;;;;
  dTheta = 45;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;  
    ;;;;;;;;;;;;;;
    ;;;;;;;;;;;;;;
    for scores = [caseload01(x),caseload02(x),caseload03(x);('SCRmodel'),...
           ('SCRmodel'),('SCRmodel');0,0,0]
      LoadCaseFileName = scores(1,:)
      LoadCaseLineName = scores(2,:)
      PeriodFrom = scores(3,:)
    end
    ;;;;;
    ;;;;;
    D_fat = xlsread('filename.xlsx',1,'B8')
    D_total = sum(D_fat)
  end
end

% This is the second file 

function seastateload1 = case01(x)
  ;;;;;;;
  ;;;;;;;
  ;;;;;;;

  dTheta = 45;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;
    model.CalculateStatics; 
    model.RunSimulation;                    %Run the Dynamic simulation
    seastateload1=sprintf('loadcase1Hs%d.sim',environment.WaveDirection) model.SaveSimulation(sprintf(seastateload1 ,environment.WaveDirection)) 
  end 
end

[dpb]

I can't follow what it is, precisely, you're driving at.

Your scr_response function never calls the other function.

I don't see if you want data by the theta angle why you are using the second function anyway, why not just do what needs doing for each theta in the first function? Then you have access to anything/everything needed.

Most particularly, why loop over theta in the second function when you're already looping over it in the first? Do what want/need for each theta then go on to the next.

[joshua Abam]

Dear Dpb

Sorry I made a mistake in the code. The entire idea is for the scr_response(x) file to run at the same dtetha as the caseload files. that means when the scr_response is called it runs at environment.waveDirection = 0 so also the caseload files are meant to run at environment.waveDirection = 0 thereafter it procceds to the next value on the loop which is environment.waveDirection = 45 till it gets to the end of the loop.

I hope I was able to explain it better this time.

Thanks

% This is the main file
function D_fat = scr_response(x)
  dTheta = 45;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;  
    for scores = [caseload01(x),caseload02(x),caseload03(x);('SCRmodel'),...
           ('SCRmodel'),('SCRmodel');0,0,0]
      LoadCaseFileName = scores(1,:)
      LoadCaseLineName = scores(2,:)
      PeriodFrom = scores(3,:)
    end
    D_fat = xlsread('filename.xlsx',1,'B8')
    D_total = sum(D_fat)
  end
end

% This is the caseload01 file 
function seastateload1 = caseload01(x)
  dTheta = 45;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;
    model.CalculateStatics; 
    model.RunSimulation;                    %Run the Dynamic simulation
    seastateload1=sprintf('loadcase1Hs%d.sim',environment.WaveDirection) 
    model.SaveSimulation(sprintf(seastateload1 ,environment.WaveDirection)) 
  end 
end

%This is the caseload02 file
function seastateload1 = caseload02(x)
  dTheta = 45;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;
    model.CalculateStatics; 
    model.RunSimulation;                    %Run the Dynamic simulation
    seastateload1=sprintf('loadcase1Hs%d.sim',environment.WaveDirection) 
    model.SaveSimulation(sprintf(seastateload1 ,environment.WaveDirection)) 
  end 
end

%This is the caseload03 file
function seastateload1 = caseload03(x)
  dTheta = 45;
  MaxTheta=360-dTheta;
  for environment_wavedirection = 0:dTheta:MaxTheta
    environment.WaveDirection = environment_wavedirection;
    model.CalculateStatics; 
    model.RunSimulation;                    %Run the Dynamic simulation
    seastateload1=sprintf('loadcase1Hs%d.sim',environment.WaveDirection) 
    model.SaveSimulation(sprintf(seastateload1 ,environment.WaveDirection)) 
  end 
end

[dpb]

I fail to see the point of all this yet...all three caseloadNN functions are identical and simply

for environment_wavedirection = 0:dTheta:MaxTheta
    model.CalculateStatics; 
    model.RunSimulation;                    %Run the Dynamic simulation
    LoadCaseFileName = sprintf('loadcase1Hs%d.sim',environment_wavedirection) 
    model.SaveSimulation(sprintf(seastateload1 ,environment.WaveDirect
    D_fat = xlsread('filename.xlsx',1,'B8')
    D_total = sum(D_fat)
end

does the basic thing...excepting for what you call LoadCaseLineName and PeriodFrom which afaict will simply return the string as LoadCaseFileName because the functions are identical.

At the end of each pass through the loop you've got the simulation run for the given theta; do whatever it is needed with it for that theta value before ending the loop and then the next pass will be the next theta incremented by dTheta.

I really don't know what you think

    for scores = [caseload01(x),caseload02(x),caseload03(x);('SCRmodel'),...
           ('SCRmodel'),('SCRmodel');0,0,0]

is doing...you can't catenate the different level strings and the row of doubles in an array anyway...it would have to be cell array.

And, as a side note, the dummy argument variable x is never used...what's the point in having it?

[joshua Abam]

Dear Dpd

I have take another look at the whole script and have decided to make environmental.waveDirection global as such what ever is used in the main file is equally used in the caseload files.

thanks for all your help, however one last question I am trying to xlsread data from a saved sprintf line saved in excel and it is prooven difficult this is the line.

%The Idea is for the xlsread to read and extract data from each of the xlsx files generated by %the sprintf

Calculate(sprintf('fatigue_analysis%d.xlsx',environment_wavedirection))

D_fat = xlsread('filename.xlsx',1,'B8')

Thanks Kind regards

Joshua Abam

[dpb]

"...each of the xlsx files generated by %the sprintf"

sprintf does nothing of the sort. All it does is return the character array built by the arguments given to it into a variable. Something else will have to use that filename to create the actual file... xlswrite (deprecated) or writetable would be candidates.

Perhaps that's buried in whatever Calculate is, but certainly as written the above makes no sense at all.

Meanwhile in

D_fat = xlsread('filename.xlsx',1,'B8')

you use the literal string 'filename.xlsx' so whatever string variable you created for a file name is totally immaterial.

Overall, I still fail to see why you think you need the auxilary function(s) at all...plural particularly since in the code you've shown so far they're all identical. And even if they're supposed to be operating on different data, that should be part of the argument list passed to them, not have a separate function for every data set.

[joshua Abam]

Dear Dpb,

I have tried using the filename to represent the output of the sprintf files however it came back with an error message stating 'too many input'

do you pls have any other suggestion as xlswrite did not produced the required information as what is needed

filename = Calculate(sprintf('analysis%d.xlsx',environment_wavedirection));

D_fat = xlsread(filename,1,'B8')

[dpb]

What is this Calculate function? You're passing a filename to it; what it returns is anybody's guess without seeing the code.

[joshua Abam]

Dear Dpd

Thanks for all your help it was indeed appreciated.

Thanks

kind regards

Joshua Abam

[joshua Abam]

Dear Dpd

I have one more challenge

I wish to know how to sum all the output generated from the for loop. in these instance the generate results from D_fat

Thanks

function D_fat = response(x)

;;;;;;;;

;;;;;;;;

;;;;;;

;;;;;;;;

dTheta = 45

MaxTheta = 360-dTheta

D_total = 0.0;

for strength = 0:dTheta:MaxTheta

WaveDirection = strength;

;;;;;;;;;;;;;;

;;;;;;;;;;;;;;

;;;;;

;;;;;

D_fat = xlsread('filename.xlsx',1,'B8')

%The idea is to sum all the generated D_fat from the for loop above

D_total = sum(D_fat)

[dpb]

Just keep a running sum...