extUtil.h:37:12: fatal error: GL/glut_cgx.h: No such file or directory
#include <GL/glut_cgx.h>
если в файле extUtil.h поменять #include <GL/glut_cgx.h> на #include «GL/glut_cgx.h» ошибка пропадет.
Вопрос как переписать
%.o: %.c
g++ -c -MD $<
#include <GL/glut_cgx.h>
#include "GL/glut_cgx.h"
g++ userFunction.o AsplitA.o AsplitL.o adjustMidsideNode.o compareStrings.o XFunktions.o badelems.o bodyMesh2.o bsort.o bsortf.o bsorti.o calcNormalen.o calcPrinc.o calcPvector.o calcWeight.o checkIfNumber.o cgx.o compare.o copyEntity.o corrad.o dataGeo.o dataMesh.o defineEntity.o dispLists.o elemChecker.o extFunktions.o extGL.o fillBody2.o frecord.o foamFaces.o fuss.o gl3grades.o getGeoDataTria.o graph.o ifind.o iinsert.o improveBadTr3.o improveMesh.o iniMeshData.o intpol.o intpol2.o intpol3.o iremove.o linelength.o m_copy.o m_sub.o mergEntity.o mesh2d.o meshSet.o meshSurf.o messages.o near3d.o nurbl2seq.o orient.o p_angle.o parser.o pickFunktions.o plotFunktions.o readEdges.o readfrd.o readFoam.o readStdCmap.o readStl.o readccx.o readDuns.o readNastran.o readNG.o readIsaac.o readstep.o readWf.o rectcyl.o renumberfrd.o selectDisplayFaces.o sendMpc.o sendSet.o setFunktions.o shapeFunctions.o spline.o splitElementsToTets.o stof.o stoi.o stopClock.o stos.o strfind.o strsplt.o sins.o surfMesh2.o sword.o trackball.o v_distA.o v_rec2cyl.o v_add.o v_angle.o v_angle_ref.o v_betrag.o v_matmult.o v_norm.o v_prod.o v_result.o v_scal.o v_rot.o v_sprod.o m_prod.o m_prodtr.o write2aba.o write2ansys.o write2aster.o write2darwin.o write2dolfyn.o write2isaac.o write2duns.o write2frd.o write2foam.o write2nas.o write2samcef.o write2tochnog.o writefbd.o writebp.o contact.o makeTriFromElems.o ../../glut-3.5/src/layerutil.o ../../glut-3.5/src/glut_dials.o ../../glut-3.5/src/glut_ext.o ../../glut-3.5/src/glut_mesa.o ../../glut-3.5/src/glut_menu.o ../../glut-3.5/src/glut_modifier.o ../../glut-3.5/src/glut_space.o ../../glut-3.5/src/glut_shapes.o ../../glut-3.5/src/glut_tablet.o ../../glut-3.5/src/glut_warp.o ../../glut-3.5/src/glut_8x13.o ../../glut-3.5/src/glut_9x15.o ../../glut-3.5/src/glut_hel10.o ../../glut-3.5/src/glut_hel18.o ../../glut-3.5/src/glut_hel12.o ../../glut-3.5/src/glut_tr10.o ../../glut-3.5/src/glut_tr24.o ../../glut-3.5/src/glut_bitmap.o ../../glut-3.5/src/glut_cursor.o ../../glut-3.5/src/glut_event.o ../../glut-3.5/src/glut_get.o ../../glut-3.5/src/glut_overlay.o ../../glut-3.5/src/glut_input.o ../../glut-3.5/src/glut_stroke.o ../../glut-3.5/src/glut_teapot.o ../../glut-3.5/src/glut_winmisc.o ../../glut-3.5/src/glut_fullscrn.o ../../glut-3.5/src/glut_cindex.o ../../glut-3.5/src/glut_roman.o ../../glut-3.5/src/glut_mroman.o ../../glut-3.5/src/glut_swidth.o ../../glut-3.5/src/glut_bwidth.o ../../glut-3.5/src/glut_win.o ../../glut-3.5/src/glut_init.o ../../glut-3.5/src/glut_util.o uselibSNL.cpp generateTet.cpp ../../libSNL/src/*.cpp -O3 -Wall -I./ -I/usr/include -I/usr/include/GL -I../../libSNL/src -I../../glut-3.5/src -I/usr/include/X11 -lGL -lGLU -lxcb -lX11 -lICE -lSM -lXt -lXext -lXfixes -lm -lpthread -lrt -o cgx
Makefile:69: recipe for target 'cgx' failed
А в терминале творится такое, мама не горюй, сюда не влезет.
$ sudo make -f ccx_Makefile > temp.txt
/usr/bin/ld: ccx_2.11.a(dgesv.o): undefined reference to symbol 'xerbla_'
//usr/lib/x86_64-linux-gnu/liblapack.so.3: error adding symbols: DSO missing from command line
collect2: error: ld returned 1 exit status
make: *** [ccx_2.11] Error 1
Как это можно обрулить?
cc -Wall -O3 -fopenmp -I ../../../spooles-master -DARCH="Linux" -DSPOOLES -DARPACK -DMATRIXSTORAGE -DNETWORKOUT -c arpack.c
In file included from arpack.c:26:0:
spooles.h:26:10: fatal error: misc.h: No such file or directory
#include <misc.h> В КАКОМ ФАЙДЕ ОН УПОМИНАЕТСЯ arpack.c или ??
^~~~~~~~
compilation terminated.
Makefile:9: recipe for target 'arpack.o' failed
make: *** [arpack.o] Error 1
CFLAGS = -Wall -O3 -fopenmp -I ../../../SPOOLES.2.2 ;; у меня установлен libspooles2.2
-DARCH="Linux" -DSPOOLES -DARPACK -DMATRIXSTORAGE -DNETWORKOUT
FFLAGS = -Wall -O3 -fopenmp
CC=cc
FC=gfortran-8 -ffixed-line-length-none -fdefault-real-8 -std=legacy -fd-lines-as-comments
.c.o :
$(CC) $(CFLAGS) -c $<
.f.o :
$(FC) $(FFLAGS) -c $<
include Makefile.inc
SCCXMAIN = ccx_2.11.c
OCCXF = $(SCCXF:.f=.o)
OCCXC = $(SCCXC:.c=.o)
OCCXMAIN = $(SCCXMAIN:.c=.o)
DIR=../../../SPOOLES.2.2
LIBS = \
$(DIR)/spooles.a \
../../../ARPACK/libarpack_INTEL.a \ ;; установлен libarpack2
-lpthread -lm -lc
ccx_2.11: $(OCCXMAIN) ccx_2.11.a $(LIBS)
./date.pl; $(CC) $(CFLAGS) -c ccx_2.11.c; $(FC) -fopenmp -Wall -O3 -o $@ $(OCCXMAIN) ccx_2.11.a $(LIBS)
ccx_2.11.a: $(OCCXF) $(OCCXC)
ar vr $@ $?
PROGRAM BALANCE
! Calculates balance after interest compounded
REAL BALANCE, INTEREST, RATE
BALANCE = 1000
RATE = 0.09
INTEREST = RATE * BALANCE
BALANCE = BALANCE + INTEREST
PRINT*, 'New balance:', BALANCE
END PROGRAM BALANCE
выдает такое
gfortran-7 balance.f
balance.f:1:1:
PROGRAM BALANCE
1
Error: Non-numeric character in statement label at (1)
balance.f:1:1:
PROGRAM BALANCE
1
Error: Unclassifiable statement at (1)
balance.f:3:1:
REAL BALANCE, INTEREST, RATE
1
Error: Non-numeric character in statement label at (1)
balance.f:3:1:
REAL BALANCE, INTEREST, RATE
1
Error: Unclassifiable statement at (1)
balance.f:4:1:
BALANCE = 1000
1
Error: Non-numeric character in statement label at (1)
balance.f:5:1:
RATE = 0.09
1
Error: Non-numeric character in statement label at (1)
balance.f:6:1:
INTEREST = RATE * BALANCE
1
Error: Non-numeric character in statement label at (1)
balance.f:7:1:
BALANCE = BALANCE + INTEREST
1
Error: Non-numeric character in statement label at (1)
balance.f:8:1:
PRINT*, 'New balance:', BALANCE
1
Error: Non-numeric character in statement label at (1)
balance.f:8:1:
PRINT*, 'New balance:', BALANCE
1
Error: Unclassifiable statement at (1)
balance.f:9:1:
END PROGRAM BALANCE
1
Error: Non-numeric character in statement label at (1)
balance.f:9:1:
END PROGRAM BALANCE
1
Error: Unclassifiable statement at (1)
f951: Error: Unexpected end of file in ‘balance.f’
какой fortran установить чтоб оно работало?
Скомпилировать сам не могу. Нужно скомпилировать CalculiX с замененным одноименным файлом dflux.f файлом dflux.f с сайта https://www.researchgate.net/project/Welding-tool-for-CalculiX. CalculiX version 2.11.
!
! Heat source for Laser welds
!
! CalculiX - A 3-dimensional finite element program
! Copyright (C) 1998-2015 Guido Dhondt
!
! This program is free software; you can redistribute it and/or
! modify it under the terms of the GNU General Public License as
! published by the Free Software Foundation(version 2);
!
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program; if not, write to the Free Software
! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
!
subroutine dflux(flux,sol,kstep,kinc,time,noel,npt,coords,
& jltyp,temp,press,loadtype,area,vold,co,lakonl,konl,
& ipompc,nodempc,coefmpc,nmpc,ikmpc,ilmpc,iscale,mi)
!
! user subroutine dflux
!
!
! INPUT:
!
! sol current temperature value
! kstep step number
! kinc increment number
! time(1) current step time
! time(2) current total time
! noel element number
! npt integration point number
! coords(1..3) global coordinates of the integration point
! jltyp loading face kode:
! 1 = body flux
! 11 = face 1
! 12 = face 2
! 13 = face 3
! 14 = face 4
! 15 = face 5
! 16 = face 6
! temp currently not used
! press currently not used
! loadtype load type label
! area for surface flux: area covered by the
! integration point
! for body flux: volume covered by the
! integration point
! vold(0..4,1..nk) solution field in all nodes
! 0: temperature
! 1: displacement in global x-direction
! 2: displacement in global y-direction
! 3: displacement in global z-direction
! 4: static pressure
! co(3,1..nk) coordinates of all nodes
! 1: coordinate in global x-direction
! 2: coordinate in global y-direction
! 3: coordinate in global z-direction
! lakonl element label
! konl(1..20) nodes belonging to the element
! ipompc(1..nmpc)) ipompc(i) points to the first term of
! MPC i in field nodempc
! nodempc(1,*) node number of a MPC term
! nodempc(2,*) coordinate direction of a MPC term
! nodempc(3,*) if not 0: points towards the next term
! of the MPC in field nodempc
! if 0: MPC definition is finished
! coefmpc(*) coefficient of a MPC term
! nmpc number of MPC's
! ikmpc(1..nmpc) ordered global degrees of freedom of the MPC's
! the global degree of freedom is
! 8*(node-1)+direction of the dependent term of
! the MPC (direction = 0: temperature;
! 1-3: displacements; 4: static pressure;
! 5-7: rotations)
! ilmpc(1..nmpc) ilmpc(i) is the MPC number corresponding
! to the reference number in ikmpc(i)
! mi(1) max # of integration points per element (max
! over all elements)
! mi(2) max degree of freedomm per node (max over all
! nodes) in fields like v(0:mi(2))...
!
! OUTPUT:
!
! flux(1) magnitude of the flux
! flux(2) not used; please do NOT assign any value
! iscale determines whether the flux has to be
! scaled for increments smaller than the
! step time in static calculations
! 0: no scaling
! 1: scaling (default)
!
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
! This version of DFLUX was made by Ossama Dreibati
!Ingenieurbeüro Dreibati
!Further Information at
!https://www.researchgate.net/publication/309385851_Low_cost
!_welding_simulation_with_open_source_codes_Crash_course-_CalculiX
!
! This Subroutine is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
!/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
implicit none
!
character*8 lakonl
character*20 loadtype
!
integer kstep,kinc,noel,npt,jltyp,konl(20),ipompc(*),
& nodempc(3,*),nmpc,ikmpc(*),ilmpc(*),node,idof,id,iscale,mi(*)
!
real*8 flux(2),time(2),coords(3),sol,temp,press,
& vold(0:mi(2),*),area,co(3,*),coefmpc(*)
integer M1
real*8 Q0,RE,RI,ZE,ZI,X0,Y0,Z0,VY,AY,PIDEG,XD,YD,
& ZD,SA,CA,A1,A2,A3,AF,XL,YL,ZL,DE,DI,R0,R02,XX,YY,ZZ,TT,R2,
& XL0,YL0,ZL0,cox,coy,coz,X_,Y_,Z_,dist
!
character*50 textt,textt1
character*2 textt2
character*1 ctext
real*8 effi,power,speed, overall_dist,weld_loc,weld_dist,
& dist_x, dist_y, dist_z,weld_time, alpha, beta,gama,
& sub_dist,start_time,y_degree,z_degree,part_dist(1000),d,
& angle(3)
integer i, j,jj,ii,n,k, io_r,method, t_node,r_node,iso_10,
& ierror, node_temp
integer node_trj(1000), node_ref(1000)
C
real*8 COND,QC,QR,QF,AC, AR,V1,V2,V3,B,C,FL_MZ,RL_MZ,
& H_BW,P_B,T_R, L_R,P_D
C
C
! INITIALIZATION
io_r = 0
i = 1
j = 1
jj = 1
textt = ''
t_node = 2
r_node = 2
sub_dist = 0.0
overall_dist = 0.0
start_time = 0.025
weld_time = time(2) - start_time
if(weld_time .lt. 0.0) return
C==========WELD PARAMETER WHICH MUST BE INPUT BY THE USER
speed = 60 !WELDING SPEED
method = 1 ! 1 LASER (Goldak) 2 ARC (DOUBLE ELLIPSOIDAL)
power = 3500.0! LASER POWER or I x U
effi = 0.225 ! EFFECIENCY
C----PARAMETER FOR ARC WELDING --> DOUBLE ELLIPSOIDAL
FL_MZ = 1.2 ! Front length of the molten zone
RL_MZ = 1.7 ! Rear length of the molten zone
C
C Width and depth
C
H_BW = 1.900 ! Half of the width of the bead
P_B = 2.700 ! Penetration of the bead
C-----PARAMETER FOR LASER WELDING --> Goldak HEAT SOURCE WITH LINEAR DECREASE OF
C HEAT INPUT WITH PENETRATION DEPTH
T_R = 0.857 ! Top RADIO
L_R = 0.738 ! LOWER RADIUS
P_D = -2.25 ! PENETRATION DEPTH. NOTE THAT UPPER LIMITE OF THE HEAT SOURCE IS ZERO
C==========TRAJECTORY NODES // NUMBER of NODES = t_node
C node_trj(number) = node number of the centre of the weld trajectory in successive order
C node_trj(first node)=number of first node of the trajectory = welding start point
C node_trj(second node)=number of second node of the trajectory
C ....
C node_trj(nth node)=number of nth node of the trajectory
C node_trj(last node)=number of last node of the trajectory = welding end point
C hereafter the trajectory is a line and therefor only the start and end point are necessary
node_trj(1)=4
node_trj(2)=82
C==========END TRAJECTORY
C==========REFERENCE NODES// NUMBER of NODES = r_node
C same as for the trajectory
node_ref(1)=8
node_ref(2)=84
C==========END REFERENCE
C=========CONDITIONS FOR TRAJECTORY
if(t_node .ne. r_node) then
write(*,*) ' the number of nodes in the trajectory and
&reference groups must be the same'
call exit(201)
endif
C==========END CONDITIONS
C==========WELDED DISTANCE
C
weld_dist = speed * weld_time
C
C========== LENGTH OF WELD TRAJECTORY
do ii= 1, t_node - 1
dist_x = co(1,node_trj(ii+1)) - co(1,node_trj(ii))
dist_x = dist_x * dist_x
dist_y = co(2,node_trj(ii+1)) - co(2,node_trj(ii))
dist_y = dist_y * dist_y
dist_z = co(3,node_trj(ii+1)) - co(3,node_trj(ii))
dist_z = dist_z * dist_z
C
dist = sqrt(dist_x + dist_y + dist_z)
part_dist(ii) = dist
overall_dist = overall_dist + dist
end do
if((kstep .eq. 1).and. (kinc .eq. 1))
& write(*,*) 'length of weld trajectory:',overall_dist
C==========DECIDING WETHER THE WELD IS COMPLETLY PERFORMED
if(weld_dist .gt. overall_dist) return
dist = 0.0
C========== LOCAL COORDINATE SYSTEM (X0,Y0,Z0)
X0 = co(1,node_trj(1))
Y0 = co(2,node_trj(1))
Z0 = co(3,node_trj(1))
do ii= 1, t_node - 1
C---- Calculate the distace between two adjecent nodes
C
sub_dist = sub_dist + part_dist(ii)
C
if(sub_dist .lt. weld_dist) then
X0 = co(1,node_trj(ii+1))
Y0 = co(2,node_trj(ii+1))
Z0 = co(3,node_trj(ii+1))
C write(*,*) 'ii', ii
else
C-------gama rotation about z axis
if(co(2,node_trj(ii+1)) .eq. co(2,node_trj(ii))) then
gama = 1.570796327
else if(co(1,node_trj(ii+1)) .eq. co(1,node_trj(ii))) then
gama = 0.d0
else
gama = atan((co(2,node_trj(ii+1))-co(2,node_trj(ii)))/
& (co(1,node_trj(ii+1)) - co(1,node_trj(ii))))
gama= 1.570796327-gama
endif
C-------beta rotation about y axis
if(co(3,node_trj(ii+1)) .eq. co(3,node_trj(ii))) then
beta = 1.570796327
else if(co(1,node_trj(ii+1)) .eq. co(1,node_trj(ii))) then
beta = 0.d0
else
beta = atan((co(3,node_trj(ii+1)) - co(3,node_trj(ii)))/
& (co(1,node_trj(ii+1)) - co(1,node_trj(ii))))
beta= 1.570796327-beta
endif
C-------alpha rotation about x axis
if(co(2,node_trj(ii+1)) .eq. co(2,node_trj(ii))) then
alpha = 1.570796327
else if(co(3,node_trj(ii+1)) .eq. co(3,node_trj(ii))) then
alpha = 0.d0
else
alpha = atan((co(3,node_trj(ii+1)) - co(3,node_trj(ii)))/
& (co(2,node_trj(ii+1)) - co(2,node_trj(ii))))
alpha= 1.570796327-alpha
endif
! end do
if(abs(gama) .lt. 0.001) gama = 0.0
if(abs(beta) .lt. 0.001) beta = 0.0
if(abs(alpha) .lt. 0.001) alpha = 0.0
! gama = -gama
! beta = -beta
! alpha = -alpha
!............TRANSLATION
cox=co(1,node_trj(ii+1))-(co(1,node_trj(ii)))
coy=co(2,node_trj(ii+1))-(co(2,node_trj(ii)))
coz=co(3,node_trj(ii+1))-(co(3,node_trj(ii)))
!............ROTATION ABOUT Z
cox=(cox*cos(gama))+(coy*sin(gama))
coy=(-1*cox*sin(gama))+(coy*cos(gama))
if(coy.lt. 0d0) gama=3.1415d0 + gama
!............ROTATION ABOUT Y
! cox = (cox*cos(beta))+(coz*sin(beta))
! coz = (-1 *cox*sin(beta))+(xoz*cos(beta))
! if(coz.lt. 0d0) beta=-beta
!...........ROTATION ABOUT X
coy=(coy*cos(alpha))+(coz*sin(alpha))
coz= (-1 *coy*sin(alpha))+(coz*cos(alpha))
if(coy.ne.(co(2,node_trj(ii+1))-(co(2,node_trj(ii)))))
& alpha=-alpha
!
gama = -gama
beta = -beta
alpha = -alpha
exit
endif
end do
C=========END LOCAL COORDINATE SYSTEM (X0,Y0,Z0)
C
C=========TRANSFORMATION FROM GLOBAL TO LOCAL COORDINATE SYSTEM
C
XX = coords(1)
YY = coords(2)
ZZ = coords(3)
C
X_ = XX - X0
Y_ = YY - Y0
Z_ = ZZ - Z0
!............ROTATION ABOUT Z
XX = (X_ * cos(gama)) + (Y_ * sin(gama))
YY = (-1 * X_ * sin(gama)) + (Y_ * cos(gama))
ZZ = Z_
!............ROTATION ABOUT Y
! XX = (XX * cos(beta)) + (ZZ * sin(beta))
! ZZ = (-1 * XX * sin(beta)) + (ZZ * cos(beta))
!...........ROTATION ABOUT X
! YY = (YY * cos(alpha)) + (ZZ * sin(alpha))
! ZZ = (-1 * YY * sin(alpha)) + (ZZ * cos(alpha))
C=========END TRANSFORMATION FROM GLOBAL TO LOCAL COORDINATE SYSTEM
C
C 122 continue
C
C==========WELD LOCATION
d = 0.0
do ii = 1, t_node - 1
d = d + part_dist(ii)
if(weld_dist .lt. d) then
if(ii .eq. 1) then
d = 0.0
exit
endif
d = d - part_dist(ii)
exit
endif
end do
TT = weld_time - (d / speed)
C==========WELDING ANGLE
do ii = 1, t_node - 1
if(X0.eq.co(1,node_trj(ii))) then
if(Y0.eq.co(2,node_trj(ii))) then
if(Z0 .eq. co(3,node_trj(ii))) then
if(co(1,node_trj(ii)).eq.co(1,node_ref(ii))) then
angle(1) = 1.570796327
else
angle(1) = co(3,node_trj(ii))-co(3,node_ref(ii))
angle(1) = angle(1)/
& (co(1,node_trj(ii))-co(1,node_ref(ii)))
angle(1) = atan(angle(1))
end if
if(co(1,node_trj(ii+1)).eq.co(1,node_ref(ii+1))) then
angle(1) = 1.570796327
else
angle(2) = co(3,node_trj(ii+1))-co(3,node_ref(ii+1))
angle(2) = angle(2)/
& (co(1,node_trj(ii+1))-co(1,node_ref(ii+1)))
angle(1) = atan(angle(1))
endif
endif
endif
endif
angle(3) = part_dist(ii)-d
angle(3) = angle(3) * (angle(1)-angle(2))
angle(3) = angle(3) / part_dist(ii)
angle(3) = angle(3) + angle(2)
end do
C========== END WELDING ANGLE
C==========END WELD LOCATION
C=========LASER WELDING
if(method.eq. 1) then
C
C
C FLUX = Q0 * exp( - R^2 / R0^2 ) with
C R^2 = ( XX-X0 )^2 + ( YY-Y0-VY*T )^2
C R0 = RE - ( RE-RI )*( ZE-ZZ+Z0 )/( ZE-ZI )
C IF R0 < RI , R0 = 0. and return
C IF R0 > RE , R0 = 0. and return
C Variables
C
Q0 = power * effi * 1000.0
RE = T_R ! Top RADIO
RI = L_R ! LOWER RADIUS
ZE = ZL0 - ZL0 ! UPPER LIMIT OF THE SOURCE
ZI = P_D ! PENETRATION DEPTH
X0 = 0.0
Y0 = 0.0
Z0 = 0.0
VY = speed
AY = -angle(3) * 180 / 3.141592654
C
C Constant
C
M1 = -1
PIDEG = ATAN(1.)
PIDEG = PIDEG / 45.0
AY = AY*PIDEG
C
C Transformation of global to local coordinates
C
XD = XX - X0
YD = VY * TT
YD = YD + Y0
ZD = ZZ - Z0
C Source rotation about Y axis
C
SA = SIN(AY)
SA = - SA
CA = COS(AY)
A1 = XD * CA
A2 = ZD * SA
XL = A1 + A2
YL = YY - YD
A1 = ZD * CA
A2 = XD * SA
ZL = A1 - A2
C write(*,*)'Xl Yl Zl', XL, YL, ZL
C
C Deciding whether the point is out of source's volume
C
DE = ZL - ZE
DI = ZL - ZI
IF( DE .GT. 0.00001 ) RETURN
A1 = DI + RI
C write(*,*) 'LASER',A1
IF( A1 .LT. 0.00001 ) RETURN
C
C write(*,*) 'LASER'
C R^2 computation
C
A1 = XL * XL
A2 = YL * YL
R2 = A1 + A2
A3 = DI * DI
IF( ZL .LE. ZI ) R2 = R2 + A3 + A3
C
C R0^2 computation
C
A1 = RE - RI
A2 = ZE - ZI
A3 = ZE - ZL
R0 = A3 / A2
R0 = R0 * A1
R0 = RE - R0
IF( ZL .LE. ZI ) R0 = RI
R02 = R0 * R0
C
C F computation
C
C write(*,*) 'R',R2,R02
IF( R2 .GT. R02 ) RETURN
A1 = R2 / R02
A2 = M1 * A1
A2 = EXP( A2 )
C
C F computation
C
FLUX(1) = Q0 * A2
C
C write(*,*) XL,YL,ZL,FLUX(1)
return
end if
C===========================================END LASER
C===========================================ARC WELDING
C Standard ARC Power source
C Power source dimensions see variables below
C
if(method.eq. 2) then ! 2 for arc welding
C
C
C Coordinates of the gauss point treated and time
C
C Parameters of the Goldak power source
C
C The absorbed power is defined within an ellipsoid
C
C Definition of the maximum front and rear power intensity
C
QF = 1.0 * power * effi * 1000.0 ! Normalized maximum front power source intensity
QR = 0.833 * power * effi * 1000.0 ! Normalized maximum rear power source intensity
C
C Definition of the measures of the Goldak ellipsoid
C They should be inside the molten zone
C
AF = FL_MZ ! Front length of the molten zone
AR = RL_MZ ! Rear length of the molten zone
C
C Width and depth
C
B = H_BW ! Half of the width of the bead
C = P_B ! Penetration of the bead
C
C Position in space - completely handled by
C the welding wizzard - weldline
C
X0 = 0.000 ! X initial location of source center
Y0 = 0.000 ! Y initial location of source center
Z0 = 0.000 ! Z initial location of source center
VY = speed ! Source displacement velocity
AY = -angle(3) * 180 / 3.141592654 !0.000 Angle of torch [deg.]
C
C Computation of the absorbed power
C
C F = QC * V1 * V2 * V3 with
C V1 = exp( -( YY-Y0-VY*TT )^2/AC^2 )
C V2 = exp( -( XX-X0 )^2/B^2 )
C V3 = exp( -( ZZ-Z0 )^2/C^2 )
C if ( -YY + Y0 +VY*TT ) greater than 0
C QC = QF et AC = AF
C else
C QC = QR et AC = AR
C
C Constant
C
M1 = -1
PIDEG = ATAN(1.)
PIDEG = PIDEG / 45.
AY = AY * PIDEG
C
C Transformation of global to local coordinates
C
XD = XX - X0
YD = VY * TT
YD = YD + Y0
ZD = ZZ - Z0
C
C Source rotation about Y axis
C
SA = SIN(AY)
SA = - SA
CA = COS(AY)
A1 = XD * CA
A2 = ZD * SA
XL = A1 + A2
YL = YY - YD
A1 = ZD * CA
A2 = XD * SA
ZL = A1 - A2
C
C Condition computation, QC and AC initialisation
C
COND = VY * YL
IF (VY .EQ. 0.) COND = YL
QC = QR
AC = AR
IF(COND .GT. 0. ) QC = QF
IF(COND .GT. 0. ) AC = AF
C
C Vi computation
C
A1 = YL * YL
A2 = AC * AC
A2 = A1 / A2
A2 = M1 * A2
V1 = EXP(A2)
C
C V2 computation
C
A1 = XL * XL
A2 = B * B
A2 = A1 / A2
A2 = M1 * A2
V2 = EXP(A2)
C
C V3 computation
C
A1 = ZL * ZL
A2 = C * C
A2 = A1 / A2
A2 = M1 * A2
V3 = EXP(A2)
C
C F computation
C
flux(1) = QC * V1 * V2 * V3
C
RETURN
endif
END
C===========================================END ARC
C
Включаю макрос на запись, черчу отрезок, ставлю размер, ставлю допуски вызвав окно ДОПУСК. Итог записана шапка подключения API к Компас и всё. Как программно проставлять допуски?
Существует такой код
(defun make-test-bar-chart-drawing-object
(&key (pane-title "Example of a bar chart")
title-position
(values *example-bar-chart-values*)
(orientation :upward)
(title-color :purple)
ellipses
(colors '(:red :green :blue :black :purple))).....)
(defun make-test-bar-chart-drawing-object
(&key ((:pane-title pane-title) "Example of a bar chart")
((:title-position title-position))
((:values values) "*example-bar-chart-values*")
((:orientation orientation) ":upward")
((:title-color title-color) ":purple")
((:ellipses ellipses))
((:colors colors) "'(:red :green :blue :black :purple)"))
(values
(format nil "pane-title ~A" pane-title)
(format nil "title-position ~A" title-position)
(format nil "values ~A" values)
(format nil "orientation ~A" orientation)
(format nil "title-color ~A" title-color)
(format nil "ellipses ~A" ellipses)
(format nil "colors ~A" colors)))
(make-test-bar-chart-drawing-object
:pane-title pane-title
:title-position "title-position"
:values :values
:orientation :orientation
:title-color :title-color
:ellipses "ellipses"
:colors :colors)
"pane-title PANE-TITLE" здесь не норм
"title-position title-position"
"values VALUES" здесь не норм
"orientation ORIENTATION" здесь не норм
"title-color TITLE-COLOR" здесь не норм
"ellipses ellipses"
"colors COLORS" здесь не норм
Какова методология разработки кода для следующей задачи
нужно создать окно
поводив в окне мышкой
получаем кривую
окно сжимаем -> раскрываем, кривая не исчезла
окно сворачиваем -> разворачиваем, кривая не исчезла
не могу прочитать
(defun cached-display-draw-an-arrow (pane x y width height)
(declare (ignore x y width height))
(when-let (dragging-info (capi:output-pane-cached-display-user-info pane))
(destructuring-bind (center-x center-y end-x end-y)
dragging-info -------- >>>>>>>это место просто вышибает
(let* ((diff-x (- end-x center-x))
(diff-y (- end-y center-y))
(len-square (+ (* diff-x diff-x) (* diff-y diff-y))))
(when (> len-square 5) ;; don't try too short
(let* ((len (sqrt len-square))
(angle (atan diff-y diff-x))
(scale (/ len 100)))
(gp:with-graphics-translation (pane center-x center-y)
(gp:with-graphics-rotation (pane angle)
(gp:with-graphics-scale (pane scale 1)
(cached-display-internal-draw-an-arrow pane))))))))))
(destructuring-bind (parameter*) list
body-form*) здесь всё понятно, а в
примере приведённом выше не могу понять где vars list form
(destructuring-bind (x y z) (list 1 2 3)
(list :x x :y y :z z)) ==> (:X 1 :Y 2 :Z 3)
The form
(when-let (position (search string1 string2))
(print position))
macroexpands to
(let ((position (search string1 string2)))
(when position
(print position)))
LispWorks CAPI. Надо отрисовать линию в output-pane посредством функции make-draw-line. Иерархия классов
lw-gt:objects-displayer
capi:pinboard-layout
capi:output-pane
^
|_________________________
|
lw-gt:apply-drawing-object |
lw-gt:drawing-object |
^ |
| |
standard-object ------------------
^
|
T
Не понимаю где начинать и что начинать.
Существует 3 листа
A (1 2 3 4 5)
B (6 7 8 9 10)
C (11 -12 13 -14 -15)
D (1 7 3 9 10)
перекрутил в голове nth, elt, do. Подкиньте идею с кодом.
(setq po
(capi:contain
(make-instance 'capi:arrow-pinboard-object
:start-x 50 :end-x 50
:start-y 50 :end-y 20
:graphics-args
'(:thickness 1 :foreground :red))))
;;I tried accessors
(pinboard-object-pinboard po)
;;=> #<CAPI:PINBOARD-LAYOUT 402001FCC3>
(pinboard-object-activep po)
;;=> T
(pinboard-object-graphics-args po)
;;=> (:THICKNESS 1 :FOREGROUND :RED)
but I would like to access, say :start-x, change it and draw other arrow in the same window, how to?
;start-x is initarg т.е. я его не могу менять налету, лезу в superclasses но там тоже не могу менять initargs.
#+clim - что это обозначает и далее #-clim-2 #+clim-2
(defvar *clim-clip-rectangle*
#-clim-2 (make-rectangle* 0 0 1 1)
#+clim-2 (make-bounding-rectangle 0 0 1 1)
"Reused by with-clipping-internal to reduce consing.")
make-rectangle* [Function]
Arguments: x1 y1 x2 y2
Summary: Returns an object of class standard-rectangle whose edges are parallel to the coordinate axes. One corner is at the point point1 (or the position (x1, y1)) and the opposite corner is at the point point2 (or the position (x2, y2)). There are no
ordering constraints among point1 and point2 (or x1 and x2, and y1 and y2).
make-bounding-rectangle [Function]
Arguments: x1 y1 x2 y2
Summary: Returns an object of the class standard-bounding-rectangle with the edges specified by x1, y1, x2, and y2, which must be real numbers.
Видать пытаюсь портировать вот это
https://github.com/McCLIM/McCLIM/blob/master/Apps/Scigraph/scigraph/draw.lisp
на LispWorks CAPI
(LAMBDA (L)
(PROG (L1 C)
(SETQ L1 L)
(SETQ C 0)
A (COND ((NULL L1) (RETURN C)))
(SETQ C (ADD1 C))
(SETQ L1 (REST L1))
(GO A)) '(1 2 3 4))
как посмотреть в С?
Переделал
(PROG ((L1 '(1 2 3 4)) (C 0))
A (COND ((NULL L1) (RETURN C)))
(SETQ C (+ 1 C))
(SETQ L1 (REST L1))
(GO A))
4
(funcall (LAMBDA (L)
(PROG (L1 C)
(SETQ L1 L)
(SETQ C 0)
A (COND ((NULL L1) (RETURN C)))
(SETQ C (+ 1 C))
(SETQ L1 (REST L1))
(GO A))) '(1 2 3 4))
4
Вот, наверное, и ответ на вопрос.
Существует ли свободный или не совсем, софт по определению местоположения по номеру мобильника?
Существует проект ОС «Фантом» - https://dz.ru/os-phantom
Ваше мнение.
Кто-нибудь знает какую-нибудь программку которая может отсортировать по списку по регионам ОКАТО, потом по отраслям ОКВЭД.
В Росстат такая точно есть, только друзей у меня там нет.
| ← назад | следующие → |