                IMPLICIT NONE
                REAL::x,vel,F,m,H
                INTEGER::ixx
                INTEGER,PARAMETER::itime = 10000
		REAL,PARAMETER::xposeq = 0.0, val=1.000,delta=0.02

		x = 0.00
		vel  = 1
                F= -val*(x-xposeq)
                m=1.00
                ixx=0

                !write(16,*) ixx*delta,x
                !write(18,*) ixx*delta,vel
                H = 0.5 *(vel**2) + 0.5 * val *(x**2)
                !write(20,*) x,vel
                write(28,*) ixx*delta,H

		do ixx = 1, itime

		vel = vel + (0.5*(F/m)*delta)
                x=x + (vel*delta)
                F = -val*(x-xposeq)
                vel = vel + (0.5*(F/m)*delta)

		!write(16,*) ixx*delta,x
                !write(18,*) ixx*delta,vel
                H = 0.5 *(vel**2) + 0.5 * val *(x**2)
                !write(20,*) x,vel
                write(28,*) ixx*delta,H

		end do

		stop
		end


MODULE crslbclnk
  USE par_kind, ONLY: wp
  USE dom_oce
  USE crs
  USE lbclnk
  USE in_out_manager
  INTERFACE crs_lbc_lnk
    MODULE PROCEDURE crs_lbc_lnk_3d, crs_lbc_lnk_2d
  END INTERFACE
  PUBLIC :: crs_lbc_lnk
  CONTAINS
  SUBROUTINE crs_lbc_lnk_3d(pt3d1, cd_type1, psgn, cd_mpp, pval)
    CHARACTER(LEN = 1), INTENT(IN   ) :: cd_type1
    REAL(KIND = wp), INTENT(IN   ) :: psgn
    REAL(KIND = wp), DIMENSION(jpi_crs, jpj_crs, jpk), INTENT(INOUT) :: pt3d1
    REAL(KIND = wp), OPTIONAL, INTENT(IN   ) :: pval
    CHARACTER(LEN = 3), OPTIONAL, INTENT(IN   ) :: cd_mpp
    LOGICAL :: ll_grid_crs
    REAL(KIND = wp) :: zval
    ll_grid_crs = (jpi == jpi_crs)
    IF (PRESENT(pval)) THEN
      zval = pval
    ELSE
      zval = 0._wp
    END IF
    IF (.NOT. ll_grid_crs) CALL dom_grid_crs
    IF (PRESENT(cd_mpp)) THEN
      CALL lbc_lnk(pt3d1, cd_type1, psgn, cd_mpp, pval = zval)
    ELSE
      CALL lbc_lnk(pt3d1, cd_type1, psgn, pval = zval)
    END IF
    IF (.NOT. ll_grid_crs) CALL dom_grid_glo
  END SUBROUTINE crs_lbc_lnk_3d
  SUBROUTINE crs_lbc_lnk_2d(pt2d, cd_type, psgn, cd_mpp, pval)
    CHARACTER(LEN = 1), INTENT(IN   ) :: cd_type
    REAL(KIND = wp), INTENT(IN   ) :: psgn
    REAL(KIND = wp), DIMENSION(jpi_crs, jpj_crs), INTENT(INOUT) :: pt2d
    REAL(KIND = wp), OPTIONAL, INTENT(IN   ) :: pval
    CHARACTER(LEN = 3), OPTIONAL, INTENT(IN   ) :: cd_mpp
    LOGICAL :: ll_grid_crs
    REAL(KIND = wp) :: zval
    ll_grid_crs = (jpi == jpi_crs)
    IF (PRESENT(pval)) THEN
      zval = pval
    ELSE
      zval = 0._wp
    END IF
    IF (.NOT. ll_grid_crs) CALL dom_grid_crs
    IF (PRESENT(cd_mpp)) THEN
      CALL lbc_lnk(pt2d, cd_type, psgn, cd_mpp, pval = zval)
    ELSE
      CALL lbc_lnk(pt2d, cd_type, psgn, pval = zval)
    END IF
    IF (.NOT. ll_grid_crs) CALL dom_grid_glo
  END SUBROUTINE crs_lbc_lnk_2d
END MODULE crslbclnk


text:
!==========================================================================
elemental function gsw_helmholtz_energy_ice (t, p)
!==========================================================================
!
!  Calculates the Helmholtz energy of ice.
!
!  t  =  in-situ temperature (ITS-90)                             [ deg C ]
!  p  =  sea pressure                                              [ dbar ]
!         ( i.e. absolute pressure - 10.1325 dbar )
!
!  Helmholtz_energy_ice  =  Helmholtz energy of ice                [ J/kg ]
!--------------------------------------------------------------------------

use gsw_mod_toolbox, only : gsw_gibbs_ice

use gsw_mod_teos10_constants, only : gsw_p0, db2pa

use gsw_mod_kinds

implicit none

real (r8), intent(in) :: t, p

real (r8) :: gsw_helmholtz_energy_ice

gsw_helmholtz_energy_ice = gsw_gibbs_ice(0,0,t,p) &
                           - (db2pa*p + gsw_p0)*gsw_gibbs_ice(0,1,t,p)

return
end function

!--------------------------------------------------------------------------

max_stars_repo_name: minsukji/ci-debug

max_issues_repo_name: minsukji/ci-debug

max_forks_repo_name: minsukji/ci-debug


! Copyright (c) 2019, NVIDIA CORPORATION.  All rights reserved.
!
! Licensed under the Apache License, Version 2.0 (the "License");
! you may not use this file except in compliance with the License.
! You may obtain a copy of the License at
!
!     http://www.apache.org/licenses/LICENSE-2.0
!
! Unless required by applicable law or agreed to in writing, software
! distributed under the License is distributed on an "AS IS" BASIS,
! WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
! See the License for the specific language governing permissions and
! limitations under the License.

! Test folding of combined array references and structure component
! references.

module m1
  type :: t1
    integer :: ia1(2)
  end type t1
  type(t1), parameter :: t1x1(*) = [t1::t1(ia1=[1,2]),t1(ia1=[3,4])]
  logical, parameter :: t1check1 = t1x1(2)%ia1(1) == 3
  logical, parameter :: t1check2 = all(t1x1(1)%ia1 == [1,2])
  logical, parameter :: t1check3 = all(t1x1(:)%ia1(1) == [1,3])
  type :: t2
    type(t1) :: dta1(2)
  end type t2
  type(t2), parameter :: t2x1(*) = &
    [t2 :: t2(dta1=[t1::t1x1]), &
           t2(dta1=[t1::t1(ia1=[5,6]),t1(ia1=[7,8])])]
  logical, parameter :: t2check1 = t2x1(1)%dta1(2)%ia1(2) == 4
  logical, parameter :: t2check2 = &
    all(t2x1(2)%dta1(2)%ia1(:) == [7,8])
  logical, parameter :: t2check3 = &
    all(t2x1(1)%dta1(:)%ia1(2) == [2,4])
  logical, parameter :: t2check4 = &
    all(t2x1(:)%dta1(1)%ia1(2) == [2,6])
end module m1
!Expect: m1.mod
!module m1
!type::t1
!integer(4)::ia1(1_8:2_8)
!end type
!type(t1),parameter::t1x1(1_8:*)=[t1::t1(ia1=[Integer(4)::1_4,2_4]),t1(ia1=[Integer(4)::3_4,4_4])]
!logical(4),parameter::t1check1=.true._4
!logical(4),parameter::t1check2=.true._4
!logical(4),parameter::t1check3=.true._4
!type::t2
!type(t1)::dta1(1_8:2_8)
!end type
!type(t2),parameter::t2x1(1_8:*)=[t2::t2(dta1=[t1::t1(ia1=[Integer(4)::1_4,2_4]),t1(ia1=[Integer(4)::3_4,4_4])]),t2(dta1=[t1::t1(ia1=[Integer(4)::5_4,6_4]),t1(ia1=[Integer(4)::7_4,8_4])])]
!logical(4),parameter::t2check1=.true._4
!logical(4),parameter::t2check2=.true._4
!logical(4),parameter::t2check3=.true._4
!logical(4),parameter::t2check4=.true._4
!end


text:
C REINTERPOLATE: Re-interpolate a spectrum onto a sample defined by another
C
       implicit    NONE
C
C The input spectrum is interpolated onto the sampling of a second spectrum, the
C model. The output has exactly the same shape as the model, when data are not
C present in the input then zeros are placed in the output.
C
C Multi column data are handled.  The x-column must be specifed for the
C re-interpolation to proceed, the default is to assume column 1.  Where
C the output data lies in the range of the input data linear interpolation
C between adjacent data values is used.
C
C Parameters:
C    outfile   ==  output file name
C    infile    ==  input file name
C    model     ==  model file against which re-interpolation is done
C    xc        ==  X-column in the infile for use in reinterpolation
C    xcm       ==  X-column in the model file
C    title     ==  new title for output file
C    comment   ==  user defined comment text
C
C Include information on standard array sizes etc.
       include '../include/spec_global.inc'
       include '../include/spec_pars.inc'
