Do not copy paste: AF

want view factors from all the surfaces $.$ combine my two codes so I can have view factors from all surfaces

Using your Monte Carlo code, determine the view factor between a control surface on the left side of

the V$-$groove with every control surface on the opposite side of the V$-$groove. Do this for every control

surface on the left side of the V$-$groove. This means you should know the view factor between any two

control surfaces on opposing sides of the V$-$groove.

combine the codes below

import random

import numpy as np

import math

import matplotlib.pyplot as plt

# Numerical solution

def rotate$\_$vector$($vector$,$ angle$)$:

x$,$ y $=$ vector

theta $=$ math.radians$($angle$)$

x$\_$rotated $=$ math.cos$($theta$)*$ x $-$ math.sin$($theta$)*$ y

y$\_$rotated $=$ math.sin$($theta$)*$ x $+$ math.cos$($theta$)*$ y

return $[$x$\_$rotated, y$\_$rotated$]$

def find$\_$intersection$($x$,$ y$,$ xhat, yhat, surface$1,$ surface$2)$:

xi$\_$s $=$ surface$1[\text{'}$X$0\text{'}]+$ surface$1[\text{'}$XHAT$\text{'}]*$ x

yi$\_$s $=$ surface$1[\text{'}$Y$0\text{'}]+$ surface$1[\text{'}$YHAT$\text{'}]*$ y

xi$\_$r $=$ surface$2[\text{'}$X$0\text{'}]+$ surface$2[\text{'}$XHAT$\text{'}]*$ x

yi$\_$r $=$ surface$2[\text{'}$Y$0\text{'}]+$ surface$2[\text{'}$YHAT$\text{'}]*$ y

return xi$\_$s$,$ yi$\_$s$,$ xi$\_$r$,$ yi$\_$r

def monte$\_$carlo$\_$ray$\_$tracing$($num$\_$rays$)$:

# Step $1$: Define surfaces

surface$\_$length $=1$ # Length of the surfaces

surface$\_$distance $=1$ # Distance between the surfaces

surface$1=\{\text{'}$X$0\text{'}$: $0,\text{'}$Y$0\text{'}$: $0,$ 'XHAT': surface$\_$length, 'YHAT': $0\}$

surface$2=\{\text{'}$X$0\text{'}$: $0,\text{'}$Y$0\text{'}$: surface$\_$distance, 'XHAT': surface$\_$length, 'YHAT': $0\}$

rays$\_$surface$1=[]$

for $\_$ in range$($num$\_$rays$)$:

# Step $3$: Find a random spot for each ray

ts $=$ random.random$()$

tr $=$ random.random$()$

# Step $4$: Assign emission locations

x$\_$ray$\_$s$1,$ y$\_$ray$\_$s$1=$ surface$1[\text{'}$X$0\text{'}]+$ surface$1[\text{'}$XHAT$\text{'}]*$ ts$,$ surface$1[\text{'}$Y$0\text{'}]+$ surface$1[\text{'}$YHAT$\text{'}]*$ ts

surface$1\_$normal $=$ rotate$\_$vector$([$surface$1[\text{'}$XHAT$\text{'}],$ surface$1[\text{'}$YHAT$\text{'}]],90)$

R$\_$theta$\_$s$1,$ R$\_$phi$\_$s$1=$ random.random$(),$ random.random$()$

theta$\_$s$1=$ math.asin$($R$\_$theta$\_$s$1)$

if R$\_$phi$\_$s$1<0\mathrm{.}5$:

theta$\_$s$1=-$theta$\_$s$1$

xhat$\_$ray$\_$s$1=$ math.cos$($theta$\_$s$1)*$ surface$1\_$normal$[0]-$ math.sin$($theta$\_$s$1)*$ surface$1\_$normal$[1]$

yhat$\_$ray$\_$s$1=$ math.sin$($theta$\_$s$1)*$ surface$1\_$normal$[0]+$ math.cos$($theta$\_$s$1)*$ surface$1\_$normal$[1]$

rays$\_$surface$1.$append$(\{\text{'}$X$\text{'}$: x$\_$ray$\_$s$1,\text{'}$Y$\text{'}$: y$\_$ray$\_$s$1,$ 'XHAT': xhat$\_$ray$\_$s$1,$ 'YHAT': yhat$\_$ray$\_$s$1,$ 'Reflections': $0,\text{'}$ts$\text{'}$: None, $\text{'}$tr$\text{'}$: None$\})$

# Step $6$: Find intersection points and save ts and tr for each surface

for ray in rays$\_$surface$1$:

xi$\_$s$,$ yi$\_$s$,$ xi$\_$r$,$ yi$\_$r $=$ find$\_$intersection$($ray$[\text{'}$X$\text{'}],$ ray$[\text{'}$Y$\text{'}],$ ray$[\text{'}$XHAT$\text{'}],$ ray$[\text{'}$YHAT$\text{'}],$ surface$1,$ surface$2)$

