src/tensors/dbcsr_tensor_test.F

192       INTEGER                                     :: io_unit
202       io_unit = 0
203       IF (mynode .EQ. 0) io_unit = unit_nr
209       IF (io_unit > 0) THEN
210          WRITE (io_unit, *)
211          WRITE (io_unit, '(A)') repeat("-", 80)
212          WRITE (io_unit, '(A,1X,I1)') "Testing matrix representations of tensor rank", ndims
213          WRITE (io_unit, '(A)') repeat("-", 80)
214          WRITE (io_unit, '(A)') "Block sizes:"
218             WRITE (io_unit, '(T4,A,1X,I1,A,1X)', advance='no') 'Dim', ${dim}$, ':'
220                WRITE (io_unit, '(I2,1X)', advance='no') blk_size_${dim}$ (iblk)
222             WRITE (io_unit, *)
226          WRITE (io_unit, '(A)') "Non-zero blocks:"
230                WRITE (io_unit, '(T4,A, I3, A, ${ndim}$I3, 1X, A)') &
236          WRITE (io_unit, *)
237          WRITE (io_unit, '(A,1X)', advance='no') "Reference map:"
238          WRITE (io_unit, '(A1,1X)', advance='no') "("
240             WRITE (io_unit, '(I1,1X)', advance='no') map1_ref(imap)
242          WRITE (io_unit, '(A1,1X)', advance='no') "|"
244             WRITE (io_unit, '(I1,1X)', advance='no') map2_ref(imap)
246          WRITE (io_unit, '(A1)') ")"
273             IF (io_unit > 0) THEN
274                WRITE (io_unit, *)
275                WRITE (io_unit, '(A,A,1X)', advance='no') TRIM(tensor_name), ':'
276                WRITE (io_unit, '(A1,1X)', advance='no') "("
278                   WRITE (io_unit, '(I1,1X)', advance='no') map1(imap)
280                WRITE (io_unit, '(A1,1X)', advance='no') "|"
282                   WRITE (io_unit, '(I1,1X)', advance='no') map2(imap)
284                WRITE (io_unit, '(A1)') ")"
286                WRITE (io_unit, '(T4,A)') "Reference distribution:"
289                   WRITE (io_unit, '(T7,A,1X)', advance='no') "Dist vec ${dim}$:"
291                      WRITE (io_unit, '(I2,1X)', advance='no') dist2_${dim}$ (idist)
293                   WRITE (io_unit, *)
297                WRITE (io_unit, '(T4,A)') "Test distribution:"
300                   WRITE (io_unit, '(T7,A,1X)', advance='no') "Dist vec ${dim}$:"
302                      WRITE (io_unit, '(I2,1X)', advance='no') dist1_${dim}$ (idist)
304                   WRITE (io_unit, *)
318             IF (verbose) CALL dbcsr_t_write_blocks(tensor2, io_unit, unit_nr)
329             IF (verbose) CALL dbcsr_t_write_blocks(tensor1, io_unit, unit_nr)
334                IF (io_unit > 0) WRITE (io_unit, '(A,1X,A)') TRIM(tensor_name), 'Test failed!'
337                IF (io_unit > 0) WRITE (io_unit, '(A,1X,A)') TRIM(tensor_name), 'Test passed!'
563       INTEGER                              :: io_unit, mynode, numnodes, mp_comm  local
594       io_unit = -1
595       IF (mynode .EQ. 0) io_unit = unit_nr
601       IF (io_unit > 0) THEN
602          WRITE (io_unit, *)
603          WRITE (io_unit, '(A)') repeat("-", 80)
604          WRITE (io_unit, '(A,1X,A,1X,A,1X,A,1X,A,1X,A)') "Testing tensor contraction", &
606          WRITE (io_unit, '(A)') repeat("-", 80)
610          IF (io_unit > 0) THEN
611             WRITE (io_unit, "(A, E9.2)") "checksum ", TRIM(tensor_1%name), cs_1
612             WRITE (io_unit, "(A, E9.2)") "checksum ", TRIM(tensor_2%name), cs_2
613             WRITE (io_unit, "(A, E9.2)") "checksum ", TRIM(tensor_3%name), cs_3
618          CALL dbcsr_t_write_block_indices(tensor_1, io_unit, unit_nr)
619          CALL dbcsr_t_write_blocks(tensor_1, io_unit, unit_nr, write_int)
635                             unit_nr=io_unit, log_verbose=log_verbose)
640          IF (io_unit > 0) THEN
641             WRITE (io_unit, "(A, E9.2)") "checksum ", TRIM(tensor_3%name), cs_3
742          IF (io_unit > 0) WRITE (io_unit, *) 'Test failed!', eql_diff
745          IF (io_unit > 0) WRITE (io_unit, *) 'Test passed!', eql_diff