C Include standard functions
       include '/mrao/include/chrlib_functions.inc'

C
C Local variables
C   file identifiers
       integer    id_to, id_from, id_model
C   status return
       integer    status
C   number of data points
       integer    nd_from, nd_model
C   array to hold data values
       real*4     x(max_ndata), x_model(max_ndata),
     *            y(max_ndata), y_out(max_ndata)
C   local variables used to hold information and do the interpolation
       real*4     xmin, xmax
C   file and directory name(s) and strings
       character  file*100, dir*100, string*80
C   x-column
       integer    xc, xcm
C   logical to test for setting of input parameters
       logical    test
C   number of columns
       integer    ncols, ncols_model
C   counters
       integer    n, nc, n1, n2

C setup for input
       call io_initio
       call cmd_extpars( 'SPECTRUM', dir, status )
       call cmd_getiline( 'reinterpolate',
     *      dir(1:chr_lenb(dir))//tool_definitions, status)
C find input file names and allocate files
       call cmd_items( 'infile', file, status)
       call spec_allocate( file(1:chr_lenb(file)),
     *                     'read', id_from, status)
       call cmd_items( 'model', file, status)
       call spec_allocate( file(1:chr_lenb(file)),
     *                     'read', id_model, status)
C find output file name and allocate
       call cmd_items( 'outfile', file, status)
       call spec_allocate( file(1:chr_lenb(file)),
     *                     'write', id_to, status)
C copy output descriptors; do this from the model to the output
       call spec_hd_copy( id_model, id_to, status )

C read control data
C .. get the x-column
       call cmd_itemi( 'xc', 1, xc, status )
C .. find model x-column if necessary
       call cmd_itemset( 'xcm', test, status )
       if (test) then
         call cmd_itemi( 'xcm', 1, xcm, status )
       else
         xcm = xc
       endif
C .. read the x-column from the model and input files
       call spec_get_data( id_from, xc, x, nd_from, status )
       call spec_get_data( id_model, xcm, x_model, nd_model, status )
C .. determine and store range of input data; check also that the input
C    and model data are in ascending order -- needed for the interpolation
       xmin = 1.0E+30
       xmax = -1.0E+30
       test = .true.
       n = 1
       xmin = x(1) - 1.0
       do while (test .and. n.le.nd_from)
         test = x(n).gt.xmin
         xmin = x(n)
         n = n + 1
       enddo
       n = 1
       xmin = x_model(1) - 1.0
       do while (test .and. n.le.nd_model)
         test = x_model(n).ge.xmin
         xmin = x_model(n)
         n = n + 1
       enddo
       if (.not.test) then
         print *,'*** Input data and/or model not in ascending order'
         print *,'*** X-columns specified as : ',xc,xcm
         goto 999
       endif
       xmin = x(1)
       xmax = x(nd_from)
C .. add user text to header
       call cmd_items( 'title', string, status )
       if (chr_lenb(string).gt.0) then
         call spec_hd_set( id_to, 'title',
     *                     string(1:chr_lenb(string)), status )
       end if
       call cmd_items( 'comment', string, status )
       if (chr_lenb(string).gt.0) then
         call spec_hd_set( id_to, 'comment',
     *                     string(1:chr_lenb(string)), status )
       end if
C .. find number of columns
       call spec_hd_enqi( id_from, 'ncols', ncols, status )
       call spec_hd_enqi( id_model, 'ncols', ncols_model, status )
C loop for each column in the model
       if (status.ne.0) goto 999
       do nc=1,ncols_model
         if (nc.eq.xc) then
           call spec_get_data( id_model, xcm, y_out, nd_model, status )
         elseif (nc.gt.ncols) then
           call spec_get_data( id_model, xcm, y_out, nd_model, status )
           do n=1,nd_model
             y_out(n) = 0.0
           enddo
         elseif (nc.le.ncols) then
           call spec_get_data( id_from, nc, y, nd_from, status)
           n1 = 1
           do n=1,nd_model
             if (x_model(n).lt.xmin .or. x_model(n).gt.xmax) then
               y_out(n) = 0.0
             else
               do while (x(n1).lt.x_model(n))
                 n1 = n1 + 1
               enddo
               n1 = n1 - 1
               n2 = n1
               do while (x(n2).lt.x_model(n))
                 n2 = n2 + 1
               enddo
               y_out(n) = y(n1) +
     *              (x_model(n)-x(n1))*(y(n2)-y(n1))/(x(n2)-x(n1))
             endif
           enddo
         endif
         if (status.ne.0) goto 999

C .. copy this column to the output
         call spec_put_data( id_to, nc, y_out, nd_model, status)

       end do

C add history item
       string = ' '
       write(string,'(A,2I2)')
     *       'REINTERPOLATE: xc/xcm = ',xc,xcm
       call spec_hd_set( id_to, 'history',
     *                   string(1:chr_lenb(string)),
     *                   status )

C deallocate files
999    continue
       call spec_deallocate( id_from, status )
       call spec_deallocate( id_model, status )
       call spec_deallocate( id_to, status )

       end

ext: f

max_stars_repo_name: CavendishAstrophysics/anmap

max_issues_repo_name: CavendishAstrophysics/anmap

max_forks_repo_name: CavendishAstrophysics/anmap


text:
      PROGRAM xei
C     driver for routine ei
      INTEGER i,nval
      REAL ei,value,x
      CHARACTER text*23
      open(7,file='FNCVAL.DAT',status='OLD')
10    read(7,'(a)') text
      if (text.ne.'Exponential Integral Ei') goto 10
      read(7,*) nval
      write(*,*) text
      write(*,'(1x,t5,a,t13,a,t28,a)')
     *     'X','Actual','EI(X)'
      do 11 i=1,nval
        read(7,*) x,value
        write(*,'(f6.2,2e16.6)') x,value,ei(x)
11    continue
      close(7)
      END

ext: for

max_stars_repo_path: NR-Functions/Numerical Recipes- Example & Functions/Examples/xei.for


text: Lim & Schrader were optometrists based out of Davisville Professional Center. Wayne Schrader OD Wayne C. Schrader continues to practice, but Vic Lim is no longer active.

max_stars_repo_name: voflo/Search

max_issues_repo_name: voflo/Search

max_forks_repo_name: voflo/Search


text:
file fpga_icebreaker.v
file ../libfpga/common/reset_sync.v
file ../libfpga/common/fpga_reset.v
list ../libfpga/common/activity_led.f

list ../soc/soc.f


!----------------------------------------------------------------
!
!> @brief     Reads POE data from files in SP3-c format
!!
!> @author    Vitali Braun
!> @date

!!
 06.02.2014: initial implementation

!!
 14.07.2014: adapted for Sentinel 1b

!!

!!
!> @param[in] reduction_model - used for ECEF -> ECI conversion
!!
!> @param[in] epoch - Epoch for which data has to be provided
!!
!> @param[inout] stateGCRF - Holds the POE states
!!
!> @param[inout] ndata - Defines number of states available
!!
!----------------------------------------------------------------
subroutine readPOE(reduction_model, epoch, stateGCRF, npoe)

    use slam_io
    use slam_math,              only: mag
    use slam_orbit_types,       only: state_t
    use slam_reduction_class,   only: Reduction_type
    use slam_rframes
    use slam_time,              only: gd2mjd, mjd2gd, delta_AT, time_t, mjd2yyddd
    use slam_types

    implicit none

    type(Reduction_type),intent(in)             :: reduction_model
    type(time_t), dimension(2),  intent(in)     :: epoch
    type(state_t), dimension(*), intent(inout)  :: stateGCRF
    integer, intent(out)                        :: npoe

    character(len=255)                          :: cbuf, cbuf2      ! character buffer
    character(len=255), dimension(100)          :: cfile            ! file to read data from

    integer             :: p                    ! counter
    integer             :: ich                  ! channel to read data from
    integer             :: nfiles               ! number of different files to read
    real(dp)            :: dummy, x, y, z, vx, vy, vz
    character(len=35)   :: string_epoch
    logical             :: osv_detected
    integer             :: start_nbr
    integer             :: stop_nbr

    type(state_t) :: state

    write (*,*) "Reading POE"

    osv_detected = .false.

    !=================================================================
    !
    ! Find out data file(s) to read data from based on passed epoch
    !
    !----------------------------------------------------------------

    if(epoch(2)%mjd - epoch(1)%mjd > 10) then
        write(*,*) "Not designed for epochs > 10 days"
        write(*,*) "+++ Program stopped!+++"
        stop
    end if

    ! starting file
    nfiles = 0
    ich = openFile('./poeData/available_files.txt', SEQUENTIAL, IN_FORMATTED)
    do
        call nxtbuf('#', 0, ich, cbuf)
        if (.not. cbuf .eq. '') then
            nfiles = nfiles + 1
            cfile(nfiles) = trim(cbuf)
        else
            exit
        end if
    end do
    close(ich)

    npoe = 0
    do p=1,nfiles

        write (*,*) "Reading file ", trim(cfile(p))
        ich = openFile("./poeData/"//trim(cfile(p)), SEQUENTIAL, IN_FORMATTED)
        do
            call nxtbuf('#', 0, ich, cbuf)
            if (.not. cbuf .eq. '') then

                ! Scan for XML-Tag
                if (index(cbuf,'<OSV>')/=0) then
                    osv_detected = .true.
                elseif (index(cbuf,'</OSV>')/=0) then
                    osv_detected = .false.
                    if (npoe < 2) then
                        ! Make sure that duplicate epoch entries are ignored
                        npoe = npoe + 1
                        call reduction_model%earthFixed2inertial(state%r, state%v, state%epoch%mjd ,stateGCRF(npoe)%r, stateGCRF(npoe)%v)
                        stateGCRF(npoe)%epoch = state%epoch
                    else if (state%epoch%mjd > stateGCRF(npoe)%epoch%mjd) then
                        ! Make sure that duplicate epoch entries are ignored
                        npoe = npoe + 1
                        call reduction_model%earthFixed2inertial(state%r, state%v, state%epoch%mjd ,stateGCRF(npoe)%r, stateGCRF(npoe)%v)
                        stateGCRF(npoe)%epoch = state%epoch
                    end if
                endif

                if (osv_detected .eqv. .true.) then
                    if(index(cbuf,'<X unit="m">')/=0 ) then
                        start_nbr = index(cbuf,'<X unit="m">',.true.)+12
                        stop_nbr = index(cbuf,'</X>',.true.) -1
                        read(cbuf(start_nbr:stop_nbr),"(1f20.6)") dummy
                        read(cbuf(start_nbr:stop_nbr),*) x
                        state%r(1) = x/1000.0

                    elseif(index(cbuf,'<Y unit="m">')/=0 ) then
                        start_nbr = index(cbuf,'<Y unit="m">',.true.)+12
                        stop_nbr = index(cbuf,'</Y>',.true.) -1
                        read(cbuf(start_nbr:stop_nbr),"(1f20.6)") dummy
                        read(cbuf(start_nbr:stop_nbr),*) y
                        state%r(2) = y/1000.0

                    elseif(index(cbuf,'<Z unit="m">')/=0 ) then

                        start_nbr = index(cbuf,'<Z unit="m">',.true.)+12
                        stop_nbr = index(cbuf,'</Z>',.true.) -1
                        read(cbuf(start_nbr:stop_nbr),"(1f20.6)") dummy
                        read(cbuf(start_nbr:stop_nbr),*) z
                        state%r(3) = z/1000.0

                    elseif(index(cbuf,'<VX unit="m/s">')/=0 ) then

                        start_nbr = index(cbuf,'<VX unit="m/s">',.true.)+15
                        stop_nbr = index(cbuf,'</VX>',.true.) -1
                        read(cbuf(start_nbr:stop_nbr),"(1f20.6)") dummy
                        read(cbuf(start_nbr:stop_nbr),*) vx
                        state%v(1) = vx/1000.0

                    elseif(index(cbuf,'<VY unit="m/s">')/=0 ) then

                        start_nbr = index(cbuf,'<VY unit="m/s">',.true.)+15
                        stop_nbr = index(cbuf,'</VY>',.true.) -1
                        read(cbuf(start_nbr:stop_nbr),"(1f20.6)") dummy
                        read(cbuf(start_nbr:stop_nbr),*) vy
                        state%v(2) = vy/1000.0

                    elseif(index(cbuf,'<VZ unit="m/s">')/=0 ) then

                        start_nbr = index(cbuf,'<VZ unit="m/s">',.true.)+15
                        stop_nbr = index(cbuf,'</VZ>',.true.)-1
                        read(cbuf(start_nbr:stop_nbr),"(1f20.6)") dummy
                        read(cbuf(start_nbr:stop_nbr),*) vz
                        state%v(3) = vz/1000

                    elseif(index(cbuf,'<UTC>')/=0 ) then

                        start_nbr = index(cbuf,'<UTC>',.true.)+9
                        stop_nbr = index(cbuf,'</UTC>',.true.)-1
                        read(cbuf(start_nbr:stop_nbr),*) string_epoch
                        read(cbuf(start_nbr:stop_nbr),*) string_epoch
                        read(string_epoch(1:4),*)   state%epoch%year
                        read(string_epoch(6:7),*)   state%epoch%month
                        read(string_epoch(9:10),*)  state%epoch%day
                        read(string_epoch(12:13),*) state%epoch%hour
                        read(string_epoch(15:16),*) state%epoch%minute
                        read(string_epoch(18:19),*) state%epoch%second
                        call gd2mjd(state%epoch)

                    elseif(index(cbuf,'<Quality>')/=0 ) then

                        start_nbr = index(cbuf,'<Quality>',.true.)+9
                        stop_nbr = index(cbuf,'</Quality>',.true.)-1
                        read(cbuf(start_nbr:stop_nbr),*) cbuf2
                    endif
                endif

                if (index(cbuf,'</List_of_OSVs>')/=0) then
                    exit
                endif
            else
                exit
            end if
        end do

        close(ich)

    end do

    return

end subroutine readPOE


text:
Debats du Senat (hansard)
1ere Session, 36e Legislature,
Volume 137, Numero 17
Le jeudi 6 novembre 1997
L'honorable Gildas L. Molgat, President
L'accord sur les revendications territoriales des Inuit
Presentation du premier rapport du comite mixte-Ajournement du debat
La loi sur le Parlement du Canada
Projet de loi modificatif-Rapport du comite
Privileges, Reglement et procedure
Huissier du Senat-Presentation du rapport du comite
Avis de motion portant autorisation au comite d'engager du personnel et de se deplacer
La protection de la conscience dans les procedures medicales-Presentation d'une petition
L'etude des effets sur la sante de l'additif MTBE-La position du gouvernement
Reponses differees a des questions orales
Depot des reponses a des questions au Feuilleton
L'energie-Peches et Oceans-La conformite avec la Loi sur les carburants de remplacement
Adoption du premier rapport du comite mixte permanent
Privileges, Reglement et procedure
L'huissier du Senat-Presentation et adoption du rapport modifie du comite
Recours au Reglement-Decision de la presidence
Sous-amendement a la motion d'amendement
Le jeudi 6 novembre 1997
La seance est ouverte a 14 heures, le President etant au fauteuil.
L'honorable B. Alasdair Graham (leader du gouvernement) :
Nos murs font echo a leurs sacrifices.
Le pays a pleure sa mort.
Des temoignages sont venus de partout dans le Dominion.
En cette enceinte particuliere, nous nous souvenons d'eux.
Nous entendons leurs voix et leurs paroles.
L'honorable Orville H. Phillips :
En tant que parlementaires, nous ne devons jamais oublier l'histoire de nos guerres.
Ils en ont garde les cicatrices toute leur vie.
Il faut egalement prendre garde a ne pas minimiser l'importance de ces guerres mondiales.
Que nous oublions vite.
De nombreuses universites ne celebrent pas le jour du Souvenir.
Le 11 novembre est un jour sacre et doit toujours le demeurer.
Ils ne vieilliront pas comme nous qui leur avons survecu;
Ils ne connaitront plus l'outrage ni le poids des annees.
Quand viendra l'heure du crepuscule et celle de l'aurore,
Nous nous souviendrons d'eux.
L'honorable Richard J. Doyle :
Que la premiere annee est la pire.
Puis, que l'on s'y fait.
Que l'on s'y fait?
La derniere se termine par ces phrases:
L'armistice a ete signe a minuit par les representants allemands.
On nous dit que la grosse sirene d'Eaton a reveille la ville tout entiere.
Dominant la fureur de la guerre
Font retentir les cieux de leur chant:
Paix sur terre aux combattants.
Des blagues et un autre poeme sur l'achat d'obligations constituent l'editorial.
L'article principal, intitule La fin de l'autocratie, contient cette prophetie:
L'empire allemand n'est plus, mais le peuple allemand demeure.
C'est avec lui que les allies doivent conclure un reglement permanent.
Et voici une citation de Lloyd George:
L'insouciance qui a plonge le monde dans cette agonie atroce doit etre severement jugee.
La manchette disait This is Victory, en caractere Jesus, comme disent les imprimeurs.
C'etait la manchette.
Les cymbales retentissent et les danseuses s'avancent,
Avec de longs bas de soie et des bras d'albatre,
Et les ombres d'hommes disparus les regardent.
Une ombre dit: Que pensiez-vous trouver,
Quand le dernier coup fut tire et que la paix fut faite?
Son ami repond, de son sourire absent: Diable,
Je pensais qu'ils prieraient pour la reconciliation.
L'honorable Nicholas W. Taylor :
J'ai servi dans la marine, ou j'avais le matricule B-3372.
Pourtant, ils furent parmi nos meilleurs soldats.
Ils furent parmi les plus decores.
Nous avons de profondes excuses a faire.
La motion a ete adoptee et le gouvernement a maintenant cree un comite des bourses.
J'ai assume les fonctions de president du comite du Monument des anciens combattants autochtones.
Nous avons verse l'argent necessaire a sa conception.
Il sera haut de trois etages environ.
Un compte a ete ouvert a la banque.
Je vais donner l'exemple.
Nous erigerons un monument.
Ils sont partis la-bas et ils se sont fort bien acquittes de leur mission.
Son Honneur le President :
Les honorables senateurs observent un moment de silence.
Son Honneur le President :
Nous accueillons egalement Son Excellence Bienvenido Garcia Negrin, ambassadeur de Cuba.
Au nom de tous les honorables senateurs, je vous souhaite la bienvenue au Senat.
L'honorable Eymard G. Corbin :
L'article dit egalement:
Ce serait une arithmetique assez simple.
J'applaudirais d'un cote et je pleurerais peut-etre de l'autre.
Il y a la blessure qui ne s'est pas tout a fait refermee.
Ce n'est pas mon cas.
J'ai un vif souvenir de l'epoque et des evenements.
Andre Laurendeau etait decede.
Louis J. Robichaud venait d'etre reconduit au pouvoir pour un troisieme mandat.
On n'oublie pas de sitot!
J'avais suivi a la tele de Radio-Canada le deroulement des etats generaux.
L'affaire etait paquetee d'avance.
Pourtant, nous sommes toujours debout et nous consolidons nos forces.
J'ai eu mes hauts et mes bas.
Il y a du fiel dans mon propos, c'est evident.
Mais je marche la tete haute neanmoins.
Je suis un citoyen a part entiere avec des droits.
Je m'attends au respect reciproque.
Que voulez-vous, j'ai vecu cette periode.
Je me souviens, mais j'ai toujours espoir.
L'honorable Philippe Deane Gigantes :
Il a repondu: Mais les Quebecois parlent le francais.
Ils veulent se rattacher a la France et quitter la Belgique.
On a demande au president Jacques Chirac de France de commenter ce phenomene.
Ma reponse a M. Domljan permettra d'eclairer en partie une question posee precedemment.
Nous preferons la notion de souverainete.
Nous ne sommes pas un Etat federal.
La question des langues regionales se pose en termes differents ici.
L'accord sur les revendications territoriales des Inuit
Cette entente sera la suite de cette tradition.
Il y avait beaucoup d'emotions hier, pour deux raisons.
Son Honneur le President :
Je remercie encore une fois tous les senateurs de leur participation.
Presentation du premier rapport du
L'honorable Celine Hervieux-Payette :
(Le texte du rapport figure dans les Journaux du Senat d'aujourd'hui.)
Son Honneur le President :
Honorables senateurs, quand etudierons-nous ce rapport?
Son Honneur le President :
Vous plait-il, honorables senateurs, d'accorder la permission?
Son Honneur le President :
Il est propose par l'honorable senateur Hervieux-Payette que le rapport soit etudie maintenant.
Vous plait-il, honorables senateurs, d'adopter la motion?
Son Honneur le President :
Le Senat demande une explication si vous le voulez bien, senateur Hervieux-Payette.
L'honorable Noel A. Kinsella (chef adjoint suppleant de l'opposition) :
A-t-on demande qu'on s'eloigne de cette facon de faire?
Son Honneur le President :
Le senateur Hervieux-Payette en a demande la permission.
Je l'ai dit en francais.
Tenons-nous actuellement le debat?
Son Honneur le President :
Je propose l'ajournement du debat.
Son Honneur le President :
Plait-il aux honorables senateurs d'adopter la motion?
(Sur la motion du senateur Kinsella, le debat est ajourne.)
La loi sur le Parlement du Canada
Projet de loi modificatif-Rapport du comite
Le jeudi 6 novembre 1997
Le comite senatorial permanent des affaires juridiques et constitutionnelles a l'honneur de presenter son
Son Honneur le President :
Honorables senateurs, quand lirons-nous le projet de loi une deuxieme fois?
Privileges, Reglement et procedure
Huissier du Senat-Presentation du rapport du comite
Le jeudi 6 novembre 1997
Son Honneur le President :
Honorables senateurs, quand etudierons-nous le rapport?
Son Honneur le President :
Son Honneur le President :
Vous plait-il, honorables senateurs, d'adopter la motion?
Son Honneur le President :
Son Honneur le President :
Je dois donc demander qui est en faveur et qui est contre.
Que tous ceux qui sont en faveur veuillent bien dire oui.
Son Honneur le President :
Que tous ceux qui sont contre veuillent bien dire non.
Son Honneur le President :
L'etude du rapport sera inscrite au Feuilleton d'aujourd'hui.
L'honorable Sharon Carstairs (leader adjoint du gouvernement) :
Son Honneur le President :
La permission est-elle accordee, honorables senateurs?
L'honorable M. Lorne Bonnell :
Honorables senateurs, je donne avis que, le mardi 18 novembre 1997, je proposerai:
L'honorable John B. Stewart :
Honorables senateurs, je donne avis que, le mercredi 19 novembre 1997, je proposerai:
L'honorable Leonard J. Gustafson :
Honorables senateurs, je donne avis que, le mardi 18 novembre 1997, je proposerai:
Que le comite presente son rapport au plus tardle 15 decembre 1998.
L'honorable Leonard J. Gustafson :
Honorables senateurs, je donne avis que, le mardi 18 novembre 1997, je proposerai:
Que le comite presente son rapport au plus tardle 15 decembre 1998.
Avis de motion portant autorisation au comite d'engager du personnel et de se deplacer
L'honorable Leonard J. Gustafson :
Honorables senateurs, je donne avis que, le mardi 18 novembre 1997, je proposerai:
L'honorable Leonard J. Gustafson :
Honorables senateurs, je donne avis que, le mardi 18 novembre 1997, je proposerai:
La protection de la conscience dans les procedures medicales-Presentation d'une petition
L'honorable John Lynch-Staunton (chef de l'opposition) :
Est-elle encore en cours et, dans l'affirmative, ou en est-elle?
Par exemple, des mandats de perquisition ont-ils ete demandes?
Est-on sur le point de porter des accusations?
Quel est l'etat actuel de l'enquete Airbus?
L'honorable B. Alasdair Graham (leader du gouvernement) :
Ils se sont arroge le merite de l'enquete.
Ils ont parle des mandats de perquisition.
Ce sont toutes deux des enquetes de la GRC.
Dans ce cas-ci, c'est le silence complet.
L'honorable J. Michael Forrestall:
Essayez-vous de cacher quelque chose dont nous n'avons pas entendu parler?
L'honorable B. Alasdair Graham (leader du gouvernement):
Le voila maintenant en avance sur le gouvernement.
Le senateur Forrestall siege peut-etre du mauvais cote.
Je ne sais cependant pas exactement a quel moment ou quelle semaine cela se fera.
Nous ne pouvons que nous interroger ensemble a ce sujet.
Aucune decision n'a encore ete prise touchant l'helicoptere de remplacement.
On en fera un examen complet a ce moment-la.
L'etude des effets sur la sante de l'additif MTBE-La position du gouvernement
L'honorable Erminie J. Cohen :
L'honorable B. Alasdair Graham (leader du gouvernement) :
Honorables senateurs, je ne sais pas si une garantie sera donnee.
C'est une question importante, et je tenterai d'obtenir tres bientot une reponse.
L'honorable Erminie J. Cohen :
L'honorable B. Alasdair Graham (leader du gouvernement) :
Je vais prendre note de cette question, honorables senateurs.
L'honorable Donald H. Oliver :
Combien de Canadiens devront payer le prix de cette augmentation prevue des cotisations?
L'honorable B. Alasdair Graham (leader du gouvernement) :
A ma connaissance, il n'est pas question de perte d'emplois, honorables senateurs.
Le leader du gouvernement a-t-il dit qu'il deposerait des etudes?
J'ai dit que, a ma connaissance, il n'etait pas...
Une reponse de ce genre merite d'etre murmuree.
Le senateur Lynch-Staunton voudrait peut-etre parler.
A ma connaissance, il n'est pas question de perte d'emplois.
Elle portait sur la vente proposee de reacteurs CANDU a la Turquie.
L'honorable B. Alasdair Graham (leader du gouvernement) :
Honorables senateurs, la technologie du reacteur CANDU est sure.
Ce reacteur est le meilleur et le plus efficace au monde.
Nous devrions en etre fier.
La question s'adressait au ministre.
N'etant jamais alle en Roumanie, il m'est difficile de comparer.
L'honorable John Lynch-Staunton :
Dans ce cas, pourquoi procede-t-on a une evaluation environnementale?
L'honorable A. Raynell Andreychuk :
Autrement dit, nous n'avons absolument aucune garantie que des evaluations environnementales seront faites.
L'honorable B. Alasdair Graham (leader du gouvernement) :
L'honorable John Lynch-Staunton :
Pouvons-nous voir les resultats?
S'ils sont disponibles, oui.
L'honorable B. Alasdair Graham (leader du gouvernement) :
Honorables senateurs, la reponse a cette question est oui.
Je voudrais souligner que nous voulons eviter une solution venant du gouvernement federal seulement.
Ce sera une solution basee sur un consensus au Canada.
Il faut esperer que toutes les provinces et tous les interesses participeront au processus.
En toute deference, je dois dire que cette reponse n'a aucun sens.
On met la charrue devant les boeufs.
Ce serait presomptueux de ma part.
Je considere la question tres importante pour le Canada.
Ma reponse a l'honorable senateur est affirmative.
L'honorable B. Alasdair Graham (leader du gouvernement):
Honorables senateurs, des discussions sont en cours a cet egard.
J'ignore si je pourrai repondre avant que la decision soit prise.
Des que j'aurai des renseignements, je les communiquerai au Senat.
L'honorable B. Alasdair Graham (leader du gouvernement) :
Honorables senateurs, je devrai prendre note de la question.
Autrement dit, nous avons eu deux inondations.
L'hiver est encore precoce, cette annee, au Manitoba, pour la troisieme annee d'affilee.
Son Honneur le President :
On ne peut pas adresser les questions a d'autres senateurs.
Je pourrais peut-etre la poser au leader du gouvernement au Senat.
Le Dakota du Nord a eu des precipitations au moins aussi abondantes.
Le gouvernement surveille-t-il le moindrement la situation?
L'honorable B. Alasdair Graham (leader du gouvernement) :
Honorables senateurs, la reponse est oui.
Je remercie le senateur pour ses observations tres legitimes.
Le gouvernement suit-il les travaux de la commission mixte?
Beaucoup de gens attendent la publication d'un rapport.
Honorables senateurs, le gouvernement suit les travaux de la commission mixte.
Je ne sais pas quand le rapport sera pret.
C'est avec plaisir que je ferai cela pour le senateur.
J'essaierai d'obtenir une reponse pour notre retour de conge.
L'honorable Eymard G. Corbin :
Je crois que personne ne voulait cela.
Certes, nous lui devons cela apres toutes ces annees.
A un moment donne, les peintures vont commencer a peler et a s'effriter.
L'honorable B. Alasdair Graham (leader du gouvernement) :
Son Honneur le President :
Toutefois, pour votre gouverne, cette question a ete portee a mon attention l'an dernier.
Il a declare qu'elles etaient en bon etat.
Vous remarquerez des lignes noires.
On nous a dit qu'elles etaient causees par l'eclairage.
L'honorable J. Michael Forrestall :
Honorables senateurs, ce n'est pas tout a fait ce que j'ai compris.
Les peintures doivent etre restaurees.
Il serait tres difficile de les enlever.
Il se peut que l'on soit oblige de les restaurer sur place.
Le President pourrait peut-etre nous eclairer davantage sur ce sujet.
Son Honneur le President :
Honorables senateurs, nous vous ferons un rapport complet.
L'honorable Donald H. Oliver :
L'honorable B. Alasdair Graham (leader du gouvernement) :
Reponses differees a des questions orales
L'honorable Sharon Carstairs (leader adjoint du gouvernement) :
( Reponse a la question posee par l'honorable Mira Spivak le 21 octobre 1997 )
Le perfectionnement professionnel, la reconnaissance et la formation sont au nombre des questions visees.
On n'envisage pas pour l'instant que la CGC obtienne un statut special.
L'energie-Revenu national-La conformite avec la Loi sur les carburants de remplacement
L'energie-Peches et Oceans-La conformite avec la Loi sur les carburants de remplacement
Adoption du premier rapport du comite mixte permanent
L'honorable Sharon Carstairs (leader adjoint du gouvernement) :
Son Honneur le President :
Plait-il aux honorables senateurs d'adopter la motion?
(La motion est adoptee, et le rapport est adopte.)
Privileges, Reglement et procedure
Le jeudi 6 novembre 1997
Son Honneur le President :
Plait-il aux honorables senateurs d'adopter la motion?
L'honorable Eymard G. Corbin :
Je crois qu'il est de notre devoir de nous poser la question.
Son Honneur le President :
On n'emploie pas toujours le reste de l'appellation de la Verge noire.
C'est tellement stupide.
On parle du huissier a maintes reprises sans utiliser l'expression de la Verge noire.
On dit simplement huissier.
Personnellement, je ne suis pas prete a faire marche arriere.
Nous devrions dire huissier du Senat.
Le Senat est autonome.
C'est pourquoi le rapport est redige ainsi.
Le senateur Corbin desire-t-il plus d'information sur la question?
Honorables senateurs, je remercie ma collegue de ses explications.
Je pense qu'il y a de l'hypocrisie dans cette designation.
Une verge est une verge.
Si nous voulons remonter dans l'histoire et parler des dates...
J'ai l'impression de manquer une plaisanterie la-bas.
J'entends beaucoup de ricanements, mais il n'y a pas matiere a rire.
Je vous vois tous sourire.
J'aurais aime que le senateur Gigantes soit ici.
A l'epoque de l'empire romain...
Le premier sens que l'on donne a verge c'est:
On ne l'a jamais questionnee.
Il y en a qui ont des idees de travers.
Ils voient un mot et puis ils pensent a autre chose.
Ils voient le mot verge et puis ils pensent a autre chose.
Ils ne pensent pas a la baguette.
C'est bien dommage.
Nous aurons sans doute l'occasion de lui rendre l'hommage de circonstance qui convient.
C'est faire des detours, des soubresauts, des jambettes inutiles, je crois.
Je n'irai pas jusqu'a vouloir bloquer cette nouvelle designation.
L'hypocrisie est la peur du sens veritable des mots.
Je trouve cela absolument regrettable.
Ce sont de petits esprits qui ont pense a cela et je le regrette.
L'honorable Anne C. Cools :
Son Honneur le President :
Senateur Cools, vous ne pouvez poser une question qu'au dernier orateur.
Je peux peut-etre poser ma question au dernier orateur.
Tout d'abord, ce sujet me prend totalement par surprise.
Je souhaite participer au debat, mais je demande un eclaircissement.
Nous avons entendu mentionner le discours du senateur Lang.
Habituellement, les rapports de comite disent: Le comite, conformement a un renvoi...
S'il y avait un renvoi du Senat, cela m'etonne.
Cela m'a prise par surprise au cours de la derniere heure ou deux.
Y avait-il un renvoi?
Le rapport a-t-il ete renvoye precedemment au comite?
Je pourrais peut-etre soulever cette question dans un rappel au Reglement.
Son Honneur le President :
Je suis desole, honorable senateur Cools.
Je parcourais mon Reglement.
Honorables senateurs, je devrais peut-etre soulever cette question dans un rappel au Reglement.
J'essayais d'eviter cela.
Son Honneur le President :
Le senateur Maheu, presidente du comite, peut repondre a cette question.
L'honorable Noel A. Kinsella (chef adjoint suppleant de l'opposition) :
Honorables senateurs, une reflexion est en train de se faire.
Le decret dont il est question dans le rapport est date du17 octobre.
Le senateur Corbin a mis le doigt sur le point nevralgique.
En trois jours, on a tente de modifier une tradition de longue date.
Il me semble que cette question est du ressort des traducteurs.
Puis-je poser une question au senateur Kinsella?
Le senateur Kinsella avait acces a un document que je n'ai pas vu.
Peut-il me dire si Sa Majeste a ete consultee a cet egard?
Le Gentilhomme huissier de la Verge noire est un de ses serviteurs au Canada.
Je remercie le senateur de sa question.
Je crains toutefois de ne pouvoir lui fournir une reponse satisfaisante.
Je n'ai pas de connaissances dans ce domaine.
Son Honneur le President :
La permission est-elle accordee de deposer ces decrets?
Honorables senateurs, il me semble que nous soyons saisis de deux questions cet apres-midi.
Je vais d'abord parler de la deuxieme question.
C'est le probleme que le comite du Reglement a tente de regler ce matin.
(iii) etudier les ordres et coutumes du Senat...
Je voudrais poser une question au leader adjoint du gouvernement.
Donc, l'objet reste, mais le titre est appele a disparaitre.
Qu'avons-nous gagne?
Qu'avons-nous perdu?
Je le dis avec politesse.
Senateur Corbin, je conviens que la demarche a, a mon avis, quelque chose d'incongru.
J'ai une question pour l'honorable leader adjoint.
La proposition de l'honorable senateur est interessante.
Evidemment, on ne nous a pas repondu.
Ce sont eux qui ont etabli la commission et publie les deux decrets.
Il y a surement une explication pour leurs gestes.
Comment pouvons-nous les connaitre?
L'honorable Nicholas W. Taylor :
Honorables senateurs, brievement, la chose me derange moi aussi.
La modernisation du langage peut parfois contraindre a des acrobaties.
Je deviens tres sceptique lorsqu'il est question de changer les traditions.
Cela devrait etre au Senat de le faire.
Pourquoi ne pas envisager par exemple le titre d'orateur?
Son Honneur le President :
Plait-il aux honorables senateurs d'adopter la motion?
Que les honorables senateurs en faveur de la motion veuillent bien dire oui.
Son Honneur le President :
Que les honorables senateurs qui sont contre la motion veuillent bien dire non.
Son Honneur le President :
L'honorable John Lynch-Staunton (chef de l'opposition) :
Honorables senateurs, j'invoque le Reglement.
Le rapport ne nous est parvenu qu'il y a peu de temps.
Je juge cela tout a fait inacceptable.
Si cela est autorise, c'est une forme de cloture.
Lorsque le debat prend fin, la motion est mise aux voix.
J'insiste, cela est hautement irregulier.
Cela va a l'encontre des coutumes de cet endroit.
Son Honneur le President :
L'honorable Anne C. Cools :
Honorables senateurs, j'appuie le rappel au Reglement du senateur Lynch-Staunton.
J'aimerais dire une ou deux choses.
Toutefois, j'aimerais aujourd'hui porter quelques points a l'attention des senateurs.
Un certain nombre de questions ont ete soulevees.
Le senateur Carstairs essaye, avec beaucoup de zele je trouve, de faire un excellent travail.
Honorables senateurs, nous devrions reflechir comme il faut a ces questions.
Je crois que le rappel au Reglement du senateur Lynch-Staunton est fonde.
Son Honneur le President :
D'autres senateurs veulent-ils prendre la parole sur ce rappel au Reglement?
Honorables senateurs, je suppose que M. Kinsella ou moi-meme pouvons eclaircir la question.
Son Honneur le President :
C'est l'objet du rappel au Reglement dont nous sommes saisis.
Il s'agit uniquement de refuser ou d'accepter l'ajournement du debat.
L'honorable Noel A. Kinsella (chef adjoint suppleant de l'opposition) :
Le Reglement definit le fonctionnement des comites.
C'est la que le probleme a commence.
C'est ce qui s'est produit dans le cas present.
Le rapport du comite a ete presente.
Voila ou nous en sommes.
Cette motion d'ajournement a ete entouree d'une serie d'irregularites.
Or, ce n'est pas le leader adjoint qui a refuse l'ajournement du debat.
Je croyais necessaire d'apporter cette precision.
Recours au Reglement-Decision de la presidence
Son Honneur le President :
Le senateur Lynch-Staunton disait que je sympathiserais avec lui.
Il est toujours possible.
Par consequent, je m'en remets au Senat.
Le Reglement ne dit rien a ce sujet.
C'est peut-etre une question que le comite du Reglement pourrait considerer.
C'est peut-etre une question qui a ete oubliee.
Pour le moment, toutefois, il n'y a rien dans le Reglement.
Par suite, on peut accepter ou refuser une demande d'ajournement.
Je ne peux rien dire de plus.
Je voudrais parler du point souleve par le senateur Cools.
Franchement, je n'ai rien entendu.
A mon avis, le rappel au Reglement est irrecevable.
Le debat peut se poursuivre.
L'honorable John B. Stewart :
Honorables senateurs, je me demande si nous examinons la bonne question.
Il n'appartient pas au Senat de repondre a cette question.
De tout temps, les divers hauts fonctionnaires du palais avaient differents symboles de fonctions.
Il est triste de voir notre passe abandonne.
Nous sommes peut-etre horrifies de la designation huissier du Senat.
Que signifie-t-elle?
Cette designation est certes tres inexacte.
Le haut fonctionnaire exerce bien d'autres fonctions hormis celle de huissier du Senat.
Neanmoins, je crois qu'il a ete choisi par la filiere appropriee.
L'honorable Nicholas W. Taylor :
Aurais-je votre permission pour ce faire?
Son Honneur le President :
Si les honorables senateurs consentent a recevoir une question, c'est d'accord.
Je pense que c'est tout a fait le cas qui nous occupe.
Je suis membre du comite permanent des privileges, du Reglement et de la procedure.
Qu'a fait le comite?
La solution a ete presentee cet apres-midi par le senateur Maheu.
Aujourd'hui, la meme expression a un sens tout a fait different.
Elle a acquis un sens completement nouveau au fil des annees.
Honorables senateurs, que s'est-il passe aujourd'hui?
Que nous a dit madame le leader adjoint du gouvernement aujourd'hui?
Cette approche me semble raisonnable.
J'exhorte le senateur Lynch-Staunton a essayer de comprendre la situation.
Ce qui s'est passe aujourd'hui n'a rien de malicieux.
Nous essayons tout simplement de nous occuper d'une situation qui vient de se manifester.
J'exhorte l'honorable chef de l'opposition a accepter la proposition.
Je ne savais pas quel etait l'objet de cette reunion speciale.
C'est ce que nous avons fait de bonne foi.
L'honorable John Lynch-Staunton (chef de l'opposition) :
Honorables senateurs, je trouve cet echange non seulement incroyable, mais surrealiste.
Je ne comprends pas cela.
Il y a ici deux dictionnaires, senateur, si vous voulez verifier.
Pourquoi ne pas aller plus loin?
Jusqu'ou allons-nous pousser cet echange ridicule?
Honorables senateurs, je maintiens que nous devrions etre tout a fait satisfaits du titre actuel.
Il est juste et il respecte une tradition de longue date.
L'honorable John Lynch-Staunton :
En consequence, dans cette optique, je propose un amendement, appuye par le senateur Kinsella:
sous ce titre en qualite de haut fonctionnaire
et en les remplacant par:
sous le titre de huissier de la Verge noire.
En anglais, le rapport du comite serait alors libelle en ces termes:
En francais, le rapport serait ainsi libelle:
Cela respecterait la tradition et denouerait l'impasse actuelle.
Son Honneur le President :
Plait-il aux honorables senateurs d'adopter la proposition d'amendement?
L'honorable John B. Stewart :
Nous avons un decret du gouverneur en conseil.
Son Honneur le President :
On me demande si la motion est recevable.
Je l'ai lue aussi attentivement que j'ai pu.
Est-ce d'accord?
Son Honneur le President :
La seance reprendra a 17 heures.
(Le Senat s'ajourne a loisir.)
Reprise de la seance a 17 h 06.
Son Honneur le President:
Je le declare donc recevable.
Je vais entendre la suite du debat.
L'honorable Philippe Deane Gigantes :
Honorables senateurs, je voudrais revenir sur l'amendement du senateur Lynch-Staunton.
Sous-amendement a la motion d'amendement
L'honorable Philippe Deane Gigantes :
que l'on modifie la motion d'amendement, en francais, en supprimant les mots:
huissier de la Verge noire
et en les remplacant par:
huissier du Baton noir
L'honorable Anne Cools :
Son Honneur le President :
Son Honneur le President :
Vous plait-il, honorables senateurs, d'adopter les sous-amendements?
Voulez-vous que je cite a nouveau l'amendement?
Son Honneur le President :
Plait-il aux honorables senateurs d'adopter l'amendement?
Son Honneur le President :
Je vais la citer.
Honorables senateurs, est-ce que cela veut dire que nous gardons le mot gentilhomme?
Je vais lire la motion pour m'assurer qu'elle est claire.
La motion executoire telle que modifiee et remodifiee se lit comme ceci:
Et en francais, cela se lirait:
Son Honneur le President :
Plait-il aux senateurs d'adopter la motion telle que modifiee?
Honorables senateurs, je tiens a preciser que je m'abstiens de voter.
Honorables senateurs, je tiens a le preciser aussi.
Son Honneur le President :
Je ne m'oppose pas.
Ce n'est pas du tout la meme chose.
Je le precise officiellement, un point c'est tout.
Son Honneur le President :
Je devrais sans doute preciser que je ne vous entendrai pas.
(La motion est adoptee avec dissidence, et le rapport modifie est adopte.)
Le Senat s'ajourne a 14 heures, le mardi 18 novembre 1997.


text:
!-----collection of functions--------------------------------g.guidati--
!
!
!-----------------------------------------------------------------------
!
!     franco guidati    IAG
!                       Uni Stuttgart
!
!
!     scope             functions for odeint
!
!.......................................................................
!               declarations
!.......................................................................

      subroutine SIG(s,y,dydx)

      USE TINoiseGeneric
      USE TINoiseGeo
      USE TINoiseDDD
      USE TIPrecision

IMPLICIT                        NONE

   ! Local variables.

INTEGER(4)                   :: k,khi,klo,i

REAL(DbKi)                   :: s,y1,y2,d1s,d2s,r2,r4,r6
REAL(DbKi)                   :: n1, n2,d1,d2,d1y1,d1y2,adds,dump
REAL(DbKi)                   :: y(6), dydx (6)

      call SPL_EX  (swork,n+1,s,khi,klo)
      call SPL_EX1 (swork,yc1,d2yc1,n+1,s,y1,d1y1,khi,klo)
      call SPL_EX1 (swork,yc2,d2yc2,n+1,s,y2,d1y2,khi,klo)
      CALL SPL_EX1 (swork, pots, d2pots, n+1, s, adds, dump, khi, klo)

      n1 =  d1y2
      n2 = -d1y1
      d1 = x1-y1
      d2 = x2-y2
      d1s = d1**2
      d2s = d2**2
      r2 = d1s + d2s
      r4 = r2 ** 2
      r6 = r2 ** 3

      dydx(1) =  (n1*d1+n2*d2)/r2* adds * pi2i
      dydx(2) =  (n1*(d2s-d1s)-2.0d0*n2*d2*d1)/r4 * adds * pi2i
      dydx(3) =  (n2*(d1s-d2s)-2.0d0*n1*d2*d1)/r4 * adds * pi2i
      dydx(4) = +(d1s*(2.0d0*n1*d1+6.0d0*n2*d2)- d2s*(6.0d0*n1*d1+2.0d0*n2*d2))/r6* adds *pi2i
      dydx(5) = +(d1s*(6.0d0*n1*d2-2.0d0*n2*d1)+ d2s*(6.0d0*n2*d1-2.0d0*n1*d2))/r6* adds *pi2i
      dydx(6) = +(d2s*(2.0d0*n2*d2+6.0d0*n1*d1)- d1s*(6.0d0*n2*d2+2.0d0*n1*d1))/r6* adds *pi2i

      return
      end

!***********************************************************************

!***********************************************************************

      subroutine SIGSPL(s,y,dydx)
      USE TINoiseGeneric
      USE TINoiseDDD
      USE TINoiseRHSin
      USE TIPrecision
      IMPLICIT                        NONE

   ! Local variables.

INTEGER(4)                   :: k,khi,klo,i

REAL(DbKi)                   :: s,y1,y2,d1s,d2s,r2,r4,r6
REAL(DbKi)                   :: n1, n2,d1,d2,d1y1,d1y2,dsdt,arg
REAL(DbKi)                   :: y(6), dydx (6)

COMPLEX(DbKi)                :: pott, d1pott, green, HANK1, HANK0, HA0,HA1
COMPLEX(DbKi)                :: dkern1, dkern2, dkern3, dkern4, dkern5, dkern6
COMPLEX(DbKi)                :: gr, gr1, gr2, gr11, gr12, gr22, dumpa
COMPLEX(DbKi)                :: poadd,poads,d1azzz

!      call SPL_EX  (sworkat,na+1,s,khi,klo)
      khi = min(int(na+1),int(s+2))
!      khi = s + 2
      klo = khi-1
      call SPL_EX1 (sworkat,yc1at,d2yc1at,na+1,s,y1,d1y1,khi,klo)
      call SPL_EX1 (sworkat,yc2at,d2yc2at,na+1,s,y2,d1y2,khi,klo)
      CALL SPL_E1A (sworkat,potsat(1,ipath),d2potsat(1,ipath),na+1,s,pott,d1pott,khi,klo)

      dsdt = sqrt( d1y1**2 + d1y2**2)
      n1 =  d1y2 / dsdt
      n2 = -d1y1 / dsdt
      d1 = x1-y1
      d2 = x2-y2
      d1s = d1**2
      d2s = d2**2
      r2 = d1s + d2s
      arg = r2 * kwave2
!
!     HANK1(x**2) = HANK1 / x !!!!!!
!
!

      HA0  = HANK0(arg)
      HA1  = HANK1(arg)
      gr    = -0.25d0*imag*HA0
      gr1   = abb1*d1*HA1
      gr2   = abb1*d2*HA1
      dkern1 = ((n1*gr1 + n2*gr2) * pott) * dsdt
      dydx( 1) =  dble(dkern1)
      dydx( 2) =  dble(-imag*dkern1)

      if(lderiv) then
        gr11  = abb1*(HA0*d1s/r2+(1.0d0-2.0d0*d1s/r2)*HA1)
        gr12  = abb1*d1*d2/r2*(HA0-2.0d0*HA1)
        gr22  = abb1*(HA0*d2s/r2+(1.0d0-2.0d0*d2s/r2)*HA1)
        dkern2 = ((n1 * gr11 + n2 * gr12) * pott ) * dsdt
        dkern3 = ((n1 * gr12 + n2 * gr22) * pott ) * dsdt
        dydx( 3) =  dble(dkern2)
        dydx( 4) =  dble(-imag*dkern2)
        dydx( 5) =  dble(dkern3)
        dydx( 6) =  dble(-imag*dkern3)
      else
        dydx( 3) = 0.0
        dydx( 4) = 0.0
        dydx( 5) = 0.0
        dydx( 6) = 0.0
      endif

      return
      end

!***********************************************************************
!***********************************************************************

      subroutine SIGSUM(s,y,dydx)
      USE TINoiseGeneric
      USE TINoiseDDD
      USE TINoiseRHSin
      USE TIPrecision

IMPLICIT                        NONE

   ! Local variables.

INTEGER(4)                   :: k,khi,klo,i

REAL(DbKi)                   :: s,y1,y2,d1s,d2s,r2,r4,r6
REAL(DbKi)                   :: n1, n2,d1,d2,d1y1,d1y2,dsdt,arg
REAL(DbKi)                   :: y(10), dydx (10)

COMPLEX(DbKi)                :: pott, d1pott, green, HANK1, HANK0, HA0,HA1
COMPLEX(DbKi)                :: dkern1, dkern2, dkern3, dkern4, dkern5, dkern6
COMPLEX(DbKi)                :: gr, gr1, gr2, gr11, gr12, gr22, gr111, gr112
COMPLEX(DbKi)                :: gr122, gr222, dipstr, dumpa
COMPLEX(DbKi)                :: poadd,poads,d1azzz

!      call SPL_EX  (sworkat,na+1,s,khi,klo)
      khi = min(int(na+1),int(s+2))
      klo = khi-1
      call SPL_EX1 (sworkat,yc1at,d2yc1at,na+1,s,y1,d1y1,khi,klo)
      call SPL_EX1 (sworkat,yc2at,d2yc2at,na+1,s,y2,d1y2,khi,klo)
      CALL SPL_E1A (sworkat,potsum,d2potsum,na+1,s,pott,d1pott,khi,klo)

!      a = (sworkat(khi) - s)
!      b = (s - sworkat(klo))
!      a2 = a*a
!      b2 = b*b
!      a3 = a2*a
!      b3 = b2*b
!      a3ma = (a3-a) /6.
!      b3mb = (b3-b) /6.
!      a2p1 = (-3.0d0*a2 + 1.0d0)/6.0d0
!c      b2p1 = (+3.0d0*b2 - 1.0d0)/6.0d0
!      y1 = a * yc1at(klo) + b * yc1at(khi)
!     >     + (a3ma * d2yc1at(klo) + b3mb * d2yc1at(khi))
!      y2 = a * yc2at(klo) + b * yc2at(khi)
!     >     + (a3ma * d2yc2at(klo) + b3mb * d2yc2at(khi))
!      pott = a * potsat(klo,ipath) + b * potsat(khi,ipath)
!     >     + (a3ma * d2potsat(klo,ipath) +
!     +        b3mb * d2potsat(khi,ipath))
!      d1y1 = (-yc1at(klo) + yc1at(khi))
!     >       + (a2p1 * d2yc1at(klo)+b2p1 * d2yc1at(khi))
!      d1y2 = (-yc2at(klo) + yc2at(khi))
!     >       + (a2p1 * d2yc2at(klo)+b2p1 * d2yc2at(khi))

      dsdt = sqrt( d1y1**2 + d1y2**2)
      n1 =  d1y2 / dsdt
      n2 = -d1y1 / dsdt
      d1 = x1-y1
      d2 = x2-y2
      d1s = d1**2
      d2s = d2**2
      r2 = d1s + d2s
      arg = r2 * kwave2
!
!     HANK1(x**2) = HANK1 / x !!!!!!
!
!

      HA0  = HANK0(arg)
      HA1  = HANK1(arg)
      gr    = -0.25d0*imag*HA0
      gr1   = abb1*d1*HA1
      gr2   = abb1*d2*HA1
      dkern1 = ((n1*gr1 + n2*gr2) * pott ) * dsdt
      dydx( 1) =  dble(dkern1)
      dydx( 2) =  dble(-imag*dkern1)

      if(lderiv) then
        gr11  = abb1*(HA0*d1s/r2+(1.0d0-2.0d0*d1s/r2)*HA1)
        gr12  = abb1*d1*d2/r2*(HA0-2.0d0*HA1)
        gr22  = abb1*(HA0*d2s/r2+(1.0d0-2.0d0*d2s/r2)*HA1)
        gr111 = abb1*d1/r2*(-kwave2*d1s+8.0d0*d1s/r2-6.0d0)*HA1+abb1*d1/r2*(-4.0d0*d1s/r2+3.0d0)*HA0
        gr222 = abb1*d2/r2*(-kwave2*d2s+8.0d0*d2s/r2-6.0d0)*HA1+abb1*d2/r2*(-4.0d0*d2s/r2+3.0d0)*HA0
        gr112 = -kwave2 * gr2 - gr222
        gr122 = -kwave2 * gr1 - gr111

        dkern2 = ((n1 * gr11 + n2 * gr12) * pott ) * dsdt
        dkern3 = ((n1 * gr12 + n2 * gr22) * pott ) * dsdt
        dkern4 = ((n1 * gr111+ n2 *gr112) * pott ) * dsdt
        dkern5 = ((n1 * gr112+ n2 *gr122) * pott ) * dsdt
        dydx( 3) =  dble(dkern2)
        dydx( 4) =  dble(-imag*dkern2)
        dydx( 5) =  dble(dkern3)
        dydx( 6) =  dble(-imag*dkern3)
        dydx( 7) =  dble(dkern4)
        dydx( 8) =  dble(-imag*dkern4)
        dydx( 9) =  dble(dkern5)
        dydx(10) =  dble(-imag*dkern5)
      else
        dydx( 3) = 0.0
        dydx( 4) = 0.0
        dydx( 5) = 0.0
        dydx( 6) = 0.0
	dydx( 7) = 0.0
        dydx( 8) = 0.0
      	dydx( 9) = 0.0
      	dydx(10) = 0.0
      endif

      return
      end

!***********************************************************************

      subroutine CFIE_s(s,y,dydx)
      USE TINoiseGeneric
      USE TINoiseGeo
      USE TINoiseDDD
      USE TIPrecision

IMPLICIT                        NONE

   ! Local variables.

INTEGER(4)                   :: k,khi,klo,i

REAL(DbKi)                   :: s,y1,y2,d1s,d2s,r2,r4,r6,ppp,tang
REAL(DbKi)                   :: n1, n2,d1,d2,d1y1,d1y2,dsdt,arg
REAL(DbKi)                   :: y(3), dydx (3)

      call SPL_EX  (swork,n+1,s,khi,klo)
      call SPL_EX1 (swork,yc1,d2yc1,n+1,s,y1,d1y1,khi,klo)
      call SPL_EX1 (swork,yc2,d2yc2,n+1,s,y2,d1y2,khi,klo)
      CALL SPL_EX1 (swork, pots, d2pots, n+1, s, ppp, tang, khi, klo)

      n1 =  d1y2 !/ sqrt(d1y1**2+d1y2**2)
      n2 = -d1y1 !/ sqrt(d1y1**2+d1y2**2)
      d1 = x1-y1
      d2 = x2-y2
      d1s = d1**2
      d2s = d2**2
      r2 = d1s + d2s
      r4 = r2 ** 2

      dydx(1) = (n1*d1+n2*d2)/r2* ppp * pi2i
      dydx(2) = (n1*(d2s-d1s)-2.0d0*n2*d2*d1)/r4 * ppp * pi2i
      dydx(3) = (n2*(d1s-d2s)-2.0d0*n1*d2*d1)/r4 * ppp * pi2i

      return
      end

!***********************************************************************
REAL(DbKi) function ansatz(r)
      USE TIPrecision
      IMPLICIT                        NONE
      REAL(DbKi)                   :: r
      ansatz = r**2/4.0d0 + r**3/9.0d0
      return
      end function
!***********************************************************************
REAL(DbKi)  function ansatd(r)
      USE TIPrecision
      IMPLICIT                        NONE
      REAL(DbKi)                   :: r
      ansatd = 1.0d0 + r
      return
      end function
!***********************************************************************
REAL(DbKi) function dsatzo(r)
      USE TINoiseGeneric
      USE TIPrecision
      IMPLICIT                        NONE
      REAL(DbKi)                   :: r
      dsatzo = 0.5+r/3.0d0
      return
      end function

max_stars_repo_name: CarloSucameli/Openfast-noise

max_issues_repo_name: CarloSucameli/Openfast-noise

max_forks_repo_name: CarloSucameli/Openfast-noise


text:
! solve b=Ax using parallel log-ent mart.
! Matlab by Joshua Semeter
! Fortran by Michael Hirsch

program test_logmart

use, intrinsic:: iso_fortran_env, only: stderr=>error_unit
use art, only: logmart
use random_utils, only: init_random_seed, randn, wp

implicit none

integer, parameter :: N=4, M=4
real(wp), dimension(N,M) :: A1, A2
real(wp) :: x(M)
logical :: ok=.true., add_noise = .true.

call init_random_seed()

A1 = reshape([5,0,0,0, &
              0,5,0,0, &
              0,0,5,0, &
              0,0,0,5], shape(A1), order=[2,1])

A2 = reshape([0,1,2,3, &
              1,0,1,2, &
              2,1,0,1, &
              3,2,1,0], shape(A2), order=[2,1])

x = [1._wp, 3._wp, 0.5_wp, 2._wp]

if (.not. run_test(A1, x, 20,add_noise)) then
  ok = .false.
  write(stderr,*) 'failed on identity test'
endif

if (.not. run_test(A2, x, 2000, add_noise)) then
  ok = .false.
  write(stderr,*) 'failed on Fiedler test'
endif

if (.not. ok) error stop

print *, 'OK: logmart'

contains

logical function run_test(A, x, max_iter, add_noise)

real(wp), intent(in) :: A(:,:), x(:)
logical, intent(in) :: add_noise
integer, intent(in) :: max_iter

real(wp), parameter :: errtol=0.05_wp
real(wp), dimension(size(A,1)) :: noise, x_est, b

run_test = .true.

block
integer :: i
print *, 'A ='
do i = 1,size(A,1)
  print '(3F10.3)', A(i,:)
enddo
end block

! ---- noisy observation
b = matmul(A,x)

if (add_noise) then
  call randn(noise)
  noise = 0.01_wp * noise
  print '(/,A,3F10.3)', 'noise',noise
  b = b + noise
endif

! ---- inversion
call logmart(A,b, max_iter=max_iter, x=x_est)

! --- check estimate
if (any(abs(x_est-x) > errtol*maxval(x))) then
  print *,x_est
  write (stderr,*) 'larger than',errtol*100,' % error'
  run_test = .false.
endif

end function run_test

end program

max_stars_repo_name: scienceopen/pyAIRtools

max_issues_repo_name: scienceopen/pyAIRtools

max_forks_repo_name: scienceopen/pyAIRtools


text: (cond true true)

max_stars_repo_name: cappuccinotm/flangc

max_issues_repo_name: cappuccinotm/flangc

max_forks_repo_name: cappuccinotm/flangc