denominator $=$ ray$[\text{'}$YHAT$\text{'}]*$ surface$2[\text{'}$XHAT$\text{'}]-$ ray$[\text{'}$XHAT$\text{'}]*$ surface$2[\text{'}$YHAT$\text{'}]$

if denominator $!=0$:

tr $=($xi$\_$s $*$ surface$2[\text{'}$YHAT$\text{'}]-$ yi$\_$s $*$ surface$2[\text{'}$XHAT$\text{'}]+$ surface$2[\text{'}$X$0\text{'}]*$ surface$2[\text{'}$YHAT$\text{'}]-$ surface$2[\text{'}$Y$0\text{'}]*$ surface$2[\text{'}$XHAT$\text{'}])/$ denominator

ts $=($ray$[\text{'}$X$\text{'}]+$ ray$[\text{'}$XHAT$\text{'}]*$ tr $-$ surface$1[\text{'}$X$0\text{'}])/$ surface$1[\text{'}$XHAT$\text{'}]$

ray$[\text{'}$ts$\text{'}]=$ ts

ray$[\text{'}$tr$\text{'}]=$ tr

# Step $7$: Determine impact surface

for ray in rays$\_$surface$1$:

ts $=$ ray$[\text{'}$ts$\text{'}]$

tr $=$ ray$[\text{'}$tr$\text{'}]$

if $0<=$ ts $<=1$:

if ray$[\text{'}$Reflections$\text{'}]==0$ or ts $<$ rays$\_$surface$1[$ray$[\text{'}$Reflections$\text{'}]-1][\text{'}$ts$\text{'}]$:

ray$[\text{'}$ImpactedSurface$\text{'}]=1$

else:

ray$[\text{'}$ImpactedSurface$\text{'}]=2$

else:

ray$[\text{'}$ImpactedSurface$\text{'}]=0$

# Print the results

hit$\_$count $=$ sum$($ray$[\text{'}$ImpactedSurface$\text{'}]==1$ for ray in rays$\_$surface$1)$

miss$\_$count $=$ num$\_$rays $-$ hit$\_$count

#print$($f$"$Number of rays that hit the surface: $\{$hit$\_$count$\}")$

#print$($f$"$Number of rays that did not hit the surface: $\{$miss$\_$count$\}")$

# view factor

view$\_$factor $=$ hit$\_$count $/$ num$\_$rays

print$($f$"$View Factor: $\{$view$\_$factor$\}")$

# Example usage

monte$\_$carlo$\_$ray$\_$tracing$(100000)$

second code

import numpy as np

import matplotlib.pyplot as plt

# Problem Definitions

phi $=120$ # phi $($deg$)$

phi$\_$radians $=$ phi $*$ np$.$pi $/180$ # phi $($rad$)$

Nsurf $=10$ # number of surfaces PER SIDE

length $=0\mathrm{.}1$ # total length per side $($m$)$

dL $=$ length $/$ Nsurf # length of one surface

# Set Initial Conditions

X$0=0$

Y$0=0$

# Calculate Surfaces for First Side

for i in range$($Nsurf$)$:

XHAT $=$ X$0+$ dL $*$ np$.$sin$($phi$\_$radians $/2)$

YHAT $=$ Y$0-$ dL $*$ np$.$cos$($phi$\_$radians $/2)$ # MINUS for $-$y

temp $=$ "surface" $+$ str$($i$)+"=\{\text{'}$X$0\text{'}$:X$0,\text{'}$Y$0\text{'}$:Y$0,$'XHAT':XHAT,'YHAT':YHAT$\}"$

exec$($temp$)$

X$0=$ XHAT

Y$0=$ YHAT

# Calculate Surfaces for Second Side

for i in range$($Nsurf$)$:

XHAT $=$ X$0+$ dL $*$ np$.$sin$($phi$\_$radians $/2)$

YHAT $=$ Y$0+$ dL $*$ np$.$cos$($phi$\_$radians $/2)$ # PLUS for $+$Y

temp $=$ "surface" $+$ str$($i $+$ Nsurf$)+"=\{\text{'}$X$0\text{'}$:X$0,\text{'}$Y$0\text{'}$:Y$0,$'XHAT':XHAT,'YHAT':YHAT$\}"$

exec$($temp$)$

X$0=$ XHAT

Y$0=$ YHAT

# Plot Surfaces

for i in range$(2*$ Nsurf$)$:

temp $="$plt$.$plot$([$surface$"+$ str$($i$)+"[\text{'}$X$0\text{'}],$ surface" $+$ str$($i$)+"[\text{'}$XHAT$\text{'}]],[$surface$"+$ str$($i$)+"[\text{'}$Y$0\text{'}],$ surface" $+$ str$($i$)+"[\text{'}$YHAT$\text{'}]],$ marker$=\text{'}$o$\text{'})"$

exec$($temp$)$

plt$.$show$()$

the output should have similar values or I'll down vote

F $0->10=0\mathrm{.}000862$

F $0->11=0\mathrm{}$